Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,597,699 G→T 7.4% intergenic (‑81/+132) purC ← / ← bamC phosphoribosylaminoimidazole‑succinocarboxamide synthetase/BamABCDE complex OM biogenesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,597,6990GT 7.4% 177.3 / 2.2 70intergenic (‑81/+132)purC/bamCphosphoribosylaminoimidazole‑succinocarboxamide synthetase/BamABCDE complex OM biogenesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (38/26);  new base T (3/3);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGGATGATTCCAGCCCTGTATTTTTACTTGCTAAA  >  NC_000913/2597554‑2597842
                                                                                                                                                 |                                                                                                                                               
tCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAAc                                                                                                                                              <  2:647641/149‑1 (MQ=255)
     gTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTTAATTCTTTTGGTGTCTGGTGTGCCGGATGTTTTGTAGGATGCGGCGTGAACGCCtt                                                                                                                                         >  1:839571/1‑149 (MQ=255)
     gTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCtt                                                                                                                                         >  2:79412/1‑149 (MQ=255)
      ttttCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTa                                                                                                                                        <  2:381776/149‑1 (MQ=255)
        ttCCGTGGTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATc                                                                                                                                      <  1:277111/149‑1 (MQ=255)
        ttCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTCTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATc                                                                                                                                      <  1:273212/149‑1 (MQ=255)
          ccGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                    <  2:456018/149‑1 (MQ=255)
           cGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGa                                                                                                                                   >  1:342510/1‑149 (MQ=255)
            gTGCTGTATACTGTATTCGCTTTACCACGATACACCTCAGCTTGCTTATGCATCTTTATCTCTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACTCCTTATCCGaa                                                                                                                                  <  1:956062/149‑1 (MQ=255)
                    taCGGCTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTcc                                                                                                                          <  1:584437/149‑1 (MQ=255)
                     aCGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGTTGTGCCGGTTTTTTTGTCGGATGCGGCGTGAACGCCTTTTCCGAACTACGTCCg                                                                                                                         >  1:108174/1‑149 (MQ=255)
                     aCGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTAAGTCCg                                                                                                                         >  2:403103/1‑149 (MQ=255)
                      cctgtttCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTTGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCGGGGGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGACTTATCCGAACTACGCCCGa                                                                                                                        <  1:256011/145‑1 (MQ=255)
                       ggTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGGGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGac                                                                                                                       <  1:984935/149‑1 (MQ=255)
                       ggTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCGTTATCACGCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGGTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGac                                                                                                                       >  2:203149/1‑149 (MQ=255)
                            tCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGCTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACTCCTTAACCGAACTACGTCCGacaaaca                                                                                                                  >  1:932735/1‑149 (MQ=255)
                              gCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGAcaaacaaa                                                                                                                <  2:498358/149‑1 (MQ=255)
                                   aCCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCg                                                                                                           <  1:339706/149‑1 (MQ=255)
                                     cACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTg                                                                                                         >  2:592221/1‑149 (MQ=255)
                                     cACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGACTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTg                                                                                                         >  1:67468/1‑149 (MQ=255)
                                         aTACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGCCAAACAAATTTCGTGCGaa                                                                                                     >  1:231274/1‑149 (MQ=255)
                                             aaCTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTAc                                                                                                 >  1:47228/1‑149 (MQ=255)
                                              aCTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTTAACGCCTTATCCGAACTACGTCCGACAAACAAATTTAGTGCGAATTAcc                                                                                                >  2:121346/1‑149 (MQ=255)
                                              aCTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTAcc                                                                                                >  1:868005/1‑149 (MQ=255)
                                              aCTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTAcc                                                                                                >  1:397271/1‑149 (MQ=255)
                                               cTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTAccc                                                                                               >  2:391527/1‑148 (MQ=255)
                                                tCAGCATGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTGTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGc                                                                                              <  2:839571/149‑1 (MQ=255)
                                                    cttgcttTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGGCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTaaa                                                                                          <  1:969268/149‑1 (MQ=255)
                                                     ttgcttTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACAGCTaaaa                                                                                         >  2:932002/1‑149 (MQ=255)
                                                        cttTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATcg                                                                                      >  2:91772/1‑149 (MQ=255)
                                                         ttttGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTTAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATcgc                                                                                     <  1:200173/149‑1 (MQ=255)
                                                             gCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGttt                                                                                 >  2:649128/1‑149 (MQ=255)
                                                               aTCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTc                                                                               >  2:406407/1‑149 (MQ=255)
                                                               aTCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTc                                                                               >  2:724463/1‑149 (MQ=255)
                                                               aTCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTc                                                                               >  1:10709/1‑149 (MQ=255)
                                                               aTCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTc                                                                               >  2:554235/1‑149 (MQ=255)
                                                                   ttATCACTCCTGGGTGTGAATTACCGTTTTAAACTCGTTTGCTGTCTGGTGTGCCGGAGGTTTTGTCGGATGAGGCGTGAACGCATTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGc                                                                           <  1:178356/149‑1 (MQ=255)
                                                                   ttATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGc                                                                           >  2:901615/1‑149 (MQ=255)
                                                                   ttATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGc                                                                           >  1:230742/1‑149 (MQ=255)
                                                                         cTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGAGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACAGCTAAAATCGCGTTATCCTGCAGACGc                                                                     >  1:128152/1‑149 (MQ=255)
                                                                                 tgtgAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGttt                                                             <  2:293587/149‑1 (MQ=255)
                                                                                          aCGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGtatcat                                                    >  2:587999/1‑149 (MQ=255)
                                                                                            gctttAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGtatcatat                                                  <  2:858095/147‑1 (MQ=255)
                                                                                            gTTTTAAAATCTTATGCTGTCTGGTGTGCCGGATGTATTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGAGCGAATTACCGCTAAAATCGCGTGTTCCTGCAGACGCACACGTTTGCGtatcatat                                                  <  1:558455/149‑1 (MQ=255)
                                                                                             ttttAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGtatcatatc                                                 <  1:551239/149‑1 (MQ=255)
                                                                                               ttAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAg                                               <  1:45636/149‑1 (MQ=255)
                                                                                               ttAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAg                                               >  1:457406/1‑149 (MQ=255)
                                                                                                 aaaaTCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGaa                                             >  1:805190/1‑149 (MQ=255)
                                                                                                       ttttGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAgg                                       >  1:983718/1‑149 (MQ=255)
                                                                                                        tttGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAggg                                      >  2:409612/1‑149 (MQ=255)
                                                                                                             tGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGga                                 >  2:19766/1‑149 (MQ=255)
                                                                                                             ggTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCTAACTACGTCCGACAAACAAAGTTCGTGCGAATTACCGCTAAAATCGAGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCACAAAAAAGGGCCGga                                 <  2:83768/148‑1 (MQ=255)
                                                                                                               tCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGgatg                               <  1:592221/149‑1 (MQ=255)
                                                                                                                 aggTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACTAACAAATTTCGTGCGAATTACCGCTAAAATCGCGATTTCCTGCAGACGCACACGTTTGCGTCTCATATCAGAAAAAAGGGCCGgatgat                             >  1:193417/2‑149 (MQ=255)
                                                                                                                      tgCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGATCTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGGATGATTCCAg                        <  2:868741/149‑1 (MQ=255)
                                                                                                                        ccGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGGATGATTCCAGcc                      >  1:167714/1‑149 (MQ=255)
                                                                                                                        ccGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGGATGATTCCAGcc                      >  2:776286/1‑149 (MQ=255)
                                                                                                                        ccGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTGTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGGATGATTCCAGcc                      <  1:479792/149‑1 (MQ=255)
                                                                                                                        ccGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCCTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCATAAAAAAGGGCCGGATGATTCCAGcc                      >  2:819901/1‑149 (MQ=255)
                                                                                                                         cGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCACAAAAAAGGGCCGGATGATTCCAGccc                     >  2:745702/1‑149 (MQ=255)
                                                                                                                          ggATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGGATGATTCCAGCCCt                    >  1:385980/1‑149 (MQ=255)
                                                                                                                               ttttGTCGGATGCGGCGTGCACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTAGCATATCAGAAAAAAGGGCCGGATGATTCAAGCCCGGTAtt               >  1:926568/1‑149 (MQ=255)
                                                                                                                                tttGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGGTTTCCGGCAGACGCACACGTTTGCGTATCTTATCAGAAAAAAGGGCCGGATGTTTCCTGCCCTGTAttt              >  1:282159/1‑149 (MQ=39)
                                                                                                                                  tGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGGATGATTCCAGCCCTGTAttttt            <  2:966528/149‑1 (MQ=255)
                                                                                                                                  tGTCGGAAGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTGGTGCGATTTACCGCTAAAATCGCGTTTTCCTGCCGACGCACACGTTTTCGTATCATTTCAGAAAAAAGGGCCGGATGATTCCAGCCCTGTAttttt            <  2:90364/149‑1 (MQ=255)
                                                                                                                                     cGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGGATGATTCCAGCCCTGTATTTTTACt         <  2:127395/149‑1 (MQ=255)
                                                                                                                                       gataaggcgTTCACGCCGCATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGGATGATTCCAGCCCTGTATTTTTACTTg       >  1:224229/6‑149 (MQ=255)
                                                                                                                                       gataaggcgTTCACGCCGCATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGGATGATTCCAGCCCTGTATTTTTACTTg       >  2:32434/6‑149 (MQ=255)
                                                                                                                                         tggggcgTTAACGCCTTATCCGAAGTACGTACGACAAACATATTTCGTGCGAATTACCGGTAATATCTGGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGACGGATTATTCCAGCCCTGTATTTTTACTTGCt     <  2:932735/146‑1 (MQ=255)
                                                                                                                                            ggcgTTCACGCCGCATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGGATGATTCCAGCCCTGTATTTTTACTTGCTaaa  <  1:32434/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTTTTGCTGTCTGGTGTGCCGGATGTTTTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGTCCGACAAACAAATTTCGTGCGAATTACCGCTAAAATCGCGTTTTCCTGCAGACGCACACGTTTGCGTATCATATCAGAAAAAAGGGCCGGATGATTCCAGCCCTGTATTTTTACTTGCTAAA  >  NC_000913/2597554‑2597842

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