Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,453,151 C→T L141L (CTG→TTG)  unknown → Bacterial non‑heme ferritin (EC 1.16.3.2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,453,1510CT100.0% 195.7 / NA 68L141L (CTG→TTG) unknownBacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (43/25);  total (43/25)

CATTGGCTGTCATTACTGCAATTGATGAACCGCCACACCAATTCGCTTCGTTAAAACAAGTACTTGAATTAACATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAACTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCGGTGGTTGACCTGTTCCACTGGCTTCATAATAAAACGTGCGATAAAAAAAT  >  1/1452927‑1453389
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
cATTGGCTGTCATTACTGCAATTGATGAACCGCCACACCAATTCGCTTCGTTAAAACAAGTACTTGAATTAACATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGATACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGaaaa                                                                                                                                                                                                                      >  1:431927/1‑251 (MQ=255)
  ttgccTGTCATTACTGAAATTGATGAACCGCCACACCAATTCGCTTCGTTAAAACAAGTACTTGAATTAACATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAgg                                                                                                                                                                                                                    <  1:126791/247‑1 (MQ=255)
         tCATTACTGCAATTGATGAACCGCCACACCAATTCGCTTCGTTAAAACAAGTACTTGAATTAACATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTc                                                                                                                                                                                                             <  1:818486/251‑1 (MQ=255)
         tCATTACTGAAATTGATGAACCGCCACACCAATTCGCTTCGTTAAAACAAGTACTTGAATTAACATACGAACATGAAAAATGGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTc                                                                                                                                                                                                             <  1:452220/251‑1 (MQ=255)
            ttACTGCAATTGATGAACCGCCACACCAATTCGCTTCGTTAAAACAAGTACTTGAATTAACATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCtt                                                                                                                                                                                                          >  1:472796/1‑251 (MQ=255)
                      tgatgaACCGCCACACCAATTCGCTTCGTTAAAACAAGTACTTGAATTAACATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAg                                                                                                                                                                                                <  1:810329/251‑1 (MQ=255)
                       gatgaACCGCCACACCAATTCGCTTCGTTAAAACAAGTACTTGAATTAACATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGa                                                                                                                                                                                               >  1:102757/1‑251 (MQ=255)
                                   cacCAATTCGCTTCGTTAAAACAAGTACTTGAATTAACATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTcc                                                                                                                                                                                                            <  1:424766/226‑1 (MQ=255)
                                                           gTACTTGAATTAACATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGACTCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGttt                                                                                                                                                           >  1:426406/1‑251 (MQ=255)
                                                           gTACTTGAATTAACATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGttt                                                                                                                                                           >  1:83611/1‑251 (MQ=255)
                                                               ttGAATTAACATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTtaaat                                                                                                                                                       <  1:62916/251‑1 (MQ=255)
                                                                      aaCATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTc                                                                                                                                                                                                             >  1:630959/1‑190 (MQ=255)
                                                                        cATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTa                                                                                                                                              >  1:800420/1‑251 (MQ=255)
                                                                           aCGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGt                                                                                                                                           <  1:765711/251‑1 (MQ=255)
                                                                                cATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACAt                                                                                                                                      >  1:324304/1‑251 (MQ=255)
                                                                                cATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACAt                                                                                                                                      >  1:767376/1‑251 (MQ=255)
                                                                                cATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTATTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGAt                                                                                                                                                                                              >  1:684380/1‑195 (MQ=255)
                                                                                      aaaTTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGAt                                                                                                                                                                                              >  1:265647/1‑189 (MQ=255)
                                                                                      aaaTTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGAttt                                                                                                                                <  1:460642/251‑1 (MQ=255)
                                                                                        aTTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAg                                                                                                                              >  1:224036/1‑251 (MQ=255)
                                                                                          tGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGcc                                                                                                                            <  1:395489/251‑1 (MQ=255)
                                                                                                    cAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTttcttcttc                                                                                                                  <  1:559018/251‑1 (MQ=255)
                                                                                                       aaTCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCa                                                                                                               >  1:876061/1‑251 (MQ=255)
                                                                                                         tCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCAtt                                                                                                             >  1:272938/1‑251 (MQ=255)
                                                                                                         tCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCAtt                                                                                                             >  1:765217/1‑251 (MQ=255)
                                                                                                                        gtAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGtt                                                                                                                                                            <  1:494397/189‑1 (MQ=255)
                                                                                                                        gtAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAAtttt                                                                                              >  1:410668/1‑251 (MQ=255)
                                                                                                                              cATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTAtt                                                                                          >  1:12879/1‑249 (MQ=255)
                                                                                                                                tttGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAg                                                                                                                                                                                                >  1:634792/1‑145 (MQ=255)
                                                                                                                                tttGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCa                                                                                      >  1:596680/1‑251 (MQ=255)
                                                                                                                                   gAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACt                                                                                                                                                                     >  1:564521/1‑169 (MQ=255)
                                                                                                                                      gagaAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAgg                                                                                                                                                                                                                    <  1:156313/119‑1 (MQ=255)
                                                                                                                                      gagaAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGt                                                                                >  1:801525/1‑251 (MQ=255)
                                                                                                                                        gaAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGcc                                                                                                                            >  1:689739/1‑205 (MQ=255)
                                                                                                                                        gaAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAAttt                                                                                               <  1:471609/234‑1 (MQ=255)
                                                                                                                                           ggACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAg                                                                                                                                                                                                >  1:310430/1‑134 (MQ=255)
                                                                                                                                           ggACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCgg                                                                            <  1:715874/250‑1 (MQ=255)
                                                                                                                                           ggACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGgt                                                                           >  1:789780/1‑251 (MQ=255)
                                                                                                                                           ggACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGgt                                                                           >  1:676480/1‑251 (MQ=255)
                                                                                                                                            gACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCgg                                                                            >  1:257531/1‑249 (MQ=255)
                                                                                                                                              ctactCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTtaaat                                                                                                                                                       <  1:85106/172‑1 (MQ=255)
                                                                                                                                                actCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGAttt                                                                                                                                <  1:178646/193‑1 (MQ=255)
                                                                                                                                                     gCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAAtc                                                                                                                                                                                                                                          >  1:547659/1‑82 (MQ=255)
                                                                                                                                                     gCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGgc                                                                                                                                                                                        <  1:125776/132‑1 (MQ=255)
                                                                                                                                                     gCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCa                                                                 >  1:50588/1‑251 (MQ=255)
                                                                                                                                                            aCTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGttt                                                          >  1:168805/1‑251 (MQ=255)
                                                                                                                                                                 ttACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGtt                                                                                                                                                            >  1:259353/1‑148 (MQ=255)
                                                                                                                                                                    cAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCggt                                                  >  1:726615/1‑251 (MQ=255)
                                                                                                                                                                      aTGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGttt                                                                                                                                                           <  1:63176/144‑1 (MQ=255)
                                                                                                                                                                      aTGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGc                                                                                                                                                <  1:795184/155‑1 (MQ=255)
                                                                                                                                                                      aTGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCggtgg                                                >  1:56596/1‑251 (MQ=255)
                                                                                                                                                                            cGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCGGTGGTTGAcc                                          <  1:462852/251‑1 (MQ=255)
                                                                                                                                                                              ttGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAAc                                                                                                                                                                      >  1:146687/1‑125 (MQ=255)
                                                                                                                                                                              ttGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTtaaata                                                                                                                                                      <  1:42056/141‑1 (MQ=255)
                                                                                                                                                                              ttGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTc                                                                                                                                          >  1:381400/1‑153 (MQ=255)
                                                                                                                                                                              ttGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTg                                                                               >  1:141624/1‑212 (MQ=255)
                                                                                                                                                                              ttGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCGGTGGTTGAc                                           >  1:613390/1‑248 (MQ=255)
                                                                                                                                                                              ttGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCGGTGGTTGAc                                           >  1:373817/1‑248 (MQ=255)
                                                                                                                                                                              ttGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCGGTGGTTGAc                                           >  1:15738/1‑248 (MQ=255)
                                                                                                                                                                               acaccgcactttctgagaataaAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACt                                                                                                     >  1:265622/16‑189 (MQ=255)
                                                                                                                                                                                   gaaCAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCGGTGGTTGACCTGTTCCa                                   >  1:351329/1‑251 (MQ=255)
                                                                                                                                                                                    aaCAGCACGAAGAAGAAAAATTGTTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCGGTGGTTGACCTGTTcc                                    <  1:740942/249‑1 (MQ=255)
                                                                                                                                                                                    aaCAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCGGTGGTTGACCTGTTCCAc                                  <  1:3771/251‑1 (MQ=255)
                                                                                                                                                                                         cACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCGGTGGTTGACCTGTTCCACTGGc                              <  1:720431/250‑1 (MQ=255)
                                                                                                                                                                                                    aaaaTTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCGGTGGTTGACCTGTTCCACTGGCTTCATAATaaaa                  <  1:577762/251‑1 (MQ=255)
                                                                                                                                                                                                     aaaTTATTCAGCGGAATTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCttg                                                                                                                                                                                                         >  1:687540/1‑67 (MQ=255)
                                                                                                                                                                                                                    aTTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCGGTGGTTGACCTGTTCCACTGGCTTCATAATAAAACGTGCGATAAAAAAAt  >  1:465249/1‑251 (MQ=255)
                                                                                                                                                                                                                    aTTTTGGATAAATTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCGGTGGTTGACCTGTTCCACTGGCTTCATAATAAAACGTGCGATAAAAAAAt  >  1:571268/1‑251 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
CATTGGCTGTCATTACTGCAATTGATGAACCGCCACACCAATTCGCTTCGTTAAAACAAGTACTTGAATTAACATACGAACATGAAAAATTGATTACTGCCAAAATCAATGAATTGGTTGGTAGAACATTTGAAGAGAAGGACTACTCTGCTTTTAACTTCTTACAATGGTACGTTGAAGAACAGCACGAAGAAGAAAAATTATTCAGCGGAATTTTGGATAAACTGAATCTTCTTGGCGAAGATGGAAAAGGTTTATTCCTTGTTGATAAAGATTTAGGCGCTTTAGCAGAAAGCAACTAGTCCGGTTTAAATAAAAGGCTAAGTCACATTGATTTAGCCTTTCTTTTTCTCATTCCAAAACTTAATTTTTATTCTCAGAAAGTGCGGTGTTTTTTTCAAGTGTTTTTGCCGGTGGTTGACCTGTTCCACTGGCTTCATAATAAAACGTGCGATAAAAAAAT  >  1/1452927‑1453389

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑