Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 942,577 T→A 8.0% L540Q (CTG→CAG)  dmsA → dimethyl sulfoxide reductase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913942,5770TA 8.0% 217.9 / 7.7 75L540Q (CTG→CAG) dmsAdimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base A (3/3);  total (37/38)
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.20e-01

TCTGAAATCAACCGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGT  >  NC_000913/942429‑942719
                                                                                                                                                    |                                                                                                                                              
tCTGAAATCAACCGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTg                                                                                                                                               <  1:107553/150‑1 (MQ=255)
 cTGAAATCAACCGTTCCCTTGAAATCCTTCATGATGATAAGATGTCCGAGCTGATTGTGGTTCTCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGTCGTTGCGCTgg                                                                                                                                              <  2:130949/150‑1 (MQ=255)
   gAAATCAACCGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGAt                                                                                                                                            <  1:959923/150‑1 (MQ=255)
    aaaTCAACCGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATg                                                                                                                                           <  1:1517795/150‑1 (MQ=255)
      aTCAACCGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCa                                                                                                                                         <  2:429014/150‑1 (MQ=255)
        cAACCGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGGTGCATc                                                                                                                                       <  1:407911/150‑1 (MQ=255)
         aaCCGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATcc                                                                                                                                      <  2:54527/150‑1 (MQ=255)
          aCCGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCt                                                                                                                                     >  2:865627/1‑150 (MQ=255)
           ccGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTg                                                                                                                                    <  1:401484/150‑1 (MQ=255)
           ccGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCAGGATGCATccag                                                                                                                                    >  2:748412/1‑148 (MQ=255)
              tACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGGTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCAGGATGCATCCAGCgc                                                                                                                                 <  1:1536018/150‑2 (MQ=255)
               aCCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCggg                                                                                                                                <  1:1226195/150‑1 (MQ=255)
                cccATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGa                                                                                                                               <  1:1176731/150‑1 (MQ=255)
                cccATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGa                                                                                                                               >  2:115282/1‑150 (MQ=255)
                  cATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAat                                                                                                                             <  2:841798/150‑1 (MQ=255)
                  cATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCAGGATGCATCCAGCgcaaag                                                                                                                             <  1:748412/150‑6 (MQ=255)
                       aaTCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTc                                                                                                                        <  1:865627/150‑1 (MQ=255)
                         tCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGGTGACATCCTGCTGCCTGACTTGACCCCTTCCGAATAGATGTACTTTGCGCTGGATGCATCCTGCGGGAATATGTCtt                                                                                                                      <  2:134818/150‑1 (MQ=255)
                          ccTTCAGGATTATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTa                                                                                                                     <  1:1477710/150‑1 (MQ=255)
                          ccTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTa                                                                                                                     <  1:548465/150‑1 (MQ=255)
                           cTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGAACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTAc                                                                                                                    <  1:1175119/150‑1 (MQ=255)
                              aaGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTg                                                                                                                 <  2:173277/149‑1 (MQ=255)
                                   tgatAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGAtttt                                                                                                            >  1:110882/1‑150 (MQ=255)
                                      tAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCaa                                                                                                         >  1:791883/1‑150 (MQ=255)
                                      tAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCaa                                                                                                         <  1:1293054/150‑1 (MQ=255)
                                         gaagTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGa                                                                                                      >  2:1176237/1‑150 (MQ=255)
                                              gCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAgg                                                                                                 >  1:859879/1‑150 (MQ=255)
                                                            gTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCgcgc                                                                                   <  2:110882/150‑1 (MQ=255)
                                                             ttATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCt                                                                                  >  2:485515/1‑150 (MQ=255)
                                                                    cTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGGTTAAACCGCGCTTTGAATg                                                                           >  1:729216/1‑150 (MQ=255)
                                                                      gCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTa                                                                         >  2:884213/1‑150 (MQ=255)
                                                                       ccACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTaa                                                                        >  1:1342070/1‑150 (MQ=255)
                                                                        cacaTGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGa                                                                              >  2:1371328/1‑143 (MQ=255)
                                                                        cacaTGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGa                                                                              <  1:1371328/143‑1 (MQ=255)
                                                                        cacaTGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAg                                                                       >  2:1153583/1‑150 (MQ=255)
                                                                            tGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAGCGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCa                                                                   <  2:354562/150‑1 (MQ=255)
                                                                               ccTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGACCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCt                                                                >  1:1237095/1‑150 (MQ=255)
                                                                                  cATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATg                                                             <  2:1098440/150‑1 (MQ=255)
                                                                                  cATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGAGCAGGTGATTAAACCGCGCTTTGAATGTAAGCCCATATATg                                                             >  1:124955/1‑150 (MQ=255)
                                                                                     cGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAAt                                                         >  2:338190/1‑150 (MQ=255)
                                                                                       gCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATg                                                        >  1:38235/1‑150 (MQ=255)
                                                                                       gCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATg                                                        <  2:1146659/150‑1 (MQ=255)
                                                                                        cGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGa                                                       >  1:452696/1‑150 (MQ=255)
                                                                                          aaaTATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCCATGAAATGAcc                                                     <  2:859879/150‑1 (MQ=255)
                                                                                           aaTATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCAGAGGGAATATGTCTTACGTGATTTTCAACCATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCa                                                    >  2:1461682/1‑150 (MQ=255)
                                                                                               tGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGCACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAACCCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCgc                                                >  2:414002/1‑149 (MQ=255)
                                                                                                 cTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAAc                                              >  2:300428/1‑150 (MQ=255)
                                                                                                   gACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCAACTATGAAATGACCAGCGAACTg                                            >  2:1109130/1‑150 (MQ=255)
                                                                                                   gACATCCTGCTGCCTGACTGCACAGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGATTTGAATGTAATACCATCTATGAAATGACCAGCGAACTg                                            >  2:495612/1‑150 (MQ=255)
                                                                                                    aCATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGTGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTgg                                           <  1:726935/150‑1 (MQ=255)
                                                                                                       tCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGc                                                  <  2:932316/140‑1 (MQ=255)
                                                                                                       tCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGc                                                  >  1:932316/1‑140 (MQ=255)
                                                                                                       tCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGc                                          >  1:294520/1‑148 (MQ=255)
                                                                                                       tCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGc                                          <  2:294520/148‑1 (MQ=255)
                                                                                                        ccTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCaaa                                       >  2:1398361/1‑150 (MQ=255)
                                                                                                              gcCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCt                                 >  2:1080844/1‑150 (MQ=255)
                                                                                                               ccTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCtt                                >  2:40140/1‑150 (MQ=255)
                                                                                                               ccTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCtt                                >  1:200966/1‑150 (MQ=255)
                                                                                                                   aCTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCg                            <  1:115282/150‑1 (MQ=255)
                                                                                                                    cTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGt                           >  1:730724/1‑150 (MQ=255)
                                                                                                                      gCACCGCTTCCGAACAGATGGACTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGa                        <  1:338190/150‑1 (MQ=255)
                                                                                                                       cACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGa                        >  2:35173/1‑150 (MQ=255)
                                                                                                                       cACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGa                        <  1:1194381/150‑1 (MQ=255)
                                                                                                                       cACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGa                        <  1:5330/150‑1 (MQ=255)
                                                                                                                         ccGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGc                      >  1:1438758/1‑150 (MQ=255)
                                                                                                                                 gAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGCCCATCTATGAAATGACCAGCGAACTGGCAAACCGTCTTGGCGGTGAGCAACAGttt              >  1:1532515/1‑150 (MQ=255)
                                                                                                                                 gAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCtt                                <  2:1435794/132‑1 (MQ=255)
                                                                                                                                 gAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCtt                                >  1:1435794/1‑132 (MQ=255)
                                                                                                                                        tGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAg       <  1:1428400/150‑1 (MQ=255)
                                                                                                                                          gACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGc     <  1:1153583/150‑1 (MQ=255)
                                                                                                                                          gACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGc     <  1:414002/150‑1 (MQ=255)
                                                                                                                                           aCTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGcc    >  1:837405/1‑150 (MQ=255)
                                                                                                                                             ttcccgcAGGATGCATCCAGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGCCCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGt  >  2:540540/5‑150 (MQ=255)
                                                                                                                                              tcccgcAGGATGCATCCAGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGc     >  2:357209/4‑146 (MQ=255)
                                                                                                                                              tcccgcAGGATGCATCCAGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGc     <  1:357209/143‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TCTGAAATCAACCGTACCCATGAAATCCTTCAGGATGATAAGAAGTGCGAGCTGATTGTGGTTATCGACTGCCACATGACCTCATCGGCGAAATATGCTGACATCCTGCTGCCTGACTGCACCGCTTCCGAACAGATGGACTTTGCGCTGGATGCATCCTGCGGGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGT  >  NC_000913/942429‑942719

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