Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,207,576 A→G 5.2% intergenic (‑44/‑571) purH ← / → rrsE IMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase/16S ribosomal RNA of rrnE operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,207,5760AG 5.2% 348.4 / 7.0 117intergenic (‑44/‑571)purH/rrsEIMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase/16S ribosomal RNA of rrnE operon
Reads supporting (aligned to +/- strand):  ref base A (59/50);  new base G (4/2);  total (64/53)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.16e-01

Only 100 of 118 total aligned reads displayed.

GCAGCTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAAAATACAG  >  NC_000913/4207430‑4207725
                                                                                                                                                  |                                                                                                                                                     
gcagcTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAAttt                                                                                                                                                    >  2:1091246/1‑150 (MQ=255)
 cagcTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGAATATTACAGAGAGCGTTAGCTGAAtttt                                                                                                                                                   >  2:1139733/1‑150 (MQ=255)
  agcTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAAttttt                                                                                                                                                  <  1:868552/150‑1 (MQ=255)
   gcTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTc                                                                                                                                                 >  1:1061/1‑150 (MQ=255)
     tCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTcgc                                                                                                                                               <  1:483884/150‑1 (MQ=255)
     cccACCCCGCGTGCGGAAAGTGCCTGGGCAAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTcgc                                                                                                                                               <  2:1396425/149‑1 (MQ=255)
      ccACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTcgcg                                                                                                                                              >  1:353185/1‑150 (MQ=255)
        acacCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGaa                                                                                                                                            >  1:381657/1‑150 (MQ=255)
        acacCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGaa                                                                                                                                            <  2:1269961/150‑1 (MQ=255)
          acCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGaaaa                                                                                                                                          <  1:341705/150‑1 (MQ=255)
           ccGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACTCTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTTCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAGTTTTTCGCGaaaat                                                                                                                                         <  1:274287/150‑2 (MQ=255)
            cgcgTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAAc                                                                                                                                        <  2:707597/150‑1 (MQ=255)
            cgcgTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAAc                                                                                                                                        >  2:554206/1‑150 (MQ=255)
            cgcgTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAAc                                                                                                                                        >  2:546878/1‑150 (MQ=255)
                  cGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCt                                                                                                                                  <  1:163125/150‑1 (MQ=255)
                       aGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAGTTTTTCGCGAAAAATTCAGCTAACGc                                                                                                                             <  1:743622/150‑1 (MQ=255)
                          gCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGATTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGcccc                                                                                                                          >  2:1467625/1‑150 (MQ=255)
                              gggCGAAATCGTCGATAGCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATCGTAAATCCCCTGGAGTAGACTATTACAGAGAGCGTTACCTCAAGTATTCGCGAACCCCTCAGCTACCCCCCCTAAc                                                                                                                      >  2:439291/1‑150 (MQ=255)
                               ggCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACg                                                                                                                     >  2:520034/1‑150 (MQ=255)
                               ggCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACg                                                                                                                     >  1:376459/1‑150 (MQ=255)
                                 cGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACggg                                                                                                                   >  2:1861800/1‑150 (MQ=255)
                                    aTTCGACGATACCGGCTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCAt                                                                                                               >  2:585570/1‑150 (MQ=255)
                                      tCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGGTGGTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCTAAAAACTCAGCTAACGCCCCTAACGGGGCATc                                                                                                              <  2:1265067/150‑1 (MQ=255)
                                      gcgacgaACCCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGTCTGGAGGACGGTGTTGCATGGTAAATCCCCTTGATTTGACTATTACAGTGAGCGTAAGCTGAATTTTTCGCGAAAAACTCACCTAACGCCCCTATCGGGGCATc                                                                                                              <  2:526060/149‑1 (MQ=255)
                                       cgacgaTACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATcc                                                                                                             <  1:1688261/150‑1 (MQ=255)
                                       cgacgaTACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCt                                                                                                            >  2:556495/1‑150 (MQ=255)
                                        gacgaTACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCt                                                                                                            >  1:782910/1‑150 (MQ=255)
                                          cgaTACGGTATTTCTCACAAACACTCCGCCGTGCGCGGCGGACTGGACGACGTTGTTTCATGGTAAATCCCCTGGCTTTGACTATTACAGAGAGCGATAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTa                                                                                                          <  2:1644062/150‑1 (MQ=255)
                                          cgaTACCGGCTTTGTCAGAACCCCTGAGCCGAGCGCGGCGGACAGGCCCACGTCGTTGCATGGTAACTCCCCTGGATTTGACTATAACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTa                                                                                                          <  2:105169/150‑1 (MQ=255)
                                          cgaTACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTa                                                                                                          <  2:631046/150‑1 (MQ=255)
                                          cgaTACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTa                                                                                                          <  2:1280957/150‑1 (MQ=255)
                                            aTACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTAtt                                                                                                        >  2:1038452/1‑150 (MQ=255)
                                            aTACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTAtt                                                                                                        <  1:1002264/150‑1 (MQ=255)
                                             tACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTAttt                                                                                                       >  1:1003419/1‑150 (MQ=255)
                                             aaCCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTAttt                                                                                                       <  2:1170183/149‑1 (MQ=255)
                                              ttcGGCTTTGTCAGAAACACTGTGCAGAGCGCGGCGGACGGGACGACGTTGTTGCATGGTAAATCCCGTGGATTTGACTATTACAGAGTGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTATCGGGGCATCCTTAtttt                                                                                                      <  1:921614/148‑1 (MQ=255)
                                               ccGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTAttttt                                                                                                     >  2:580035/1‑150 (MQ=255)
                                               ccGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTGGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTAttttt                                                                                                     <  2:1003554/150‑1 (MQ=255)
                                                 ggCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCg                                                                                                   <  1:1029889/150‑1 (MQ=255)
                                                       gTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCa                                                                                             >  2:801756/1‑150 (MQ=255)
                                                        tCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCAt                                                                                            >  1:667900/1‑150 (MQ=255)
                                                           gAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGt                                                                                         >  2:1515700/1‑150 (MQ=255)
                                                           gAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGt                                                                                         >  2:18769/1‑150 (MQ=255)
                                                           gAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGt                                                                                         >  1:1432463/1‑150 (MQ=255)
                                                               cacTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCGGCATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCAGCCTGATTTTTCGCACGCATTGTTACg                                                                                     >  2:314918/1‑150 (MQ=255)
                                                                 cTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGaa                                                                                   <  1:626532/150‑1 (MQ=255)
                                                                  tGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCGTGCTTATTTTTCGCCCGCATTGTAACGaaa                                                                                  >  1:1280555/1‑150 (MQ=255)
                                                                  tGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAAAGaaa                                                                                  >  1:322739/1‑150 (MQ=255)
                                                                       agagCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGtt                                                                             <  2:728436/150‑1 (MQ=255)
                                                                         cgcgcGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTg                                                                           <  2:407993/149‑1 (MQ=255)
                                                                           cgcgGCGGACTGGACGACGTTGTTGCATGTTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTgcg                                                                         <  1:858825/150‑1 (MQ=255)
                                                                           cgcgGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTgcg                                                                         >  2:1496762/1‑150 (MQ=255)
                                                                             cggcggACTGGACGCCGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTgcgcc                                                                       >  2:487890/1‑149 (MQ=255)
                                                                                  gACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCt                                                                  >  1:634262/1‑150 (MQ=255)
                                                                                   aCTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTc                                                                 >  2:688867/1‑150 (MQ=255)
                                                                                    cTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTcg                                                                <  2:444704/150‑1 (MQ=255)
                                                                                      ggACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTcgcg                                                              >  1:1497573/1‑150 (MQ=255)
                                                                                      ggACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTcgcg                                                              >  2:654989/1‑150 (MQ=255)
                                                                                         cgacgTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTAATTTTCGCCCGCATTGTAACGAAAACGTTTGAGCAACGCTCGCGAAt                                                           >  1:658144/1‑150 (MQ=255)
                                                                                             gttgttGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAAttttt                                                       <  2:1127551/150‑1 (MQ=255)
                                                                                              ttgttgCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAAtttttg                                                      >  2:861809/1‑149 (MQ=255)
                                                                                                gtttCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTtctc                                                    <  1:1610311/150‑1 (MQ=255)
                                                                                                gttgCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCTAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTtctc                                                    <  2:948370/150‑1 (MQ=255)
                                                                                                gttgCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTtctc                                                    <  2:201622/150‑1 (MQ=255)
                                                                                                 ttgCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTctct                                                   >  2:1752310/1‑150 (MQ=255)
                                                                                                  tgCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCAGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCtt                                                  >  1:55144/1‑150 (MQ=255)
                                                                                                   gCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTATCGAATACGTTTGCGCACCGCTCGCGAAGTTTTCTCttt                                                 >  2:1304498/1‑150 (MQ=255)
                                                                                                   gCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGAAACGCTCGCGAATTTTTCTCttt                                                 >  2:757322/1‑150 (MQ=255)
                                                                                                       ggTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAAt                                             >  2:1828424/1‑150 (MQ=255)
                                                                                                        gTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATg                                            <  2:24729/150‑1 (MQ=255)
                                                                                                         tAAATCCCCTGGATTTGACTATTACAGAGAGGGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATgg                                           <  1:554206/150‑1 (MQ=255)
                                                                                                         aaaaTCCCCTGGATTTTACTATTACAGACAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATgg                                           <  1:546878/149‑1 (MQ=255)
                                                                                                          aaaTCCCCTGTATTTGACTATTACAGAGAGCGTTAGCTGAAATTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGt                                          <  2:1741306/150‑1 (MQ=255)
                                                                                                           aaTCCCCTGGATTTGACTATTACAGATAGCTTTAGCTGAATTTTTCGCTAAAATCTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTg                                         <  2:1522905/150‑1 (MQ=255)
                                                                                                           aaTCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTg                                         <  1:590431/150‑1 (MQ=255)
                                                                                                            aTCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTTAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGTGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGa                                        <  1:1082847/150‑1 (MQ=255)
                                                                                                               cccTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATca                                     >  1:381006/1‑150 (MQ=255)
                                                                                                                ccTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATcac                                    >  1:1873435/1‑150 (MQ=255)
                                                                                                                ccTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATcac                                    >  2:1008508/1‑150 (MQ=255)
                                                                                                                      tttGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAAtttt                              <  1:1297974/150‑1 (MQ=255)
                                                                                                                       ttGACTATTACAGAGAGCGTTAGCTGACTTTTTCGCGAAAAACTCCGATAACGCCCCTCACGTGGCAGCCTTATTTTTCGCCCGCATAGTAACGAACACGTTTGCGCCACGCTCGCGAACTTTTCTCTTTCAATGGTGATCACAATTTTg                             <  2:85606/150‑1 (MQ=255)
                                                                                                                         gACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGAc                           <  1:328745/150‑1 (MQ=255)
                                                                                                                          aCTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCATCTAACGCCCCTAACGGGGCATGCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACt                          <  2:676818/150‑1 (MQ=255)
                                                                                                                           cTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACtg                         <  1:1552235/150‑1 (MQ=255)
                                                                                                                           cTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTATTCGCCCGCATTGTAACGAAAACGTTTGAGCAAAGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACtg                         >  1:1774358/1‑150 (MQ=255)
                                                                                                                                aCAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGTCATCCTTATTTTTCGCCCGCATTGTAACGAACACGTTTGCGCAACGCTCGCGAATTTTTCGCATTCAATGGTGATCACAATTTTGACTGTGGtt                    >  1:288849/1‑150 (MQ=255)
                                                                                                                                aCAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCAGCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTATTTCAATGGTGATCACAATTTTGACTGTGGtt                    >  2:143673/1‑150 (MQ=255)
                                                                                                                                  agagagCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTGTCAATGGTGAGCACAATTTTGACTGTGGTTAc                  >  1:1481990/1‑150 (MQ=255)
                                                                                                                                    ggggCGTTAGCTGAGTTTTTCGCGAAAAATTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCg                >  1:78715/4‑150 (MQ=255)
                                                                                                                                    agagCGGTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCg                <  1:1139733/150‑1 (MQ=255)
                                                                                                                                         gTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAAATTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGc           <  2:528693/150‑1 (MQ=255)
                                                                                                                                           tAGCTGAGTTTTTCGCGAAAAATTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCaa         >  1:500597/1‑150 (MQ=255)
                                                                                                                                           tAGCTGAGTTTTTCGCGAAAAATTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCaa         >  1:973767/1‑150 (MQ=255)
                                                                                                                                            aGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCaaa        >  2:1105672/1‑150 (MQ=255)
                                                                                                                                            aGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCaaa        <  2:485991/150‑1 (MQ=255)
                                                                                                                                             gCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCaaaa       >  1:46227/1‑150 (MQ=255)
                                                                                                                                              cTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGATTTTTTCTCTTTCATTGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAAAAt      >  1:1769298/1‑150 (MQ=255)
                                                                                                                                               tGAATTTTTCGCTAAAAACTCAGCTAACGCCCCTAACGGGGCAGCCTGATTTTTCGCCCGCATTGAAACGAAAACGTTTGCGCAACGCTCGCGAATTATTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAAAATa     <  1:331496/150‑1 (MQ=255)
                                                                                                                                                gAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCTCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGGGGTTACCGGGGGCAAAATAc    >  2:1117625/1‑150 (MQ=255)
                                                                                                                                                  aTTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTTTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAAAATACAg  <  2:33974/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCAGCTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAAAATACAG  >  NC_000913/4207430‑4207725

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: