Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,374,260 G→A 5.5% A77V (GCA→GTA)  nanR ← sialic acid‑inducible nan operon repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,374,2600GA 5.5% 240.3 / 2.3 73A77V (GCA→GTA) nanRsialic acid‑inducible nan operon repressor
Reads supporting (aligned to +/- strand):  ref base G (37/32);  new base A (2/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCATTTCGGAGAGTTTTTTACGC  >  NC_000913/3374115‑3374406
                                                                                                                                                 |                                                                                                                                                  
ttCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCa                                                                                                                                                 >  4:831892/1‑149 (MQ=255)
  cGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAgc                                                                                                                                               <  1:747125/149‑1 (MQ=255)
       tGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTc                                                                                                                                          <  2:301291/149‑1 (MQ=255)
        gggCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCa                                                                                                                                         >  2:560878/1‑149 (MQ=255)
         ggCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCAc                                                                                                                                        <  2:79045/149‑1 (MQ=255)
              tCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAgcg                                                                                                                                              <  3:804611/138‑1 (MQ=255)
              tCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAgcg                                                                                                                                              >  4:804611/1‑138 (MQ=255)
                 cACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGaa                                                                                                                                >  3:198779/1‑149 (MQ=255)
                       gATGAGAAAGGAAATCTTTCGCCCTGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGAACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGACAGCGCTTCACGCACCGAAGGACGc                                                                                                                          >  2:151834/1‑149 (MQ=255)
                        aTGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGcc                                                                                                                         >  1:71151/1‑149 (MQ=255)
                             aaaGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACg                                                                                                                    >  2:472658/1‑149 (MQ=255)
                              aaGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGt                                                                                                                   >  4:317807/1‑149 (MQ=255)
                               aGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGTGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGtt                                                                                                                  <  2:616053/149‑1 (MQ=255)
                               aGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGtt                                                                                                                  >  2:33136/1‑149 (MQ=255)
                                ggAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTACCAGCGCTTCACGCACCGAAGGACGCCCGACGTTa                                                                                                                 >  4:688633/1‑149 (MQ=255)
                                ggAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTACCAGCGCTTCACGCACCGAAGGACGCCCGACGTTa                                                                                                                 >  2:604149/1‑149 (MQ=255)
                                 gAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTaa                                                                                                                <  1:195661/149‑1 (MQ=255)
                                 gAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTaa                                                                                                                >  1:434167/1‑149 (MQ=255)
                                   aaTCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAg                                                                                                              >  1:282205/1‑149 (MQ=255)
                                     tCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGaa                                                                                                            <  1:472658/149‑1 (MQ=255)
                                      cTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAAc                                                                                                           <  2:71151/149‑1 (MQ=255)
                                       tttCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACg                                                                                                          >  3:633062/1‑149 (MQ=255)
                                        ttCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGc                                                                                                         <  2:771216/149‑1 (MQ=255)
                                         tCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGcc                                                                                                        <  3:440893/149‑1 (MQ=255)
                                                gCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTACCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGtt                                                                                                 <  1:604149/149‑1 (MQ=255)
                                                gCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTACCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGtt                                                                                                 <  3:688633/149‑1 (MQ=255)
                                                     aaaGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGt                                                                                            <  1:846734/149‑1 (MQ=255)
                                                        gCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCa                                                                                         <  3:318168/149‑1 (MQ=255)
                                                          tCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGTACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGa                                                                                       >  1:484902/1‑149 (MQ=255)
                                                          tCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGa                                                                                       >  3:372896/1‑149 (MQ=255)
                                                            aCCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACg                                                                                     >  3:684517/1‑149 (MQ=255)
                                                             ccGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACgg                                                                                    <  3:391745/149‑1 (MQ=255)
                                                              cGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGt                                                                                   >  2:340812/1‑149 (MQ=255)
                                                               gatgatAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTa                                                                                  <  1:560878/149‑1 (MQ=255)
                                                                  gatAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAAtt                                                                               >  2:599428/1‑149 (MQ=255)
                                                                  gatAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAAtt                                                                               <  1:504061/149‑1 (MQ=255)
                                                                  gatAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAAtt                                                                               >  4:506070/1‑149 (MQ=255)
                                                                         tCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCAcc                                                                        <  4:372896/149‑1 (MQ=255)
                                                                                       cgagacgcgagCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCac                                                          <  1:32374/149‑1 (MQ=255)
                                                                                        gagacgcgagCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCacg                                                         >  1:193107/1‑149 (MQ=255)
                                                                                        gagacgcgagCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCacg                                                         >  3:271804/1‑149 (MQ=255)
                                                                                                     gTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGatcatc                                            <  3:831892/149‑1 (MQ=255)
                                                                                                     gTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGatcatc                                            >  1:78900/1‑149 (MQ=255)
                                                                                                      ttCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCt                                           >  1:161243/1‑149 (MQ=255)
                                                                                                        cGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGTGGATCATCTGt                                         <  3:506070/149‑1 (MQ=255)
                                                                                                         gCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGtt                                        >  3:827816/1‑149 (MQ=255)
                                                                                                           cGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGtt                                                                            >  4:770695/1‑111 (MQ=255)
                                                                                                           cGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGtt                                                                            <  3:770695/111‑1 (MQ=255)
                                                                                                            gttgttTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCa                                     >  4:140642/1‑149 (MQ=255)
                                                                                                             ttgttTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCac                                                          <  3:454402/127‑1 (MQ=255)
                                                                                                             ttgttTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCac                                                          >  4:454402/1‑127 (MQ=255)
                                                                                                                tttATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGctc                                 >  3:576749/1‑149 (MQ=255)
                                                                                                                tttATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGctc                                 >  3:696292/1‑149 (MQ=255)
                                                                                                                    tttGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCtcttct                             >  2:849617/1‑149 (MQ=255)
                                                                                                                      tGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCttcttc                           <  1:11011/149‑1 (MQ=255)
                                                                                                                      tGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCttcttc                           <  4:198779/149‑1 (MQ=255)
                                                                                                                        cACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTcca                         <  2:434167/149‑1 (MQ=255)
                                                                                                                        cACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCACGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTcca                         <  3:251422/149‑1 (MQ=255)
                                                                                                                          ccAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTccacc                       >  3:99238/1‑149 (MQ=255)
                                                                                                                            aGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCAt                     >  4:470258/1‑149 (MQ=255)
                                                                                                                               ccTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCATTTc                  <  2:193107/149‑1 (MQ=255)
                                                                                                                                cTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCATTTCg                 <  2:496086/149‑1 (MQ=255)
                                                                                                                                cTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCATTTCg                 >  3:126087/1‑149 (MQ=255)
                                                                                                                                 tttGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCacga                                                        <  3:36949/109‑1 (MQ=255)
                                                                                                                                 tttGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCacga                                                        >  4:36949/1‑109 (MQ=255)
                                                                                                                                 tttGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCATTTCgg                >  3:509787/1‑149 (MQ=255)
                                                                                                                                      gTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCATTTCGGAGAGt           >  2:95203/1‑149 (MQ=255)
                                                                                                                                      gTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCATTTCGGAGAGt           >  3:353108/1‑149 (MQ=255)
                                                                                                                                         ttAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCATTTCGGAGAGtttt        <  4:744667/149‑1 (MQ=255)
                                                                                                                                          tAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCATTTCGGAGAGttttt       <  1:340812/149‑1 (MQ=255)
                                                                                                                                             cGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCATTTCGGAGAGTTTTTTAc    <  2:177548/149‑1 (MQ=255)
                                                                                                                                              gCTGCCATCGCTTCACGCACCGAAGCACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCATTTCGGAGAGTTTTTTAcg   >  2:120794/1‑149 (MQ=255)
                                                                                                                                               cTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCATTTCGGAGAGTTTTTTAcgc  <  4:126087/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TTCGAAATGGGCAATCCCACCGGGATGAGAAAGGAAATCTTTCGCCATGCCGGAAAGCTCACCGATGATAGTGTCCGCAGAAGGACGCGAGACGCGAGCGCGTTCGCCGTTGTTTATTTGCACCAGACCTTTGCGTTTTAACGCTGCCAGCGCTTCACGCACCGAAGGACGCCCGACGTTAAAGAACGCCATCAGTTCGCGTTCAGACGGTAATTGTTCACCTTCGCCAAATTCACGACGGCGGATCATCTGTTCCAGCTCTTCTTCCACCATTTCGGAGAGTTTTTTACGC  >  NC_000913/3374115‑3374406

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: ‑