Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,519,557 A→G 5.3% W383R (TGG→CGG)  hofQ ← DNA catabolic putative fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,519,5570AG 5.3% 236.2 / 2.8 76W383R (TGG→CGG) hofQDNA catabolic putative fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base A (38/33);  new base G (2/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

CGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGC  >  NC_000913/3519409‑3519704
                                                                                                                                                    |                                                                                                                                                   
cgctTCAAATAGATTAAACACGTTCAGCGTAAAAACATTACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCa                                                                                                                                                     <  8:173022/149‑1 (MQ=255)
   gTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAggg                                                                                                                                                  <  4:122932/149‑1 (MQ=255)
   gTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAggg                                                                                                                                                  <  8:303923/149‑1 (MQ=255)
    tCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCTTTTACGGCTGTTTACTCACTGGATACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGACCCAGGGa                                                                                                                                                 >  3:186312/1‑149 (MQ=255)
             tGAAACACGTTCAGCGTAAAAACATTACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGAATGTCGCCaa                                                                                                                                       <  7:130071/149‑1 (MQ=255)
              gAAACACGTTCAGCGTAAAAACATTACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACAACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCACGAACCAGGGAATGTCGcc                                                                                                                                         >  3:221643/1‑147 (MQ=255)
                   acGTTCAGCGTAAAAACATTACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAAc                                                                                                                                  <  1:173616/149‑1 (MQ=255)
                   acGTTCAGCGTAAAAACATTACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAAc                                                                                                                                  <  1:168311/149‑1 (MQ=255)
                   acGTTCAGCGTAAAAACATTACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAAc                                                                                                                                  <  7:347383/149‑1 (MQ=255)
                      ttCAGCGTAAAAACATTACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACgg                                                                                                                                >  4:296789/1‑148 (MQ=255)
                             tAAAAACATTACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTAt                                                                                                                        <  5:66560/149‑1 (MQ=255)
                               aaaaCATTACCGCGTTTACGGCTGATTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATc                                                                                                                       <  4:186312/148‑1 (MQ=255)
                                    aTTACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGAcc                                                                                                                 >  6:115811/1‑149 (MQ=255)
                                     ttACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGAccc                                                                                                                >  7:26984/1‑149 (MQ=255)
                                          gCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGattt                                                                                                           <  3:296789/149‑1 (MQ=255)
                                                 cGGCTGTTTACTCACTGGAAACCAGTCGTTGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCTGGGAATGTCGCCAAGCAACGGTACGCTATCCTGTCCCGATTTAttttt                                                                                                     <  8:26984/148‑1 (MQ=255)
                                                    cTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACgg                                                                                                 <  5:115811/149‑1 (MQ=255)
                                                              aCTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATgc                                                                                       >  8:262629/1‑149 (MQ=255)
                                                              aCTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATgc                                                                                       >  6:208652/1‑149 (MQ=255)
                                                                  gAACCCAGTCGTGGCGTGATAACCACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCGGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATgccgcc                                                                                   >  3:8409/1‑149 (MQ=255)
                                                                  gAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCGGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATgccgcc                                                                                   >  3:8413/1‑149 (MQ=255)
                                                                  gAAACCAGTCGTGGCGTGATAAACACCACTAACTAGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGCACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATgccgcc                                                                                   >  5:299557/1‑149 (MQ=255)
                                                                       cAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAgcg                                                                              >  8:30760/1‑149 (MQ=255)
                                                                        aGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTc                                                                                                                                            <  6:225449/86‑1 (MQ=255)
                                                                        aGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTc                                                                                                                                            >  5:225449/1‑86 (MQ=255)
                                                                                 gTGATAAACACCACTAACTCGCGTCGTTAATCTTCTTTTACGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGACGCCCAGCGCCAACGttt                                                                     >  5:47205/1‑148 (MQ=255)
                                                                                   gATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTtctc                                                                  >  3:88083/1‑149 (MQ=255)
                                                                                    aTAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACg                                                                                                  <  8:21706/116‑1 (MQ=255)
                                                                                    aTAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACg                                                                                                  >  7:21706/1‑116 (MQ=255)
                                                                                    aTAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTcc                                                                 <  2:95179/149‑1 (MQ=255)
                                                                                        acacCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCtt                                                             >  7:44497/1‑149 (MQ=255)
                                                                                         caccacTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTTCGGGTAAAAATgcc                                                                                      <  8:149739/123‑1 (MQ=255)
                                                                                         caccacTAACTCGCGTCGTTCATCTTCTTTTACGTCATCACGAAATAATAGCCCGCACCAGGGAATGTCGACAAGCAACGGTACGCAATCCTGAACCGATATAAATTTACGCGTAAAAATgcc                                                                                      >  7:149739/1‑123 (MQ=255)
                                                                                          accacTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCtttt                                                           <  3:39497/149‑1 (MQ=255)
                                                                                          accacTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCtttt                                                           <  3:353796/149‑1 (MQ=255)
                                                                                          accacTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCtttt                                                           >  4:353796/1‑149 (MQ=255)
                                                                                            cacTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGa                                                         >  8:29599/1‑149 (MQ=255)
                                                                                                 aCTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTc                                                     >  7:367674/1‑148 (MQ=255)
                                                                                                           ttCATCTTCTTTTCCGTCATGACGAATTTATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTTTCCTGACCCTATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGtt                                           <  6:47205/148‑1 (MQ=255)
                                                                                                           ttCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGttt                                          <  5:208652/149‑1 (MQ=255)
                                                                                                             cATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGTGTTTc                                         <  7:30760/148‑1 (MQ=255)
                                                                                                              aTCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCtttt                                                           <  4:112799/129‑1 (MQ=255)
                                                                                                              aTCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCtttt                                                           >  3:112799/1‑129 (MQ=255)
                                                                                                              aTCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCg                                        >  3:264737/1‑148 (MQ=255)
                                                                                                              aTCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGa                                       >  1:185537/1‑149 (MQ=255)
                                                                                                               tcttctTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCt                                                <  1:276866/139‑1 (MQ=255)
                                                                                                               tcttctTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGACGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCt                                                >  2:276866/1‑139 (MQ=255)
                                                                                                                cttcttTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGGCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGAtc                                     >  2:273630/1‑149 (MQ=255)
                                                                                                                cttcttTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGAtc                                     >  6:59520/1‑149 (MQ=255)
                                                                                                                 ttcttTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACATGCGTTTCGAtct                                    >  3:330836/1‑149 (MQ=255)
                                                                                                                   cttTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGc                                              <  7:89136/137‑1 (MQ=255)
                                                                                                                   cttTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGc                                              >  8:89136/1‑137 (MQ=255)
                                                                                                                    ttttCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTATAAATGCCGCCCAGCGCAAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCTCTCTCCt                                 >  2:282626/1‑149 (MQ=255)
                                                                                                                        ccGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCt                              <  7:262629/148‑1 (MQ=255)
                                                                                                                         cGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTa                            <  8:367674/149‑1 (MQ=255)
                                                                                                                          gTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCTTTTCGATCTCCTGCTTAt                           <  8:44497/149‑1 (MQ=255)
                                                                                                                               gACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCaa                                                                          >  1:316427/1‑97 (MQ=255)
                                                                                                                               gACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCaa                                                                          <  2:316427/97‑1 (MQ=255)
                                                                                                                               gACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATc                      >  7:185063/1‑149 (MQ=255)
                                                                                                                               gACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATc                      >  8:371865/1‑149 (MQ=255)
                                                                                                                               gACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTACGATCTCCTGCTTATCAATc                      >  7:323142/1‑149 (MQ=255)
                                                                                                                                aCGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCg                     >  2:124331/1‑149 (MQ=255)
                                                                                                                                 cGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGc                    <  4:330836/149‑1 (MQ=255)
                                                                                                                                  gAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGc                    >  4:271499/1‑148 (MQ=255)
                                                                                                                                  gAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGc                    >  5:230063/1‑148 (MQ=255)
                                                                                                                                  gAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGc                    >  6:12013/1‑148 (MQ=255)
                                                                                                                                   aaaTAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCa                  <  4:88083/149‑1 (MQ=255)
                                                                                                                                   aaaTAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCa                  <  5:59520/149‑1 (MQ=255)
                                                                                                                                    aataatTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAg                 <  4:264737/149‑1 (MQ=255)
                                                                                                                                       aatTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCAc              >  7:101804/1‑149 (MQ=255)
                                                                                                                                           gCCCGAACCGGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTc           <  4:8413/148‑1 (MQ=255)
                                                                                                                                           gCCCGAACCGGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGGTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTc           <  4:8409/148‑1 (MQ=255)
                                                                                                                                           gCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTc           <  6:299557/148‑1 (MQ=255)
                                                                                                                                              cGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCa       >  1:112097/1‑149 (MQ=255)
                                                                                                                                               gAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCAt      >  8:363241/1‑149 (MQ=255)
                                                                                                                                                 aCCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCg    <  1:273630/149‑1 (MQ=255)
                                                                                                                                                   cAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGc  <  2:185537/149‑1 (MQ=255)
                                                                                                                                                    aGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGAtctc                                   <  2:137414/115‑1 (MQ=255)
                                                                                                                                                    aGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGAtctc                                   >  1:137414/1‑115 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCTTTTACGGCTGTTTACTCACTGGAAACCAGTCGTGGCGTGATAAACACCACTAACTCGCGTCGTTCATCTTCTTTTCCGTCATGACGAAATAATTGCCCGAACCAGGGAATGTCGCCAAGCAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGC  >  NC_000913/3519409‑3519704

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