Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,038,378 C→T 5.5% R213R (CGC→CGT 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,3780CT 5.5% 286.2 / 4.2 90R213R (CGC→CGTappCcytochrome bd‑II ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (37/48);  new base T (2/3);  total (39/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

GCGTATGGAGATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGG  >  NC_000913/1038231‑1038526
                                                                                                                                                   |                                                                                                                                                    
gCGTATGGAGATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTTATGGCGGGCTACGTGACCGCGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCaa                                                                                                                                                    >  2:1240839/1‑150 (MQ=255)
  gTATGGAGATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATg                                                                                                                                                  <  1:921662/150‑1 (MQ=255)
   tATGGAGATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGt                                                                                                                                                 >  1:21869/1‑150 (MQ=255)
      ggAGATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGc                                                                                                                                              <  2:1577393/150‑1 (MQ=255)
      ggAGATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGc                                                                                                                                              <  1:63266/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:1593704/150‑1 (MQ=255)
                              ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                      <  1:684368/150‑1 (MQ=255)
                              ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                      <  2:1020467/150‑1 (MQ=255)
                              ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                      >  2:1457188/1‑150 (MQ=255)
                              ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                      <  1:1295918/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                                                                                                                                <  1:655769/128‑1 (MQ=255)
                                          ggTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTg                                                                                                                                >  2:655769/1‑128 (MQ=255)
                                            tCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTActc                                                                                                        <  1:665888/150‑1 (MQ=255)
                                            tCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTActc                                                                                                        <  1:868755/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                                                                                               >  2:422114/1‑150 (MQ=255)
                                                     tGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTa                                                                                               >  1:1078548/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                                                             <  1:315832/150‑1 (MQ=255)
                                                                                       gACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTa                                                             >  2:1074083/1‑150 (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                                       >  2:1405325/1‑150 (MQ=255)
                                                                                                             gCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCg                                       >  1:656888/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)
                                                                                                                                                   |                                                                                                                                                    
GCGTATGGAGATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGG  >  NC_000913/1038231‑1038526

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: