Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,952,898 G→T 8.6% L1036I (CTT→ATT)  recB ← exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,952,8980GT 8.6% 321.5 / 16.0 93L1036I (CTT→ATT) recBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (40/45);  new base T (5/3);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

TTGGATTTATAGTCGAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTT  >  NC_000913/2952753‑2953044
                                                                                                                                                 |                                                                                                                                                  
ttGGATTTATAGTCGAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTg                                                                                                                                                 >  4:228947/1‑149 (MQ=255)
ttGGATTTATAGTCGAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTg                                                                                                                                                 >  3:584988/1‑149 (MQ=255)
ttGGATTTATAGTCGAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTg                                                                                                                                                 >  4:480153/1‑149 (MQ=255)
        ataGTCGAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGa                                                                                                                                         <  3:228947/149‑1 (MQ=255)
            tCGAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATaa                                                                                                                                      >  2:534994/1‑148 (MQ=255)
            tCGAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAg                                                                                                                                     >  3:940050/1‑149 (MQ=255)
             cGAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGc                                                                                                                                    >  3:16821/1‑149 (MQ=255)
              gAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGc                                                                                                                                    <  1:573582/148‑1 (MQ=255)
                gCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCgg                                                                                                                                  <  2:579605/148‑1 (MQ=255)
                 cAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGtt                                                                                                                                <  3:61277/149‑1 (MQ=255)
                 cAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGt                                                                                                                                 <  1:594296/148‑1 (MQ=255)
                     taataaCGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCg                                                                                                                                          <  2:628567/135‑1 (MQ=255)
                     taataaCGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCg                                                                                                                                          >  1:628567/1‑135 (MQ=255)
                     taataaCGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCAc                                                                                                                             >  2:903179/1‑148 (MQ=255)
                        taaCGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAAt                                                                                                                         <  3:480153/149‑1 (MQ=255)
                         aaCGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAAt                                                                                                                         <  2:211173/148‑1 (MQ=255)
                             cccTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGc                                                                                                                     >  1:381856/1‑148 (MQ=255)
                             cccTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGc                                                                                                                     >  1:234703/1‑148 (MQ=255)
                              ccTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAg                                                                                                                   >  4:906260/1‑149 (MQ=255)
                                    tGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATaaaa                                                                                                             <  3:397075/149‑1 (MQ=255)
                                    tGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATaaaa                                                                                                             <  3:165220/149‑1 (MQ=255)
                                        ggtacacCCGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTGTCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTcc                                                                                                         >  4:848493/4‑149 (MQ=255)
                                        ggAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTc                                                                                                          >  1:68926/1‑148 (MQ=255)
                                         gAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTcc                                                                                                         <  2:89988/148‑1 (MQ=255)
                                         gAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCa                                                                                                        <  3:635412/149‑1 (MQ=255)
                                          aaCACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAATCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCAt                                                                                                       >  4:380078/1‑149 (MQ=255)
                                          aaCACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCa                                                                                                        >  2:59201/1‑148 (MQ=255)
                                           acacCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCAtc                                                                                                      >  3:388694/1‑149 (MQ=255)
                                               cAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAACCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAATCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTcc                                                                                                   >  1:829276/1‑148 (MQ=255)
                                                aGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCAc                                                                                                 >  4:195843/1‑149 (MQ=255)
                                                aGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCAc                                                                                                 <  3:195843/149‑1 (MQ=255)
                                                   tCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCt                                                                                               >  2:290495/1‑148 (MQ=255)
                                                     gATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGtt                                                                                            <  3:125086/149‑1 (MQ=255)
                                                      aTAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGTGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGttt                                                                                           >  3:478554/1‑149 (MQ=255)
                                                       tAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCAc                                                                                                                             >  4:332401/1‑114 (MQ=255)
                                                       tAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCAc                                                                                                                             <  3:332401/114‑1 (MQ=255)
                                                       tAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATaaaa                                                                                                             >  4:753895/1‑130 (MQ=255)
                                                       tAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATaaaa                                                                                                             <  3:753895/130‑1 (MQ=255)
                                                          aGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTg                                                                                       >  3:185442/1‑149 (MQ=255)
                                                           gCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTg                                                                                       <  2:358928/148‑1 (MQ=255)
                                                             cTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTgcg                                                                                     >  1:342203/1‑148 (MQ=255)
                                                             cTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTgcg                                                                                     >  1:778861/1‑148 (MQ=255)
                                                              ttttAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAATCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTgcgg                                                                                    <  1:512750/148‑1 (MQ=255)
                                                              ttttAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTgcggg                                                                                   >  4:187752/1‑149 (MQ=255)
                                                               tttAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTgcggg                                                                                   <  2:342203/148‑1 (MQ=255)
                                                                    cATGCCGCGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACCGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGaa                                                                              <  1:664622/148‑1 (MQ=255)
                                                                    cATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGaaa                                                                             <  4:194309/149‑1 (MQ=255)
                                                                    cATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGaa                                                                              <  1:233106/148‑1 (MQ=255)
                                                                      tGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAg                                                                            <  2:820241/148‑1 (MQ=255)
                                                                      tGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGt                                                                           >  4:158700/1‑149 (MQ=255)
                                                                         cACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTg                                                                         <  1:731294/148‑1 (MQ=255)
                                                                              aCCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAATCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCa                                                                    <  2:712474/148‑1 (MQ=255)
                                                                               ccTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAgg                                                                  >  3:440159/1‑149 (MQ=255)
                                                                                cTGCATGAACTCCAGCGGCGGGTAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGc                                                                 >  3:942449/1‑149 (MQ=255)
                                                                                  gCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCt                                                                <  2:234703/148‑1 (MQ=255)
                                                                                      gAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGc                                                           <  4:287696/149‑1 (MQ=255)
                                                                                      gAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGc                                                           >  3:226808/1‑149 (MQ=255)
                                                                                        aCTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCg                                                         <  4:712816/149‑1 (MQ=255)
                                                                                          tCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGt                                                       >  3:470109/1‑149 (MQ=255)
                                                                                            cAGCGGCGGTCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGGTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGttt                                                     >  4:141902/1‑149 (MQ=255)
                                                                                             aGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGttt                                                     >  1:603959/1‑148 (MQ=255)
                                                                                              gcggcgGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCa                                                   <  4:470109/149‑1 (MQ=255)
                                                                                               cggcggGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCAt                                                  <  4:269176/149‑1 (MQ=255)
                                                                                                ggcggGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGc                                                                 >  4:124287/1‑133 (MQ=255)
                                                                                                ggcggGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGc                                                                 <  3:124287/133‑1 (MQ=255)
                                                                                                ggcggGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCAtt                                                 <  4:185442/149‑1 (MQ=255)
                                                                                                   gggCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTgag                                              >  3:939761/1‑149 (MQ=255)
                                                                                                    ggCAGCCTGTGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTgag                                              <  2:98476/148‑1 (MQ=255)
                                                                                                     gCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTgaga                                             <  1:290495/148‑1 (MQ=255)
                                                                                                       aGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGt                                          <  3:187752/149‑1 (MQ=255)
                                                                                                       aGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGt                                          >  4:715474/1‑149 (MQ=255)
                                                                                                         ccTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGc                                        <  3:463595/149‑1 (MQ=255)
                                                                                                            gCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCt                                      >  2:255694/1‑148 (MQ=255)
                                                                                                            gCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCt                                      <  2:174481/148‑1 (MQ=255)
                                                                                                               gATAGTGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTggag                                  <  4:478554/149‑1 (MQ=255)
                                                                                                                        tCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCa                                                                    <  4:800620/106‑1 (MQ=255)
                                                                                                                        tCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCa                                                                    >  3:800620/1‑106 (MQ=255)
                                                                                                                        tCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGt                         <  4:16821/149‑1 (MQ=255)
                                                                                                                          aaaCTGGCGGATTAACGTATCAATCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGa                       >  3:420149/1‑149 (MQ=255)
                                                                                                                          aaaCTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTg                        <  1:583413/148‑1 (MQ=255)
                                                                                                                           aaCTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGa                       <  1:534994/148‑1 (MQ=255)
                                                                                                                           aaCTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGAt                      <  3:35229/149‑1 (MQ=255)
                                                                                                                             cTGGCGGATTAACGTATCAATCTGACTGGCGATAAGCGGTTGACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATc                     >  1:186029/1‑148 (MQ=255)
                                                                                                                               ggCGGATTAACGTATCAATCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCAc                  <  3:380078/149‑1 (MQ=255)
                                                                                                                               ggCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCa                   <  1:230957/148‑1 (MQ=255)
                                                                                                                               ggCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCAc                  <  4:226808/149‑1 (MQ=255)
                                                                                                                                       aaCGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCaa          >  4:230694/1‑149 (MQ=255)
                                                                                                                                            aTCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGt                                                                                             >  1:441900/1‑61 (MQ=255)
                                                                                                                                            aTCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGt                                                                                             <  2:441900/61‑1 (MQ=255)
                                                                                                                                             tCAATCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCg     >  1:913962/1‑148 (MQ=255)
                                                                                                                                               aaGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGtt  >  4:575900/1‑149 (MQ=255)
                                                                                                                                               aaGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGtt  <  3:575900/149‑1 (MQ=255)
                                                                                                                                                aGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGtt  >  2:371766/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TTGGATTTATAGTCGAGCAGGTAATAACGCCCTTCGTGGCGGAACACCAGGTCGATAAAGCCTTTTAACATGCCACGTACCTGCATGAACTCCAGCGGCGGGCAGCCTGCGGATAGCGGGTCAAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTT  >  NC_000913/2952753‑2953044

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