Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,445,549 A→G T480A (ACC→GCC)  tamB → translocation and assembly module subunit TamB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,445,5490AG97.1% 238.2 / ‑2.3 70T480A (ACC→GCC) tamBtranslocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base G (34/33);  minor base T (1/1);  total (35/35)
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.01e-01

GTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTCGGGCCAA  >  NC_000913/4445413‑4445684
                                                                                                                                        |                                                                                                                                       
gTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCt                                                                                                                                      >  1:402584/1‑140 (MQ=255)
   gCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGc                                                                                                                                   >  1:530773/1‑140 (MQ=255)
          tGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGa                                                                                                                            <  2:404710/140‑1 (MQ=255)
          tGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGa                                                                                                                            <  2:829945/140‑1 (MQ=255)
            gCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCTCCTGGCAGATGGACg                                                                                                                          >  1:851285/1‑140 (MQ=255)
            gCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGcc                                                                                                                                       >  1:166994/1‑127 (MQ=255)
                cAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGcc                                                                                                                      <  2:530773/140‑1 (MQ=255)
                   gCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTTCGGCGGCGCCTGGCAGATgg                                                                                                                             >  2:668398/1‑130 (MQ=255)
                   gCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATgg                                                                                                                             <  1:699585/130‑1 (MQ=255)
                   gCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATgg                                                                                                                             <  1:668398/130‑1 (MQ=255)
                   gCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATgg                                                                                                                             >  2:699585/1‑130 (MQ=255)
                   gCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGcc                                                                                                                      >  1:879038/1‑137 (MQ=255)
                   gCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGcc                                                                                                                      <  2:879038/137‑1 (MQ=255)
                   gCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGa                                                                                                                   >  1:327600/1‑140 (MQ=255)
                   gCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGa                                                                                                                   >  1:34399/1‑140 (MQ=255)
                      aTTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGtt                                                                                                                >  2:667310/1‑140 (MQ=255)
                             ggCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTg                                                                                                         >  1:343600/1‑140 (MQ=255)
                                   gTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGt                                                                                                   <  2:343600/140‑1 (MQ=255)
                                    tGAGCTAACGCTTAACGGCATTAACACCGCCAAATAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTa                                                                                                  >  1:1008735/1‑140 (MQ=255)
                                       gCTAACGCTTAACGGCATTACCCCCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACg                                                                                               <  1:205400/140‑1 (MQ=255)
                                       gCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACg                                                                                               <  1:519577/140‑1 (MQ=255)
                                            cGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCgg                                                                                                    >  2:824572/1‑130 (MQ=255)
                                            cGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCgg                                                                                                    <  1:824572/130‑1 (MQ=255)
                                            cGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTaaa                                                                                          <  1:59583/140‑1 (MQ=255)
                                             gCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGACGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAAc                                                                                         <  2:327600/140‑1 (MQ=255)
                                                     gCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACaaa                                                                                 >  1:605317/1‑140 (MQ=255)
                                                       aTTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGt                                                                               >  1:786998/1‑140 (MQ=255)
                                                       aTTAACACCGCCAAAGAGTTCCAGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGt                                                                               >  2:856332/1‑140 (MQ=255)
                                                         tAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGa                                                                             >  1:798793/1‑140 (MQ=255)
                                                                gCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACATTGa                                                                      >  2:998503/1‑140 (MQ=255)
                                                                gCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGa                                                                      <  1:660041/140‑1 (MQ=255)
                                                                 ccAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGAc                                                                     >  2:932949/1‑140 (MQ=255)
                                                                 ccAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGAc                                                                     >  1:534187/1‑140 (MQ=255)
                                                                     agagTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGa                                                                      >  1:818500/1‑135 (MQ=255)
                                                                     agagTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGa                                                                      <  2:818500/135‑1 (MQ=255)
                                                                     agagTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCa                                                                 >  2:442200/1‑140 (MQ=255)
                                                                     agagTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCa                                                                 >  2:906372/1‑140 (MQ=255)
                                                                        gTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATTGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGc                                                              <  1:898957/140‑1 (MQ=255)
                                                                           cccGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGa                                                           >  1:63018/1‑140 (MQ=255)
                                                                                  tGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCTCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCaa                                                    <  2:851285/140‑1 (MQ=255)
                                                                                  tGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCaa                                                    <  2:828676/140‑1 (MQ=255)
                                                                                  tGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCaa                                                    <  2:534187/140‑1 (MQ=255)
                                                                                  tGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCaa                                                    <  1:932949/140‑1 (MQ=255)
                                                                                  tGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCaa                                                    <  2:63018/140‑1 (MQ=255)
                                                                                         cGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTtat                                             >  2:684073/1‑140 (MQ=255)
                                                                                          gAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTtata                                            <  1:684073/140‑1 (MQ=255)
                                                                                            aaCTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATg                                          >  2:818642/1‑140 (MQ=255)
                                                                                              cTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGTGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCa                                        >  2:925065/1‑140 (MQ=255)
                                                                                              cTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCa                                        <  2:798793/140‑1 (MQ=255)
                                                                                              cTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCa                                        <  1:1086098/140‑1 (MQ=255)
                                                                                                   tGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGga                                   >  1:252667/1‑140 (MQ=255)
                                                                                                   tGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGga                                   <  1:667310/140‑1 (MQ=255)
                                                                                                    ggCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGgat                                  <  1:998503/140‑1 (MQ=255)
                                                                                                    ggCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGgat                                  <  2:252667/140‑1 (MQ=255)
                                                                                                    ggCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGgat                                  <  2:450340/140‑1 (MQ=255)
                                                                                                      cTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACa                                                  >  2:371522/1‑122 (MQ=255)
                                                                                                      cTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACa                                                  <  1:371522/122‑1 (MQ=255)
                                                                                                      cTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGgatga                                >  2:385133/1‑140 (MQ=255)
                                                                                                          aTTAAAACCCGCGGTAGCCTGTACGGCGGTGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATccc                            <  1:925065/140‑1 (MQ=255)
                                                                                                           ttAAAACCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCa                           <  1:856332/140‑1 (MQ=255)
                                                                                                               aaCCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGc                       >  2:566539/1‑140 (MQ=255)
                                                                                                                aCCCGCGGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCt                      <  1:566539/140‑1 (MQ=255)
                                                                                                                     cgGTAGCCTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATc                 >  1:607100/1‑140 (MQ=255)
                                                                                                                           ccTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAAc           <  2:150439/140‑1 (MQ=255)
                                                                                                                           ccTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAAc           <  2:607100/140‑1 (MQ=255)
                                                                                                                            cTGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACt          <  2:605317/140‑1 (MQ=255)
                                                                                                                             tGTACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTc         >  2:1052079/1‑140 (MQ=255)
                                                                                                                               tACGGCGGCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTCgg       >  2:743996/1‑140 (MQ=255)
                                                                                                                                 cggcggCGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTCGGGc     <  1:791662/140‑1 (MQ=255)
                                                                                                                                    cggcGCCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTCGGGCCaa  >  2:828685/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATGGACGTGCCGGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTCGGGCCAA  >  NC_000913/4445413‑4445684

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