Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,038,383 T→C 5.4% V215A (GTC→GCC) ‡ appC → cytochrome bd‑II ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,038,3830TC 5.4% 287.8 / 5.0 93V215A (GTC→GCC) ‡appCcytochrome bd‑II ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base T (36/52);  new base C (2/3);  total (38/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

TATGGAGATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTG  >  NC_000913/1038234‑1038532
                                                                                                                                                     |                                                                                                                                                     
tATGGAGATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGt                                                                                                                                                       >  1:21869/1‑150 (MQ=255)
   ggAGATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGc                                                                                                                                                    <  1:63266/150‑1 (MQ=255)
   ggAGATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGc                                                                                                                                                    <  2:1577393/150‑1 (MQ=255)
       gtGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTa                                                                                                                                                <  1:125318/149‑1 (MQ=255)
           ccAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCt                                                                                                                                            <  1:938101/150‑1 (MQ=255)
             aGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGTAATGCGACATTGCGCTcc                                                                                                                                          <  2:13382/150‑2 (MQ=255)
             aGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGTAATGCGACATTGCGCTcc                                                                                                                                          <  2:40358/150‑2 (MQ=255)
               cTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGtt                                                                                                                                        <  2:456242/150‑1 (MQ=255)
                ttCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGttt                                                                                                                                       <  2:224810/150‑1 (MQ=255)
                  cAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGc                                                                                                                                     <  2:1342050/150‑1 (MQ=255)
                  cAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGc                                                                                                                                     <  2:1355005/150‑1 (MQ=255)
                   aGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGcc                                                                                                                                    >  2:868755/1‑150 (MQ=255)
                    gCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCa                                                                                                                                   >  1:1246597/1‑150 (MQ=255)
                     cGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCAt                                                                                                                                  >  2:821202/1‑150 (MQ=255)
                      gAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATc                                                                                                                                 <  1:1225928/150‑1 (MQ=255)
                       aGCTGGTCTTTAATCCTGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCTGTGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGTAGCGCAATGTCGCATTACGCTCGTTTGCCATCg                                                                                                                                >  2:1640823/1‑150 (MQ=255)
                       aGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCg                                                                                                                                <  2:792103/150‑1 (MQ=255)
                        gCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCgg                                                                                                                               >  2:695181/1‑150 (MQ=255)
                           ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                            <  2:1020467/150‑1 (MQ=255)
                           ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                            <  1:684368/150‑1 (MQ=255)
                           ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                            <  1:1295918/150‑1 (MQ=255)
                           ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                            >  2:1457188/1‑150 (MQ=255)
                           ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                            <  2:1593704/150‑1 (MQ=255)
                             tCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCg                                                                                                                          <  2:1641520/150‑1 (MQ=255)
                              cTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGt                                                                                                                         >  2:1206065/1‑150 (MQ=255)
                              cTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGt                                                                                                                         >  1:159521/1‑150 (MQ=255)
                               tttAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTc                                                                                                                        >  2:1252704/1‑150 (MQ=255)
                                   aTCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTTTGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCg                                                                                                                    >  1:993782/1‑150 (MQ=255)
                                    tCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCgg                                                                                                                   >  2:1395384/1‑150 (MQ=255)
                                     ccGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGt                                                                                                                  <  1:113160/150‑1 (MQ=255)
                                       ggTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTg                                                                                                                                      >  2:655769/1‑128 (MQ=255)
                                       ggTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTg                                                                                                                                      <  1:655769/128‑1 (MQ=255)
                                         tCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTActc                                                                                                              <  1:868755/150‑1 (MQ=255)
                                         tCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTActc                                                                                                              <  1:665888/150‑1 (MQ=255)
                                          cAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTActct                                                                                                             <  1:572015/150‑1 (MQ=255)
                                            gCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTgg                                                                                                           >  1:792758/1‑150 (MQ=255)
                                                ggTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCGTCGGTTCCGTCTTCGGTACTCTGGCGAt                                                                                                       >  2:599365/1‑150 (MQ=255)
                                                 gTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGAtt                                                                                                      >  2:1368548/1‑150 (MQ=255)
                                                  tGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTa                                                                                                     >  1:1078548/1‑150 (MQ=255)
                                                  tGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTa                                                                                                     >  2:422114/1‑150 (MQ=255)
                                                     aaTTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCg                                                                                                  <  1:1589361/150‑1 (MQ=255)
                                                        ttGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTa                                                                                               <  1:1206065/150‑1 (MQ=255)
                                                          gtgCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTAcc                                                                                             <  2:1084196/150‑1 (MQ=255)
                                                           tgCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTTTGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTAccc                                                                                            <  2:993782/150‑1 (MQ=255)
                                                           tgCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTAccc                                                                                            >  1:40318/1‑150 (MQ=255)
                                                            gCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCt                                                                                           >  1:993220/1‑150 (MQ=255)
                                                               cacCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCa                                                                                        <  2:1310314/150‑1 (MQ=255)
                                                                 ccGTAATGGCGGGCTACGTTACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTTCAAc                                                                                      <  2:157392/150‑1 (MQ=255)
                                                                 ccGTAATGGCGGGCTACGTGACCGGGGCAATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAAc                                                                                      >  2:748492/1‑150 (MQ=255)
                                                                     aaTGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCgg                                                                                  <  2:1364385/150‑1 (MQ=255)
                                                                      aTGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCAATTATCGGTACCCTGCAACTCGga                                                                                 <  1:528046/150‑1 (MQ=255)
                                                                         gCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGAc                                                                              <  1:695181/150‑1 (MQ=255)
                                                                          cGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACa                                                                             >  1:427810/1‑150 (MQ=255)
                                                                               tACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCt                                                                        <  1:266732/150‑1 (MQ=255)
                                                                                   ggACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAATTCGGAGACAGTTCTGCGt                                                                    <  2:1554973/149‑1 (MQ=255)
                                                                                    gACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTa                                                                   >  2:1074083/1‑150 (MQ=255)
                                                                                    gACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTa                                                                   <  1:315832/150‑1 (MQ=255)
                                                                                     cccGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTAt                                                                  <  1:203160/149‑1 (MQ=255)
                                                                                     aCCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTAt                                                                  >  1:311882/1‑150 (MQ=255)
                                                                                      ccttggCCATGTTTATTATGGCGATCAGCGCCTTGTATTTACTTCGCGTACGGTAGCGCAATGTCTCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATg                                                                 <  1:634537/146‑1 (MQ=255)
                                                                                      ccGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATg                                                                 >  1:763351/1‑150 (MQ=255)
                                                                                       cGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGa                                                                <  2:159521/150‑1 (MQ=255)
                                                                                        ggggCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGaa                                                               >  1:686314/1‑150 (MQ=255)
                                                                                         gggCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAg                                                              <  1:1046817/150‑1 (MQ=255)
                                                                                          ggCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGt                                                             <  1:1608060/150‑1 (MQ=255)
                                                                                              aTGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTcgcg                                                         <  1:1240839/150‑1 (MQ=255)
                                                                                               tGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCgcgc                                                        >  1:1125935/1‑150 (MQ=255)
                                                                                                  ttattaTGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAg                                                     >  2:1184293/1‑150 (MQ=255)
                                                                                                  ttattaTGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGt                                                    <  2:152599/150‑1 (MQ=255)
                                                                                                   tattatGGCGATCAGCGCCTGGTATTTACTGAGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGt                                                    >  2:430849/1‑150 (MQ=255)
                                                                                                     ttatGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTAc                                                  >  1:982495/1‑150 (MQ=255)
                                                                                                        tGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAAc                                               <  1:1074083/150‑1 (MQ=255)
                                                                                                        tGGCGATCAGAGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAAc                                               >  2:278926/1‑150 (MQ=255)
                                                                                                          gCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCg                                             >  1:656888/1‑150 (MQ=255)
                                                                                                          gCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCg                                             >  2:1405325/1‑150 (MQ=255)
                                                                                                           cGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCgg                                            >  2:1046762/1‑150 (MQ=255)
                                                                                                                     cTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTggc                                  >  2:286592/1‑150 (MQ=255)
                                                                                                                         tATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGgcggcg                              <  2:656888/150‑1 (MQ=255)
                                                                                                                              aCTGAGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGa                         <  1:430849/150‑1 (MQ=255)
                                                                                                                                tGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAg                       >  2:857334/1‑150 (MQ=255)
                                                                                                                                tGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAg                       <  1:794219/150‑1 (MQ=255)
                                                                                                                                      gACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGt                 <  2:1311719/150‑1 (MQ=255)
                                                                                                                                       aCGGTAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTg                <  2:263222/150‑1 (MQ=255)
                                                                                                                                        cGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTgg               <  2:993220/150‑1 (MQ=255)
                                                                                                                                        aacgAGCGTAATGCGACATTGCGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATg                           >  1:1334874/4‑138 (MQ=255)
                                                                                                                                        aacgAGCGTAATGCGACATTGCGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATg                           <  2:1334874/135‑1 (MQ=255)
                                                                                                                                          ggAGCGCAATGTCGCATTACGCTCGTTTGCCGTCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCa             <  1:599365/150‑1 (MQ=255)
                                                                                                                                          cgAGCGTAATGCGACATTGCGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCa             >  1:376302/2‑150 (MQ=255)
                                                                                                                                             gcgcAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACGGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAAc          <  1:977270/150‑1 (MQ=255)
                                                                                                                                               gcAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTAGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACgg        >  1:786324/1‑150 (MQ=255)
                                                                                                                                                aaaTGTCGCATTACGCTCGTTTTCCATCGTTTCCGTCTTCGGTACTCTTGCTATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGTAGTCGCGCAAGTACAACCGGTAAATCTGGCGGAGATGGAAGGGGAGTGGCAAACGGa       <  2:1246597/149‑1 (MQ=255)
                                                                                                                                                   tGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacc    <  2:792758/150‑1 (MQ=255)
                                                                                                                                                   tGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacc    <  2:1125935/150‑1 (MQ=255)
                                                                                                                                                   tGTCGCATTACGCTCGTTTGCCATCGGTTCAGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacc    <  1:1469199/150‑1 (MQ=255)
                                                                                                                                                     tCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacctg  <  1:748492/150‑1 (MQ=255)
                                                                                                                                                     tCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacctg  >  1:525461/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATGGAGATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTG  >  NC_000913/1038234‑1038532

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: