Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,402,074 A→C 6.4% L375R (CTG→CGG) ‡ PP_3881 ← terminase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,402,0740AC 6.4% 278.1 / 9.4 93L375R (CTG→CGG) ‡PP_3881terminase large subunit
Reads supporting (aligned to +/- strand):  ref base A (44/43);  new base C (3/3);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01

TCCACGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCTTCCTTG  >  NC_002947/4401925‑4402220
                                                                                                                                                     |                                                                                                                                                  
tCCACGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACa                                                                                                                                                    >  1:1522315/1‑150 (MQ=255)
tCCACGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACa                                                                                                                                                    >  2:1995433/1‑150 (MQ=255)
 ccACGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAg                                                                                                                                                   <  1:1452746/150‑1 (MQ=255)
 ccACGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAg                                                                                                                                                   <  2:1308424/150‑1 (MQ=255)
   aCGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGt                                                                                                                                                 >  1:640256/1‑150 (MQ=255)
       cgcgTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAg                                                                                                                                             <  2:1853988/150‑1 (MQ=255)
        gcgTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGt                                                                                                                                            <  1:2107710/150‑1 (MQ=255)
          gTATGGCACTTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGGGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTcc                                                                                                                                          <  1:973272/150‑1 (MQ=255)
                   cTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCg                                                                                                                                         >  1:1699050/1‑142 (MQ=255)
                   cTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCg                                                                                                                                         <  2:1699050/142‑1 (MQ=255)
                    ttGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTg                                                                                                                                >  2:2259712/1‑150 (MQ=255)
                    ttGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTg                                                                                                                                >  2:537735/1‑150 (MQ=255)
                    ttGTCGCGGGCGTCGCGCTCGTCGAGCTCGTGGTCGGGTACCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTg                                                                                                                                <  1:752768/150‑1 (MQ=255)
                     tGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTgg                                                                                                                               <  1:537735/150‑1 (MQ=255)
                       tCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGa                                                                                                                             <  2:2022660/150‑1 (MQ=255)
                           gggCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGAt                                                                                                                         <  2:1522315/150‑1 (MQ=255)
                            ggCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGAtt                                                                                                                        <  2:1896422/150‑1 (MQ=255)
                                    gcTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACa                                                                                                                >  2:1946429/1‑150 (MQ=255)
                                      ttGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGc                                                                                                              <  2:366360/150‑1 (MQ=255)
                                        gTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGa                                                                                                            >  1:355983/1‑150 (MQ=255)
                                         tCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGAc                                                                                                           <  2:1924468/150‑1 (MQ=255)
                                           gAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAg                                                                                                                                                   <  1:1112391/108‑1 (MQ=255)
                                           gAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAg                                                                                                                                                   >  2:1112391/1‑108 (MQ=255)
                                                       gggTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGa                                                                                                                             >  2:472506/1‑118 (MQ=255)
                                                       gggTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGa                                                                                                                             <  1:472506/118‑1 (MQ=255)
                                                            tCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAg                                                                                        >  1:1356318/1‑150 (MQ=255)
                                                               aGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGAc                                                                                                           <  1:103887/128‑1 (MQ=255)
                                                               aGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGAc                                                                                                           >  2:103887/1‑128 (MQ=255)
                                                               aGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCt                                                                                     >  2:595467/1‑150 (MQ=255)
                                                                gAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTg                                                                                    >  1:1581501/1‑150 (MQ=255)
                                                                 aaaTACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTgg                                                                                   <  1:69987/150‑1 (MQ=255)
                                                                          aaCAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACa                                                                          >  2:1589499/1‑150 (MQ=255)
                                                                            cAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATc                                                                        >  2:493776/1‑150 (MQ=255)
                                                                             aGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATcc                                                                       >  2:2007056/1‑150 (MQ=255)
                                                                               ccgccAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCa                                                                     >  2:1502633/1‑150 (MQ=255)
                                                                                 gccAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAgg                                                                   >  2:417344/1‑150 (MQ=255)
                                                                                  ccAGTGCGGATCCTCGCATGTCGGGTAGAACAGCAGGACAAATGCTGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCACCGCTGGGACGATTGCGCAACAGCGACATCGGTACGAGTCGTTCAGCCTGGCTCCAGACATCCCAGga                                                                  <  2:1483717/150‑1 (MQ=255)
                                                                                   cAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCt                                                                                     <  2:765450/130‑1 (MQ=255)
                                                                                   cAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCt                                                                                     >  1:765450/1‑130 (MQ=255)
                                                                                    aGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGgaga                                                                <  1:391774/150‑1 (MQ=255)
                                                                                    aGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGgaga                                                                <  2:640256/150‑1 (MQ=255)
                                                                                    aGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGgaga                                                                <  1:973856/150‑1 (MQ=255)
                                                                                     gTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAg                                                                                                                                                   >  1:211182/1‑66 (MQ=255)
                                                                                     gTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAg                                                                                                                                                   <  2:211182/66‑1 (MQ=255)
                                                                                          gATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACg                                                          >  2:476973/1‑150 (MQ=255)
                                                                                           aTCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGAc                                                                                                                            >  1:1861471/1‑83 (MQ=255)
                                                                                           aTCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGAc                                                                                                                            <  2:1861471/83‑1 (MQ=255)
                                                                                              cTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGa                                                      <  1:1060428/150‑1 (MQ=255)
                                                                                               tCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGAc                                                     >  1:930974/1‑150 (MQ=255)
                                                                                               tCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGAc                                                     >  1:1481223/1‑150 (MQ=255)
                                                                                                     gTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCa                                               <  1:1589499/150‑1 (MQ=255)
                                                                                                        gggTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTAcgcg                                                                                                >  1:888847/1‑98 (MQ=255)
                                                                                                        gggTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTAcgcg                                                                                                <  2:888847/98‑1 (MQ=255)
                                                                                                          gTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGt                                          >  1:1710925/1‑150 (MQ=255)
                                                                                                           tAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTc                                         <  1:1946429/150‑1 (MQ=255)
                                                                                                             gAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCa                                       <  1:2259712/150‑1 (MQ=255)
                                                                                                             gAACAGCAGGACAAATGATGTCAGGTCCGTCGGACTGGACCTGTCCAGtac                                                                                                                                          <  2:503683/51‑3 (MQ=21)
                                                                                                             gAACAGCAGGACAAATGATGTCAGGTCCGTCGGACTGGACCTGTCCAGtac                                                                                                                                          >  1:503683/1‑49 (MQ=21)
                                                                                                                 agcagGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCg                                   >  1:270128/1‑150 (MQ=255)
                                                                                                                 agcagGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCg                                   <  2:355983/150‑1 (MQ=255)
                                                                                                                   cagGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCa                                 <  1:239740/150‑1 (MQ=255)
                                                                                                                    agGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAg                                >  2:960689/1‑150 (MQ=255)
                                                                                                                      gACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGaa                              <  2:930974/150‑1 (MQ=255)
                                                                                                                         aaaTGATGTCAGGTCCGTCGTACTGGACCTGTCCAGt                                                                                                                                            <  2:179094/37‑1 (MQ=255)
                                                                                                                         aaaTGATGTCAGGTCCGTCGTACTGGACCTGTCCAGt                                                                                                                                            >  1:179094/1‑37 (MQ=255)
                                                                                                                         aaaTGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGcc                                                                                      <  2:1224507/91‑1 (MQ=255)
                                                                                                                         aaaTGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGcc                                                                                      >  1:1224507/1‑91 (MQ=255)
                                                                                                                         aaaTGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAgat                           >  1:880746/1‑148 (MQ=255)
                                                                                                                         aaaTGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGtt                           <  1:255959/150‑1 (MQ=255)
                                                                                                                         aaaTGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGtt                           >  2:255959/1‑150 (MQ=255)
                                                                                                                          aaTGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGt                                          <  2:2019223/134‑1 (MQ=255)
                                                                                                                          aaTGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGt                                          >  1:2019223/1‑134 (MQ=255)
                                                                                                                           atgatgTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGa                         <  1:1502633/150‑1 (MQ=255)
                                                                                                                           atgatgTCAGGTACGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTACACACAGCTCATGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGa                         <  2:2142656/150‑1 (MQ=255)
                                                                                                                            tgatgTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAg                        >  2:1446146/1‑150 (MQ=255)
                                                                                                                             gatgTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAgg                       <  1:960689/150‑1 (MQ=255)
                                                                                                                             gatgTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAgg                       <  2:270128/150‑1 (MQ=255)
                                                                                                                                 tCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGcgcc                   >  2:2274300/1‑150 (MQ=255)
                                                                                                                                 tCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGAAGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATACACGGCGACTCAGCACGAGTCCATTCGCAGAAGTTGAGGggcc                   >  1:1154086/1‑150 (MQ=255)
                                                                                                                                     gTCCGTCGTACTGGACAGGTCCCGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCa                                                                                         <  1:61543/76‑1 (MQ=255)
                                                                                                                                     gTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCa                                                                                         >  2:61543/1‑76 (MQ=255)
                                                                                                                                     gTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATc                                                              <  1:2279771/103‑1 (MQ=255)
                                                                                                                                     gTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATc                                                              >  2:2279771/1‑103 (MQ=255)
                                                                                                                                     gTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTc                                    <  1:2065558/129‑1 (MQ=255)
                                                                                                                                     gTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTc                                    >  2:2065558/1‑129 (MQ=255)
                                                                                                                                     gTCCGTCGTACTGGACAGGTCCAGTCCACCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGa                                                      >  1:1340561/1‑111 (MQ=255)
                                                                                                                                     gTCCGTCGTACTGGACAGGTCCAGTCCACCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGa                                                      <  2:1340561/111‑1 (MQ=255)
                                                                                                                                         gtcgtACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATg           >  1:984838/1‑150 (MQ=255)
                                                                                                                                          tcgtACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATgg          >  1:881607/1‑150 (MQ=255)
                                                                                                                                           cggaCTGGACCTGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTc                                                                                <  2:439022/76‑1 (MQ=255)
                                                                                                                                           cggaCTGGACCTGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTc                                                                                >  1:439022/4‑79 (MQ=255)
                                                                                                                                             tACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCtt       >  2:795131/1‑150 (MQ=255)
                                                                                                                                                 ggACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCTTCCtt   >  1:1491440/1‑150 (MQ=255)
                                                                                                                                                  gACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCTTCCTTg  >  1:338408/1‑150 (MQ=255)
                                                                                                                                                  gACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCTTCCTTg  >  2:1573130/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TCCACGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCTTCCTTG  >  NC_002947/4401925‑4402220

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