Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,402,075 G→T 6.9% L375M (CTG→ATG) ‡ PP_3881 ← terminase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,402,0750GT 6.9% 277.2 / 7.9 87L375M (CTG→ATG) ‡PP_3881terminase large subunit
Reads supporting (aligned to +/- strand):  ref base G (42/39);  new base T (3/3);  total (45/42)
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.06e-01

CCACGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCTTCCTTG  >  NC_002947/4401926‑4402220
                                                                                                                                                     |                                                                                                                                                 
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                                                                                                                                         <  2:1699050/142‑1 (MQ=255)
                  cTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCg                                                                                                                                         >  1:1699050/1‑142 (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                                                                                                                             <  1:472506/118‑1 (MQ=255)
                                                      gggTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGa                                                                                                                             >  2:472506/1‑118 (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:973856/150‑1 (MQ=255)
                                                                                   aGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGgaga                                                                <  2:640256/150‑1 (MQ=255)
                                                                                   aGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGgaga                                                                <  1:391774/150‑1 (MQ=255)
                                                                                    gTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAg                                                                                                                                                   <  2:211182/66‑1 (MQ=255)
                                                                                    gTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAg                                                                                                                                                   >  1:211182/1‑66 (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:1481223/1‑150 (MQ=255)
                                                                                              tCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGAc                                                     >  1:930974/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                                                                                                                                          >  1:503683/1‑49 (MQ=21)
                                                                                                            gAACAGCAGGACAAATGATGTCAGGTCCGTCGGACTGGACCTGTCCAGtac                                                                                                                                          <  2:503683/51‑3 (MQ=21)
                                                                                                                agcagGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCg                                   <  2:355983/150‑1 (MQ=255)
                                                                                                                agcagGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCg                                   >  1:270128/1‑150 (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)
                                                                                                                                                     |                                                                                                                                                 
CCACGCCGCGTAAGGCACCTTGTCGCGGGCTTCGCGCTTGTCGAGCTCGTGGTCGGGTATCCAGAAATACGGCAACAGCCGCCAGTGCGGATCCTCGTATGTCGGGTAGAACAGCAGGACAAATGATGTCAGGTCCGTCGTACTGGACAGGTCCAGTCCGCCAACGCTGGGACGATTGCGCAACAGCGACATCGGTACGCGTTCTTCAGCCTGGCTCCAGACATCCCAGGAGATCCACGGCGACTCAGCACCAGTCCATTCGCAGAAGTTGAGGCGCCGCACCATGGCTTCCTTG  >  NC_002947/4401926‑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: