Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,707,031 C→A 5.8% P38H (CCC→CAC)  rsxC → SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,707,0310CA 5.8% 287.1 / 4.6 102P38H (CCC→CAC) rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (42/52);  new base A (4/3);  total (46/56)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

ACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTA  >  NC_000913/1706883‑1707180
                                                                                                                                                    |                                                                                                                                                     
acacCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTAccc                                                                                                                                                      >  1:1590615/1‑150 (MQ=255)
 cacCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTAcccc                                                                                                                                                     <  2:1162917/150‑1 (MQ=255)
   ccATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTg                                                                                                                                                   >  1:572406/1‑150 (MQ=255)
   ccATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTg                                                                                                                                                   >  2:1085687/1‑150 (MQ=255)
      ttCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGgcc                                                                                                                                                >  1:703514/1‑149 (MQ=255)
        cccGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGTACCTGGCGCa                                                                                                                                              <  2:1520019/150‑1 (MQ=255)
         ccGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAg                                                                                                                                             <  1:1376171/150‑1 (MQ=255)
         ccGTGCGTATCATTCCCGTCGAACCCCATGCTTAAGATATCCTCAGCATTCAGAAAAAATAAAATCTCGGCTTTCAACGGCGGCATCCATCCACCGGAGATGAACACCCAGTCCAACGGAACAGCCCAGCGCCAGGGGTACCTGGCGCAg                                                                                                                                             <  1:1610949/150‑1 (MQ=255)
           gTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCg                                                                                                                                           <  1:891306/150‑1 (MQ=255)
            tGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTACGCCAGGGGTACCTGGCGCAgggg                                                                                                                                          <  2:1808945/150‑4 (MQ=255)
             gCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAAAGGCTGCATACATCCACCGGAGATGAAAACCAAGTCAAACGGTACACCCCTGCGCCAGGTACCCAGGGCGCAGCGtt                                                                                                                                         >  1:147424/1‑150 (MQ=255)
               gTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCCTCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGtttt                                                                                                                                       <  2:697243/150‑1 (MQ=255)
               gTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCCCCACGTACCCCTGGCGCAGCggtgt                                                                                                                                       >  1:1665748/1‑146 (MQ=255)
                tATCATTCCCGTGTAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTg                                                                                                                                      <  2:1450610/150‑1 (MQ=255)
                tATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAAACCAGTCAAAAGGTACACCCCTGCGCCAGGTACCCCTGTCGCACCGTTTTg                                                                                                                                      >  1:582415/1‑150 (MQ=255)
                tATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCAGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTg                                                                                                                                      <  1:126710/150‑1 (MQ=255)
                 aTCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGt                                                                                                                                     <  2:669741/150‑1 (MQ=255)
                   cATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTa                                                                                                                                   >  2:1905776/1‑150 (MQ=255)
                        ccGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCa                                                                                                                              >  1:711919/1‑150 (MQ=255)
                         cGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCAc                                                                                                                             >  1:513534/1‑150 (MQ=255)
                           tgtaaCACCATGCATAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTg                                                                                                                           <  2:1419266/147‑1 (MQ=255)
                           tGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCACCGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTg                                                                                                                           >  2:738338/1‑150 (MQ=255)
                            ggAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGa                                                                                                                          >  1:1868216/1‑150 (MQ=255)
                              aaCACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGaaa                                                                                                                        <  1:676007/150‑1 (MQ=255)
                                 acCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAg                                                                                                                     <  2:1034860/150‑1 (MQ=255)
                                 acCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTATTCCACTGAAACAGc                                                                                                                    <  2:969954/150‑1 (MQ=255)
                                  ccATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGc                                                                                                                    <  2:1004111/150‑1 (MQ=255)
                                  ccATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACAGAAACAGc                                                                                                                    >  1:1683155/1‑150 (MQ=255)
                                  ccATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCAGGCGCAGCGGTTTGTTATTACACTGAAACAGc                                                                                                                    >  2:797046/1‑150 (MQ=255)
                                      gCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACACCTGGCGCAGCGTTTTGTTATTACACTGAAACACCatat                                                                                                                >  1:190302/1‑150 (MQ=255)
                                       cTTCAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATAtt                                                                                                               <  2:29202/150‑1 (MQ=255)
                                         tAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGTGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTgg                                                                                                             <  2:610544/150‑1 (MQ=255)
                                              tATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTg                                                                                                        >  1:1779741/1‑150 (MQ=255)
                                              tATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCACGTACCCCAGGCGCACCGTATTGTTATTCCACTGAAACAGCATAATGGCGCTg                                                                                                        >  1:1852567/1‑150 (MQ=255)
                                                 tctcTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAg                                                                                                     <  1:39223/150‑1 (MQ=255)
                                                  ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAgg                                                                                                    >  2:1822507/1‑150 (MQ=255)
                                                  ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGGTTTGTTTTTGCACTGAAAGAGCAGATTGGCGCTGTAgg                                                                                                    >  2:1396556/1‑150 (MQ=255)
                                                   tctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGt                                                                                                   <  2:1655946/150‑1 (MQ=255)
                                                    ctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTg                                                                                                  >  1:514831/1‑150 (MQ=255)
                                                     tGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCCGCGGTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGa                                                                                                 >  1:761632/1‑150 (MQ=255)
                                                      gCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAg                                                                                                <  1:1802296/150‑1 (MQ=255)
                                                      gCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGACCCCGGGCGCAGCGTGTTGTTATTCCACTGAAGCAGCAGATTGGCGCGGCAggggcg                                                                                                >  1:307548/1‑146 (MQ=255)
                                                           cAGAAAAAATAACTTCGGTGATTTCAACCTCGGCATCCAGCCACCGGAGATGAAAACCCAGTCCAACGTTACATCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgtg                                                                                           <  2:1236970/150‑1 (MQ=255)
                                                           cAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGGtgtg                                                                                           >  1:785654/1‑150 (MQ=255)
                                                            aGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGc                                                                                          <  2:1323123/150‑1 (MQ=255)
                                                             taaaaaaTAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCg                                                                                         <  2:147424/149‑1 (MQ=255)
                                                             gAATAAATAAAATCTGGTATTTGTACGGCGGCATCCATCCACCGGAGATGATATCCGAGTCCTACGGTACACCCCTTCGCCAGGTACCCCTGGCGCAGCGTTTTGGTATTCCTCTGAAACAGCATATTGGCGCTGATGGTGAGTTGTGCg                                                                                         <  1:64065/150‑1 (MQ=255)
                                                             gAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCg                                                                                         >  2:35858/1‑150 (MQ=255)
                                                              aaaaaaTAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGATGAAGGTGAGTTGTGCGt                                                                                        <  1:1810463/150‑1 (MQ=255)
                                                                aaaataaaatCTGGGAATTCAACGGCGGCATCCCTCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTa                                                                                      <  2:46478/150‑1 (MQ=255)
                                                                 aaataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGATAACCCAGTCCAACGGTACCCCCCTGCGCCAGGTACGCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAg                                                                                     <  2:1819276/150‑1 (MQ=255)
                                                                 aaataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAg                                                                                     <  1:621653/150‑1 (MQ=255)
                                                                  aataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGc                                                                                    <  1:1483502/150‑1 (MQ=255)
                                                                  aataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGACCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGATGTGCGTTAGc                                                                                    >  2:806330/1‑150 (MQ=255)
                                                                     aaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTc                                                                                 >  2:464446/1‑150 (MQ=255)
                                                                        atCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGc                                                                              <  1:1494700/150‑1 (MQ=255)
                                                                        atCTGGGAATTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGGTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGc                                                                              <  2:1089373/150‑1 (MQ=255)
                                                                          cTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGa                                                                            <  2:172376/150‑1 (MQ=255)
                                                                           tGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGAt                                                                           >  2:1079601/1‑150 (MQ=255)
                                                                           tGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCCCCCGGTACACCGGGCGGACCGGTTTGTTATTACCCTGATACAGCATTTTGGCGCTGTAGGTGAGTTGTGGGGTTGCGTCGGAGAt                                                                           >  1:1597490/1‑150 (MQ=255)
                                                                              gATTTCAACGGCGGCATCCATCCACCGGAGAAGAACACCCAGTCCAACGGTACACCCCTGCGCCAGGAACCCCTGGCGCAGCGTTTTGTTATTCCACTGCAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGAGaaa                                                                        >  1:296188/1‑150 (MQ=255)
                                                                               aTTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAg                                                                       <  2:514025/150‑1 (MQ=255)
                                                                                tttCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGt                                                                      >  2:878633/1‑150 (MQ=255)
                                                                                 ttCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCCGGTACCCCAGGCGCAGCGGTTTGTTTTTCCACAGAAACAGCATATTGGCGGTGAAGGAGAGTTGGGCGGTAGCGTCGGCGAGAAAgga                                                                     >  2:1879347/1‑148 (MQ=255)
                                                                                 ttCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTa                                                                     <  2:763500/150‑1 (MQ=255)
                                                                                   cAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCCGCGTTTTGTTATTCCACTGAAACTGCATATTGGCGCTGAAGGTGAGTTGTGGGTTAGCGTCGGCGATAAAGTAtt                                                                   >  2:311811/1‑150 (MQ=255)
                                                                                       ggcggcATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTgcgc                                                               <  2:585650/150‑1 (MQ=255)
                                                                                         cggcATCCATCCACCGGAGCTGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGGTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCgg                                                             <  2:8318/150‑1 (MQ=255)
                                                                                         cggcAACCAACCACCGGAGATGAAAACCCAGTCCAACGGTACCCCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGCTATTCCACTGAAACACCATATTGGCGCTGAAGGTGAGTTGTGCGGTAGCGTCGGCGATAAAGTATTGCGCgg                                                             <  2:1176131/150‑1 (MQ=255)
                                                                                            catccatcCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGACGATAAAGTATTGCGCGGCCa                                                          <  2:803030/150‑1 (MQ=255)
                                                                                              tccatccaCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGc                                                        >  2:23006/1‑150 (MQ=255)
                                                                                               ccatccaCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGcc                                                       <  2:449358/150‑1 (MQ=255)
                                                                                                catccaCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTGTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTAGTGCGTAAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCg                                                      >  1:85989/1‑150 (MQ=255)
                                                                                                    caCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTa                                                  >  1:191423/1‑150 (MQ=255)
                                                                                                    caCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCATa                                                  >  2:1232596/1‑150 (MQ=255)
                                                                                                      ccGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTAcc                                                <  1:755711/150‑1 (MQ=255)
                                                                                                      ccGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTAcc                                                <  1:1905776/150‑1 (MQ=255)
                                                                                                        ggAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCg                                              >  2:704985/1‑150 (MQ=255)
                                                                                                        ggAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCACTCGCGCAGCGTTTTGTTATTCCGCTGAAACAGAATATTGTCGCAGAAGGAGAGTTGTGCGTTAGCGACGGCGATAAAGTAATGCGCGGCCGGCCGCTTACCCg                                              >  2:1342671/1‑150 (MQ=255)
                                                                                                            aTGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTCTTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGt                                          <  2:92540/150‑1 (MQ=255)
                                                                                                               acaaCCCAGTCCAACGCTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTcgc                                       <  1:315487/148‑1 (MQ=255)
                                                                                                                 aaCCCAGTCCAACGGTACACCCCTGCCCCAGGAACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCgg                                     <  1:1707697/150‑1 (MQ=255)
                                                                                                                     cAGTCCACCGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCaaa                                 <  2:719353/150‑1 (MQ=255)
                                                                                                                     cAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCaaa                                 <  2:1886525/150‑1 (MQ=255)
                                                                                                                           aaCGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGGTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGGCGCGGCAAAAtgctg                           <  2:1293862/150‑1 (MQ=255)
                                                                                                                               gTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGGTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTg                       >  1:1773385/1‑150 (MQ=255)
                                                                                                                                tACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACGGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGt                      >  2:261174/1‑150 (MQ=255)
                                                                                                                                  cacCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTGGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTc                    <  1:1085687/150‑1 (MQ=255)
                                                                                                                                    ccccTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCATAATGCTGCCTGTTCAc                  >  1:70531/1‑150 (MQ=255)
                                                                                                                                     cgcTGCGCCAGGGGTACCTGGCGCACCGTTTTGTTAATCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCCTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGAGGTAGCGGCAATATGCTGACTTTTCAcc                 >  2:1460767/3‑149 (MQ=255)
                                                                                                                                     cccTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAATCAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAcg                 <  1:1162481/150‑1 (MQ=255)
                                                                                                                                     cccTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTAATCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAcg                 <  2:1256859/150‑1 (MQ=255)
                                                                                                                                      gcTGCGCCAGGGGTACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAcgc                >  1:224494/2‑150 (MQ=255)
                                                                                                                                       cTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAcgcg               <  2:870974/150‑1 (MQ=255)
                                                                                                                                               ggTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCAccgc       >  2:103282/1‑148 (MQ=255)
                                                                                                                                                  aCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCggg    <  1:892961/150‑1 (MQ=255)
                                                                                                                                                  aCCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCAGTTCACGCGCACACCTCggg    >  2:790939/1‑150 (MQ=255)
                                                                                                                                                   ccccTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGt   <  2:190302/150‑1 (MQ=255)
                                                                                                                                                    cccTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTa  >  1:519911/1‑150 (MQ=255)
                                                                                                                                                    cccTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTa  >  1:709417/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTA  >  NC_000913/1706883‑1707180

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