Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 577,783 T→C C152R (TGC→CGC)  ybcT → endopeptidase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913577,7830TC99.6% 1450.7 / ‑3.9 454C152R (TGC→CGC) ybcTendopeptidase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (270/182);  total (272/182)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

Only 100 of 470 total aligned reads displayed.

ACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAG  >  NC_000913/577693‑577872
                                                                                          |                                                                                         
acCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAgcgc                                                                                         >  1:8215041/1‑93 (MQ=255)
 ccGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCa                                                                                        >  2:2301227/1‑93 (MQ=255)
 ccGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCa                                                                                        <  2:7886085/93‑1 (MQ=255)
   gCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGa                                                                                      >  1:2736682/1‑93 (MQ=255)
   gCTGAACGGGATTATTTCACACTCAGCGACAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGAATATTAATGAGCAGCGCAGa                                                                                      >  1:3163282/1‑93 (MQ=255)
     tGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATa                                                                                    <  2:5373545/93‑1 (MQ=255)
     tGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATa                                                                                    <  1:8633433/93‑1 (MQ=255)
     tGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATa                                                                                    >  2:2088822/1‑93 (MQ=255)
      gAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATag                                                                                   >  1:10004739/1‑93 (MQ=255)
         cGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGc                                                                                <  2:8163344/93‑1 (MQ=255)
         cGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGc                                                                                <  2:5768370/93‑1 (MQ=255)
           ggATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTg                                                                              >  1:709922/1‑93 (MQ=255)
            gATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGa                                                                             <  1:9373420/93‑1 (MQ=255)
              ttattTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGAcc                                                                           >  1:3342979/1‑93 (MQ=255)
               tattTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCa                                                                          <  2:2662158/93‑1 (MQ=255)
                    taCCCGCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGAAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCg                                                                     <  2:8000523/92‑1 (MQ=255)
                    cACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCg                                                                     >  1:6467749/1‑93 (MQ=255)
                     aCCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGa                                                                    >  1:4921654/1‑93 (MQ=255)
                      cccTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGAt                                                                   >  1:5127353/1‑93 (MQ=255)
                      cccTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGAt                                                                   <  2:3976321/93‑1 (MQ=255)
                        cTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAGCCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATgg                                                                 <  1:3672480/93‑1 (MQ=255)
                        cTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATgg                                                                 <  1:4819307/93‑1 (MQ=255)
                         tCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATggg                                                                <  1:135773/93‑1 (MQ=255)
                            gagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCaa                                                             <  2:6457729/93‑1 (MQ=255)
                              gagagGCTGATTACTATGCAAAAACAACTGGGAGGAAACCCGAAGTATATTAATTAGCAGCGCAGATAGAGCTGACCAAATCAGTGGGCAAcc                                                           >  2:1559002/1‑92 (MQ=255)
                              gagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACt                                                           >  1:459069/1‑93 (MQ=255)
                                  ggCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATg                                                       <  1:6231386/93‑1 (MQ=255)
                                   gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGc                                                      >  2:5054261/1‑93 (MQ=255)
                                     tGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCaa                                                    <  1:5398337/93‑1 (MQ=255)
                                     tGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCaa                                                    <  1:6024140/93‑1 (MQ=255)
                                      gATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAat                                                   >  1:8418042/1‑93 (MQ=255)
                                      gATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAat                                                   <  1:1751262/93‑1 (MQ=255)
                                      gATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAat                                                   >  2:6955923/1‑93 (MQ=255)
                                       aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACGATATCGATGGGCAACTCATGCAatt                                                  >  2:4444197/1‑93 (MQ=255)
                                       aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAatt                                                  >  1:2172273/1‑93 (MQ=255)
                                       aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAatt                                                  >  1:8957777/1‑93 (MQ=255)
                                       aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAatt                                                  <  1:8030326/93‑1 (MQ=255)
                                       aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAatt                                                  >  1:444562/1‑93 (MQ=255)
                                        tCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAatta                                                 >  2:3029266/1‑93 (MQ=255)
                                        tCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAatta                                                 >  2:1074092/1‑93 (MQ=255)
                                        tCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAatta                                                 >  1:549234/1‑93 (MQ=255)
                                          aCTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAattatt                                               >  2:7329570/1‑93 (MQ=255)
                                          aCTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAattatt                                               >  2:7539910/1‑93 (MQ=255)
                                           cTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTAttt                                              >  1:7304013/1‑93 (MQ=255)
                                           cTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTAttt                                              >  1:4246980/1‑93 (MQ=255)
                                                 aaaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCa                                        >  1:1167267/1‑93 (MQ=255)
                                                  aaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCaa                                       >  1:2384749/1‑93 (MQ=255)
                                                  aaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCaa                                       >  2:9570965/1‑93 (MQ=255)
                                                  aaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCaa                                       >  1:9514310/1‑93 (MQ=255)
                                                  aaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCaa                                       <  1:6316363/93‑1 (MQ=255)
                                                   aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAAt                                      >  1:1954093/1‑93 (MQ=255)
                                                   aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAAt                                      <  1:6881872/93‑1 (MQ=255)
                                                   aaaCAACTGGAAGGAACACAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAACTATTTTGAGCAAt                                      <  2:4230524/93‑1 (MQ=255)
                                                    aacaacTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCCATTATTTTGAGCAATa                                     >  1:5893984/1‑93 (MQ=255)
                                                       aacTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacac                                  >  2:1879275/1‑93 (MQ=255)
                                                        acTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacaca                                 <  1:5587385/93‑1 (MQ=255)
                                                        acTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacaca                                 <  1:6126257/93‑1 (MQ=255)
                                                        acTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacaca                                 <  1:8288263/93‑1 (MQ=255)
                                                        acTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacaca                                 <  1:9967471/93‑1 (MQ=255)
                                                        acTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacaca                                 <  2:6062266/93‑1 (MQ=255)
                                                         cTGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacacac                                <  2:7730122/93‑1 (MQ=255)
                                                          tGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcg                               >  2:4945669/1‑93 (MQ=255)
                                                          tGGAAGGAACCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcg                               >  2:6536191/1‑93 (MQ=255)
                                                           ggaaggaaCCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc                              <  1:7472584/93‑1 (MQ=255)
                                                           ggaaggaaCCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc                              >  1:9864838/1‑93 (MQ=255)
                                                           ggaaggaaCCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc                              >  2:2527273/1‑93 (MQ=255)
                                                           ggaaggaaCCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc                              >  1:9584425/1‑93 (MQ=255)
                                                            gaaggaaCCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgcg                             >  2:8436369/1‑93 (MQ=255)
                                                              aggaaCCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCt                           >  2:554878/1‑93 (MQ=255)
                                                              aggaaCCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCt                           >  1:5752742/1‑93 (MQ=255)
                                                              aggaaCCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCt                           <  1:1541803/93‑1 (MQ=255)
                                                               ggaaCCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCtt                          <  1:1923042/93‑1 (MQ=255)
                                                                 aaCCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTcc                        <  2:7200482/93‑1 (MQ=255)
                                                                 aaCCCAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCCTcc                        <  1:1674041/93‑1 (MQ=255)
                                                                   cccAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                      >  2:9803911/1‑93 (MQ=255)
                                                                    ccAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGc                     >  1:522329/1‑93 (MQ=255)
                                                                     cAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                    >  2:853912/1‑93 (MQ=255)
                                                                     cAGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                    >  2:9336504/1‑93 (MQ=255)
                                                                      aGAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCgg                   >  2:471102/1‑93 (MQ=255)
                                                                       gAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGa                  >  2:3314625/1‑93 (MQ=255)
                                                                       gAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGa                  >  2:625399/1‑93 (MQ=255)
                                                                       gAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGa                  <  1:3540896/93‑1 (MQ=255)
                                                                       gAAGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGa                  >  1:2938913/1‑93 (MQ=255)
                                                                        aaGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAg                 <  1:809205/93‑1 (MQ=255)
                                                                         aGTATATTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGt                >  2:593776/1‑93 (MQ=255)
                                                                           tataTTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGtat              <  1:4066284/93‑1 (MQ=255)
                                                                            atatTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCCTCCAGCGGAGtata             >  2:3710094/1‑93 (MQ=255)
                                                                              atTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATaaa           <  2:5144812/93‑1 (MQ=255)
                                                                              atTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATaaa           <  2:4188719/93‑1 (MQ=255)
                                                                              atTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATaaa           <  2:128874/93‑1 (MQ=255)
                                                                              atTAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATaaa           <  1:7748457/93‑1 (MQ=255)
                                                                              atTAATGAGCACCGCAGATAGAGCTGACCATATCGATGGGCAACTCAGGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATaaa           <  1:9766984/93‑1 (MQ=255)
                                                                              atTAAGGAGGAGCGCAGAGAGAGCTGACCATATCGATGGGCAACTCATGCAATTAGTTTGAGCAGTACACACGCTCTTCCAGCGGAGTATaaa           <  2:2792067/93‑1 (MQ=255)
                                                                               ttAATGAGGAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGAGCTTCCAGCGGAGTATAAAt          <  1:7608646/93‑1 (MQ=255)
                                                                               ttAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAAt          >  2:8777595/1‑93 (MQ=255)
                                                                               ttAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAAt          >  2:9161336/1‑93 (MQ=255)
                                                                                tAATGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGGGCAATACCCACGCGCTTCAAGCGGATAATAAAtc         >  1:7830923/1‑92 (MQ=255)
                                                                                 aaTGAGCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGc        >  1:3921590/1‑93 (MQ=255)
                                                                                      gCAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaaa   <  2:123446/93‑1 (MQ=255)
                                                                                       cAGCGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAg  <  1:9016753/93‑1 (MQ=255)
                                                                                          |                                                                                         
ACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAG  >  NC_000913/577693‑577872

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: