Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,580 G→A 8.1% intergenic (+25/+216) ybiT → / ← ybiU ABC‑F family putative regulatory ATPase/DUF1479 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913857,5800GA 8.1% 287.8 / 9.7 99intergenic (+25/+216)ybiT/ybiUABC‑F family putative regulatory ATPase/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base G (52/37);  new base A (4/4);  total (56/43)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

Only 100 of 102 total aligned reads displayed.

TTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATT  >  NC_000913/857431‑857729
                                                                                                                                                     |                                                                                                                                                      
ttGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTg                                                                                                                                                        >  2:1521249/1‑150 (MQ=255)
ttGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTg                                                                                                                                                        >  1:833804/1‑150 (MQ=255)
     tCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑ggg                                                                                                                                                  <  1:490509/150‑1 (MQ=255)
     tCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACACCGATACTGTC‑ggg                                                                                                                                                  >  1:1446402/1‑150 (MQ=255)
     tCACACGACCGTGAGTTCGTAAGCTCCCTGCCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑ggg                                                                                                                                                  <  1:1147628/150‑1 (MQ=255)
       acacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGCGt                                                                                                                                                >  1:1225083/1‑150 (MQ=255)
       acacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATC‑GGGGt                                                                                                                                                <  2:651159/150‑1 (MQ=255)
       acacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGAGCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGt                                                                                                                                                >  1:1148607/1‑150 (MQ=255)
        cacGACCGTGAGTGAGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGTAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATC‑GGGGTg                                                                                                                                               <  1:697470/150‑1 (MQ=255)
         acGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGa                                                                                                                                              >  1:1402325/1‑150 (MQ=255)
           gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATC‑GGGGTGAcg                                                                                                                                            <  1:1619522/150‑1 (MQ=255)
                tGATTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCCCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGTATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACggcggc                                                                                                                                       <  2:801888/150‑1 (MQ=255)
                      cGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTATAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGA                                                                                                                                 >  1:1039732‑M1/1‑150 (MQ=255)
                      cGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATC‑GGGGt                                                                                                                                                >  2:474778/1‑135 (MQ=255)
                      cGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATC‑GGGGt                                                                                                                                                <  1:474778/135‑1 (MQ=255)
                       gTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGAC                                                                                                                                >  1:1449106‑M1/1‑150 (MQ=255)
                       gTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGAC                                                                                                                                >  2:938440‑M1/1‑150 (MQ=255)
                        tAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACA                                                                                                                               >  2:940744‑M1/1‑150 (MQ=255)
                           gCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAAT                                                                                                                            >  1:1640540‑M1/1‑150 (MQ=255)
                            cTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATG                                                                                                                           <  2:1002413‑M1/150‑1 (MQ=255)
                              cccTGGCGCCCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACATGCGAAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCACGATAATGTC‑GGGGTGACGGCGGCTGATGACAAATGCa                                                                                                                         >  1:989530/1‑150 (MQ=255)
                               ccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCTTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAAataaa                                                                                                                        >  1:301401‑M1/1‑145 (MQ=255)
                                cTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCtaa                                                                                                                       >  1:995438‑M1/1‑147 (MQ=255)
                                       cccGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCt                                                                                                                >  2:184136/1‑150 (MQ=255)
                                           cATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATg                                                                                                            >  2:779207/1‑150 (MQ=255)
                                            aTTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgc                                                                                                           <  1:224374/150‑1 (MQ=255)
                                                 ggAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTa                                                                                                       >  2:209698/1‑150 (MQ=255)
                                                 ggAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGAGCTAc                                                                                                      >  1:306015/1‑150 (MQ=255)
                                                 ggAAATCACCCCGGAACGCGTGATCGACTTGAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTAc                                                                                                      >  1:688646/1‑150 (MQ=255)
                                                        aCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATc                                                                                               <  1:38898/150‑1 (MQ=255)
                                                            cGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGc                                                                                           <  2:1047053/150‑1 (MQ=255)
                                                                  gcgTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGc                                                                                     <  1:47508/150‑1 (MQ=255)
                                                                   cgTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTTCGCa                                                                                    >  1:209853/1‑150 (MQ=255)
                                                                      gATCGACTTTCGCGGTAATTCCTAATATTACGTGCGTAGTACTGGGAGCGAGTAAGTGATTATCCGTCCCCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTa                                                                                 <  1:574530/150‑1 (MQ=37)
                                                                         cGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCACTACGCTTATCAGGCCTACGCAGTACCt                                                                              >  1:496878/1‑150 (MQ=255)
                                                                          gACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGCCGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTg                                                                             >  1:501317/1‑150 (MQ=255)
                                                                           aCTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGc                                                                            >  1:1230607/1‑150 (MQ=255)
                                                                           aCTTAAGCGTTAATAACGATTATTACCTGCGTAGTAAAGGGATGTATTAAGTGTTTATCCGTCACCCCGATACTTTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGc                                                                            <  1:1389080/150‑1 (MQ=37)
                                                                                   ggTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt                                                                    >  1:334396/1‑150 (MQ=255)
                                                                                   ggTAATTACGAAGATTACATGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt                                                                    <  2:739400/150‑1 (MQ=255)
                                                                                      aaTTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGaa                                                                 >  2:188524/1‑150 (MQ=255)
                                                                                          aCGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTg                                                             >  1:1370026/1‑150 (MQ=255)
                                                                                          aCGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGAAATATATTGATTCTg                                                             >  2:1349786/1‑150 (MQ=255)
                                                                                           cGAACATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGGCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGc                                                            <  1:729114/150‑1 (MQ=255)
                                                                                            gAAGATTACAGGCGTAGTACAGTGATCCAGTCAGTGTATATCCGTCACCCCTAAACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCATATCAGGCCTACTCAGTACCTGCAATATATTGAATCTGCa                                                           <  1:274960/150‑1 (MQ=255)
                                                                                             aaGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCAGATCCTCTACGCTTATCAGGCCTACGCATTACCTGCAATATATTGAATCTGCat                                                          >  1:129369/1‑149 (MQ=255)
                                                                                              tttttACCTGGGTTTAAACGGCATCGCGTATGTTATTAGCGGCCACCCCAAAACTTGG‑GGGGTGAGGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAgg                                                         <  1:1343999/147‑1 (MQ=18)
                                                                                                tttACCTGCGTAGTATAGGGATCGTGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGAc                                                       <  2:102725/149‑1 (MQ=255)
                                                                                                ataacttGCGTAGTAATGGGATCGAGTATGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGAc                                                       <  2:818143/144‑1 (MQ=255)
                                                                                                aTTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGAc                                                       <  2:1208291/150‑1 (MQ=255)
                                                                                                aTTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGAc                                                       >  2:463213/1‑150 (MQ=255)
                                                                                                   aCCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACttt                                                    <  2:833804/150‑1 (MQ=255)
                                                                                                     cTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGt                                                  >  2:905395/1‑150 (MQ=255)
                                                                                                      tGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTa                                                 >  2:277242/1‑150 (MQ=255)
                                                                                                        cGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCgctttgttaatgcctgatgcggcgttaacgccttatccggcctacaaagtcctgcagattcaatatattgcaggt                                               >  1:241784‑M1/1‑75 (MQ=255)
                                                                                                         gtagtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACATGCAATATATTTAATCTGCAGGACTTTGTAGGc                                              >  2:311163/1‑150 (MQ=255)
                                                                                                          tagtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGcc                                             >  2:346886/1‑150 (MQ=255)
                                                                                                           tgtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCg                                            <  2:1112904/149‑1 (MQ=255)
                                                                                                           agtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCg                                            >  1:1330777/1‑150 (MQ=255)
                                                                                                            ttaAAGGTATCTAGTAAGGTTTTATCCGTAACCCCGATACTGTG‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCgg                                           <  2:1394431/149‑1 (MQ=25)
                                                                                                              aaaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGAt                                         >  1:88306/1‑150 (MQ=255)
                                                                                                               aaGGGATCGAGTAAGTGATTATCCGTCACCCCGAAACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCTGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATa                                        >  1:1303543/1‑150 (MQ=255)
                                                                                                                 gggATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAg                                      >  1:1530281/1‑150 (MQ=255)
                                                                                                                  ggATCGAGTAAGTGATTATCGGTCACCGCGATAGTGTC‑GGGGTGACGGCGGCTGATGACAAATGAAAAATAGCCTGATGCGCTACGCTTATCAGGCCTAGGCAGTACCTGCAATATTTTGAATCTGCAGGACTTTGTAGGCCGGATAAgg                                     <  2:920717/150‑1 (MQ=37)
                                                                                                                   gATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                                    >  1:408117/1‑150 (MQ=255)
                                                                                                                   gATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                                    >  2:1256851/1‑150 (MQ=255)
                                                                                                                    agcGAGTAAGTGATTATCCGTCACCCCGAAACTGCC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGTGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCg                                   <  1:851605/148‑1 (MQ=255)
                                                                                                                    aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCg                                   >  1:1594113/1‑150 (MQ=255)
                                                                                                                     tCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGt                                  <  2:1225083/150‑1 (MQ=255)
                                                                                                                          tAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcg                             <  2:1640540/150‑1 (MQ=255)
                                                                                                                            aGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgcc                           >  2:1131566/1‑150 (MQ=255)
                                                                                                                            aGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgcc                           >  2:1730793/1‑150 (MQ=255)
                                                                                                                             ttGATAATCCTTCACCCCGATCCTGTC‑GGGGTTACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgccg                          <  2:877224/149‑1 (MQ=255)
                                                                                                                             gggATTATCCGACCCCCCGATACTGCC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgccg                          <  1:1521249/148‑1 (MQ=255)
                                                                                                                              tGATTATCCGTCTCCCCGATACTGTC‑TGGGGGACGGCGGCTGATGACAAATTCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgccgc                         <  2:301401/150‑1 (MQ=255)
                                                                                                                               gATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCa                        >  1:791994/1‑150 (MQ=255)
                                                                                                                                aTTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCAt                       <  2:1730944/150‑1 (MQ=255)
                                                                                                                                aTTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGATTATCAGGCCTACGCAGTACCTGCAATATATTTAATCTGCCGGACTTTGTAGGCCGGATAAGTCGTTAACGACGCaa                       >  2:1044908/1‑149 (MQ=38)
                                                                                                                                 ttATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATc                      <  1:724798/150‑1 (MQ=255)
                                                                                                                                 gccgccGTCACCCCGACAGTATC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATTTCCAGATGCGCTACGCTTATCAGGCCTACGCAGTACCAGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATc                      >  2:1645115/5‑150 (MQ=21)
                                                                                                                                  aatcccgcACCCCTAAACTGCC‑GGGGGTACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCAGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCa                     <  1:1349786/143‑1 (MQ=21)
                                                                                                                                   aTCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATTTTGATTCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAg                    >  1:976201/1‑150 (MQ=255)
                                                                                                                                     cggccACCCCTATCCTGCG‑GGGGTTACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGc                  <  2:1536468/146‑1 (MQ=21)
                                                                                                                                     ccGTCACCCCTATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGc                  <  2:1682119/150‑1 (MQ=255)
                                                                                                                                      cGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCa                 >  1:1091774/1‑150 (MQ=255)
                                                                                                                                      cGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCa                 >  1:355011/1‑150 (MQ=255)
                                                                                                                                       gTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTTAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCAGCATCAGGCAt                >  2:1085094/1‑150 (MQ=255)
                                                                                                                                        tCACCCCGATACTGTCGGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAt                <  1:209698/150‑1 (MQ=255)
                                                                                                                                        cccccccGATACTGGC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCTGATGCGGCGTTAACGCCTTATCCGGccta               <  1:514788/147‑4 (MQ=21)
                                                                                                                                         cACCCCGATACTGTT‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa              <  2:892552/150‑1 (MQ=255)
                                                                                                                                         cACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa              >  2:123837/1‑150 (MQ=255)
                                                                                                                                         cACCCCGACAGTATC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa              >  2:1003967/1‑150 (MQ=255)
                                                                                                                                         cACCCCGACAGTATC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa              >  2:365530/1‑150 (MQ=255)
                                                                                                                                          aCCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCgctttgttaatgcctgatgcggcgttaacgccttatccggcctacaaagtcc                                                                      >  2:1496879‑M1/1‑41 (MQ=255)
                                                                                                                                          aCCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCgctttgttaatgcctgatgcggcgttaacgccttatccggcctacaaagtcc                                                                      <  1:1496879‑M1/93‑53 (MQ=255)
                                                                                                                                           ccccGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc            >  2:1450411/1‑150 (MQ=255)
                                                                                                                                           ccccGAAACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc            <  1:856487/150‑1 (MQ=255)
                                                                                                                                            cccGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCAACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACa           >  2:1097079/1‑150 (MQ=255)
                                                                                                                                              cGATACTGTC‑GGGGTGACGGCGGCAGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACaaa         <  1:232199/150‑1 (MQ=255)
                                                                                                                                                     gTC‑GGGGTGACGGCGGCTTATGACAAATTCAAAATAGCCTGATGCGCTACGCTTATCAGTCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCgcatt  >  2:306015‑M2/1‑145 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTC‑GGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATT  >  NC_000913/857431‑857729

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