Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,623,709 T→G I104S (ATT→AGT)  yhhH → putative NTF2 fold immunity protein for polymorphic toxin RhsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,623,7090TG98.7% 502.9 / ‑1.5 148I104S (ATT→AGT) yhhHputative NTF2 fold immunity protein for polymorphic toxin RhsB
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (72/74);  total (74/74)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

Only 100 of 148 total aligned reads displayed.

GTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATC  >  NC_000913/3623613‑3623805
                                                                                                |                                                                                                
gTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGc                                                                                               >  1:593098/1‑100 (MQ=255)
gTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACAGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGc                                                                                               >  2:2908388/1‑100 (MQ=255)
 ttATAGACAGGATGCCGCAGCAGAAGAAAAACCATACGAAATCACTAAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCt                                                                                              <  1:1018287/100‑1 (MQ=255)
   aTGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCtgg                                                                                            >  1:3114154/1‑100 (MQ=255)
    tGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCtggt                                                                                           >  2:29899/1‑100 (MQ=255)
     ggACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCtggtg                                                                                          >  2:2501039/1‑100 (MQ=255)
     ggACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCtggtg                                                                                          >  1:641789/1‑100 (MQ=255)
     ggACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCtggtg                                                                                          >  2:2573696/1‑100 (MQ=255)
       aCAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTggg                                                                                        <  2:3089379/100‑1 (MQ=255)
       aCAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTggg                                                                                        <  1:1377472/100‑1 (MQ=255)
        cAGGATGCCGCAGAAGACGAAAAACCATACGAAATCACTGAGTTTACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTgggg                                                                                       <  1:3749938/100‑1 (MQ=255)
        cAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTgggg                                                                                       <  1:1674693/100‑1 (MQ=255)
          ggATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAGGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGtt                                                                                     >  1:1224627/1‑100 (MQ=255)
          ggATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGtt                                                                                     >  1:1618385/1‑100 (MQ=255)
          ggATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGtt                                                                                     <  1:891588/100‑1 (MQ=255)
          ggATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGtt                                                                                     <  1:891275/100‑1 (MQ=255)
          ggATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGtt                                                                                     >  1:1555643/1‑100 (MQ=255)
          ggATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGtt                                                                                     >  1:764177/1‑100 (MQ=255)
          ggATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGtt                                                                                     >  1:1555946/1‑100 (MQ=255)
            aTGCCGCAGGAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGtttt                                                                                   <  2:2142441/100‑1 (MQ=255)
            aTGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGtttt                                                                                   <  2:2280977/100‑1 (MQ=255)
             ngCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGttttt                                                                                  >  1:2422806/2‑100 (MQ=255)
              gCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTtta                                                                                 <  1:1812823/100‑1 (MQ=255)
               ccGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTttat                                                                                <  1:3325569/100‑1 (MQ=255)
                cGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTttatt                                                                               >  2:2382354/1‑100 (MQ=255)
                 gCAGAAGAAGAAAGACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTttatta                                                                              <  2:169187/100‑1 (MQ=255)
                 gCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTttatta                                                                              <  2:575673/100‑1 (MQ=255)
                 gCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTttatta                                                                              >  2:1934126/1‑100 (MQ=255)
                  cAGAAGGAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGTTTGTTGAAGGTACCATGCACTCAGTTCAAAGTGCTGGTGGTGTTTTtattat                                                                             >  2:2404783/1‑100 (MQ=255)
                  cAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTtattat                                                                             <  1:2369716/100‑1 (MQ=255)
                  cAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTtattat                                                                             <  1:2316743/100‑1 (MQ=255)
                  cAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTtattat                                                                             <  2:2654414/100‑1 (MQ=255)
                   agaagaagaAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATtata                                                                            <  1:970244/100‑1 (MQ=255)
                   agaagaagaAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATtata                                                                            <  1:1711300/100‑1 (MQ=255)
                     aagaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGa                                                                          >  1:3046336/1‑100 (MQ=255)
                      agaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGaa                                                                         <  2:1240120/100‑1 (MQ=255)
                       gaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGaaa                                                                        >  1:2123558/1‑100 (MQ=255)
                       gaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGaaa                                                                        >  2:3082231/1‑100 (MQ=255)
                       gaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGaaa                                                                        >  1:1068251/1‑100 (MQ=255)
                       gaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGATGTTGAAGGTACCATTCACGCAGACCAAAGGGCTCGTGGGGTTTTTATTATAGaaa                                                                        >  2:3785183/1‑100 (MQ=255)
                       gaagaaAAACCATACGAAATCACTGAGTTAACAACAAGATGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGaaa                                                                        >  1:1148045/1‑100 (MQ=255)
                         agaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATa                                                                      <  1:1727009/100‑1 (MQ=255)
                         agaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATa                                                                      >  1:943952/1‑100 (MQ=255)
                         agaaAAACCATACGAAATCACGGAGTTAACAACAAGTTGGGTCGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATa                                                                      >  1:3479468/1‑100 (MQ=255)
                         agaaAAAACATACGAAATCACTGAGTTAACAACAAGTTAGGTTGTTGAAGGTACCATTCACTCAGAACAAAGTGCTGGTGGGGTTTTTATTATAGAAATa                                                                      <  2:1388300/100‑1 (MQ=255)
                          gaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAg                                                                     >  2:669748/1‑100 (MQ=255)
                          gaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAg                                                                     >  1:3000628/1‑100 (MQ=255)
                               aCCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCaaa                                                                <  2:3681035/100‑1 (MQ=255)
                                   tACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAatg                                                            <  1:451237/100‑1 (MQ=255)
                                     cGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAatgat                                                          <  2:2955929/100‑1 (MQ=255)
                                     cGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAatgat                                                          >  1:267166/1‑100 (MQ=255)
                                     cGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATCACAGAAATAGGCAAAAatgat                                                          <  1:3483956/100‑1 (MQ=255)
                                        aaTCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATggg                                                       >  1:27081/1‑100 (MQ=255)
                                        aaTCACTGAGTTAACAACAAGATGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTAAAAAAAAAAAAAAAAATGATggg                                                       <  2:1089558/100‑1 (MQ=255)
                                              tGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAAtt                                                 >  1:1618483/1‑100 (MQ=255)
                                               gAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTc                                                >  1:3882692/1‑100 (MQ=255)
                                                aGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCt                                               <  1:2398177/100‑1 (MQ=255)
                                                  ttAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGa                                             >  2:1050071/1‑100 (MQ=255)
                                                   tAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGCATTCTGaa                                            >  1:2902348/1‑100 (MQ=255)
                                                    aacaacAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAt                                           <  1:4059690/100‑1 (MQ=255)
                                                    aacaacAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAt                                           <  2:1526509/100‑1 (MQ=255)
                                                     acaacaAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATGATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAtt                                          <  2:3905464/100‑1 (MQ=255)
                                                        acaaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTg                                       <  2:2175674/100‑1 (MQ=255)
                                                         caaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTgg                                      <  2:2425659/100‑1 (MQ=255)
                                                         caaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTgg                                      <  2:1976234/100‑1 (MQ=255)
                                                         caaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTgg                                      <  1:3910023/100‑1 (MQ=255)
                                                         caaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTAAAAAAAAAAAAAAAAATGATGGGAGAATTCTGAATTTTgg                                      <  2:416546/100‑1 (MQ=255)
                                                          aaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGAAAAAATGATGGGAGAATTCTGAATTTTGGc                                     <  2:1692373/100‑1 (MQ=255)
                                                           aGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGcc                                    >  1:1462885/1‑100 (MQ=255)
                                                              tGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACg                                 <  2:3067521/100‑1 (MQ=255)
                                                                  ttgttgAAGGTACCATTCACTCAGACCAAAGTGGTGGTTGGGGTTTGAAAATACAAATAGGCAAAAATGATGGGAGAATTCGGAATTTTGGCCACGGaaa                             <  2:1936847/100‑1 (MQ=255)
                                                                   tgttgCAGGTACCATTCACTCAGACCAAAGTGCTGGCGGGGTTTTTATTATATAAATAGGCAACAATGATGGGAGAATTCTGAATTTTGGCCACGGaaaa                            >  2:3404180/1‑100 (MQ=255)
                                                                   tgttgAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGaaaa                            >  2:3671453/1‑100 (MQ=255)
                                                                    gttgAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAAt                           >  1:1195063/1‑100 (MQ=255)
                                                                     ttgAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATa                          <  2:1050524/100‑1 (MQ=255)
                                                                     ttgAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATa                          <  1:1846562/100‑1 (MQ=255)
                                                                       gAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAg                        >  1:2028115/1‑100 (MQ=255)
                                                                       gAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAg                        >  1:1370351/1‑100 (MQ=255)
                                                                       gAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAg                        >  1:1312536/1‑100 (MQ=255)
                                                                       gAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAAAAAAAAAAAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAg                        <  2:550648/100‑1 (MQ=255)
                                                                          ggTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTAt                     >  1:1059812/1‑100 (MQ=255)
                                                                          ggTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTAt                     >  2:2588917/1‑100 (MQ=255)
                                                                           gTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATg                    <  1:1482575/100‑1 (MQ=255)
                                                                            tACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGt                   >  1:1983709/1‑100 (MQ=255)
                                                                            tACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGt                   <  1:2603956/100‑1 (MQ=255)
                                                                              ccATTCACTCCGTCCAAAGTGCTGGTGGGGCTTTTATTCTAGAAATACGCAAAAATGATGGGAGTATTATTACTTTTTTGCACTTAAATTATGTATGTaa                 >  2:3564982/1‑100 (MQ=255)
                                                                              ccATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTaa                 <  2:3065971/100‑1 (MQ=255)
                                                                               cATTCACTCAGACCAAAGTGCTGGTGGTGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTATGAATTTTGGCCACGGAAAATAAGTATGTAAg                >  2:1276885/1‑100 (MQ=255)
                                                                               cATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAg                >  1:3077270/1‑100 (MQ=255)
                                                                                aTTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAgg               >  1:2640618/1‑100 (MQ=255)
                                                                                  tctcTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGAt             <  2:1187514/97‑1 (MQ=255)
                                                                                  tCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGAt             <  2:256892/100‑1 (MQ=255)
                                                                                   cACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATa            <  1:1649094/100‑1 (MQ=255)
                                                                                     cTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATaag          >  1:3450986/1‑99 (MQ=255)
                                                                                      taaGACGAAAGTGGTGGTGGGGTTTTTATTATAGAAATAGAAAAAAATGCTGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATc         <  2:3935509/98‑1 (MQ=255)
                                                                                      tCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATc         <  1:2076199/100‑1 (MQ=255)
                                                                                      tCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGAAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATc         <  2:2906218/100‑1 (MQ=255)
                                                                                         gACCAAAGTGCTGGTGGGGTTTTTATTATGGAAATAGGCAAAAATGATGGGAGAATTCTGAAGTTTGGCCACGGAAAAGAAGTATGTAAGGATAATCaaa      >  1:3802232/1‑100 (MQ=255)
                                                                                         gACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCaaa      >  2:2685966/1‑100 (MQ=255)
                                                                                          aCCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAAc     <  2:886298/100‑1 (MQ=255)
                                                                                                |                                                                                                
GTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATC  >  NC_000913/3623613‑3623805

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