Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 629,933 T→C W273R (TGG→CGG)  cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913629,9330TC99.5% 606.7 / ‑4.8 198W273R (TGG→CGG) cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (75/122);  total (75/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

Only 100 of 246 total aligned reads displayed.

CCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGC  >  NC_000913/629843‑630022
                                                                                          |                                                                                         
 cGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGt                                                                                        <  2:9780104/93‑1 (MQ=255)
 cGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGt                                                                                        <  2:3335605/93‑1 (MQ=255)
   gccTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTa                                                                                      <  1:2709042/92‑1 (MQ=255)
   aCCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTa                                                                                      <  2:9951172/93‑1 (MQ=255)
    ccTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTAc                                                                                     <  1:1946216/93‑1 (MQ=255)
    ccTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTAc                                                                                     <  2:843420/93‑1 (MQ=255)
    ccTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGCGGTTAc                                                                                     >  1:4432022/1‑93 (MQ=255)
     cTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTAct                                                                                    >  1:1820981/1‑93 (MQ=255)
     cTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTAct                                                                                    <  2:8524463/93‑1 (MQ=255)
     cTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTAct                                                                                    >  2:2845144/1‑93 (MQ=255)
       gggCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTActgc                                                                                  <  1:9470976/93‑1 (MQ=255)
       gggCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTActgc                                                                                  >  2:4642268/1‑93 (MQ=255)
       gggCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTActgc                                                                                  >  2:7034111/1‑93 (MQ=255)
       gggCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTActgc                                                                                  >  2:795164/1‑93 (MQ=255)
        ggCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTActgct                                                                                 <  1:3352488/93‑1 (MQ=255)
        ggCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGCGGTTActgct                                                                                 >  1:9924510/1‑93 (MQ=255)
          cACCGTACTTTGACTTTACCGGCATGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTgg                                                                               <  2:5527326/93‑1 (MQ=255)
           aCCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGAGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGc                                                                              >  1:9672872/1‑93 (MQ=255)
             cGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGccc                                                                            <  2:3988621/93‑1 (MQ=255)
               tACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCcg                                                                          >  2:1153008/1‑93 (MQ=255)
               tACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGTGCGGTTACTGCCGGGCCcg                                                                          >  2:6241887/1‑93 (MQ=255)
               tACTTTGACTTTACCGGCGTGCAGTTGACCTGGATACTAGTGGGTTACGATTTTGAGGCGGCGGTGCTGCCAGTGCGGTTACTGCTGGcacca                                                                          <  2:8942711/93‑2 (MQ=255)
                tcTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCcgc                                                                         <  2:7215483/92‑1 (MQ=255)
                aCTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCcgc                                                                         <  1:5694709/93‑1 (MQ=255)
                    tGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGa                                                                     <  1:4557797/93‑1 (MQ=255)
                     gACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGAc                                                                    <  2:2686648/93‑1 (MQ=255)
                     gACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGAc                                                                    <  2:8907652/93‑1 (MQ=255)
                      aCTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGTGCGGTTACTGGTGGTCCCGCGTGACt                                                                   >  1:392985/1‑93 (MQ=255)
                      aCTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGTGCGGGTACTGCTGGCCCCGCGTGACt                                                                   >  1:136937/1‑93 (MQ=255)
                       cTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTa                                                                  <  1:2967589/93‑1 (MQ=255)
                         ttACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTAcc                                                                <  2:2271211/93‑1 (MQ=255)
                           aCCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGGGACTTCctc                                                              >  2:5943046/1‑93 (MQ=255)
                           aCCGACGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACctc                                                              >  2:6016747/1‑93 (MQ=255)
                             cGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACctctc                                                            <  1:7322922/93‑1 (MQ=255)
                              ggCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACctctct                                                           <  2:8733449/93‑1 (MQ=255)
                               gCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACctctctc                                                          >  2:3026880/1‑92 (MQ=255)
                                cGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTAc                                                         <  2:3698047/93‑1 (MQ=255)
                                 gTGCAGTTGGCCTGGATGCTGGTGGGTTCCGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTAcc                                                        <  2:3099408/93‑1 (MQ=255)
                                 gTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTAcc                                                        >  1:804464/1‑93 (MQ=255)
                                 gTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTAcc                                                        >  2:8166488/1‑93 (MQ=255)
                                  tGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGACCGCCGTGAGTACCTCTCTACCt                                                       >  1:3508347/1‑93 (MQ=255)
                                  tGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTGACTGCTGGCCCCGCATCACTACCTGTCTatgt                                                       >  2:5319543/1‑90 (MQ=255)
                                  tGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTACTGCTGGGCCCGCGTGACTACCTCTCTACCt                                                       >  2:252615/1‑93 (MQ=255)
                                  tGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCt                                                       >  2:4178124/1‑93 (MQ=255)
                                     aGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTcc                                                    <  2:5835703/93‑1 (MQ=255)
                                      gTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGGCCCGCGCGACTACCTCTCTCCCTTCCt                                                   >  1:912351/1‑93 (MQ=255)
                                      gTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCt                                                   <  2:5138196/93‑1 (MQ=255)
                                       ttGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTg                                                  >  1:5366803/1‑93 (MQ=255)
                                       ttGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTg                                                  >  2:7855026/1‑93 (MQ=255)
                                          aCCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaaa                                               >  2:7107206/1‑93 (MQ=255)
                                          aCCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGGTACTGCTGGCCCCGCGGGACTACCGCTCTAGCTTCCTGaaa                                               >  2:6318888/1‑93 (MQ=255)
                                           ccTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaaaa                                              >  1:5433073/1‑93 (MQ=255)
                                             tGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                            <  1:9622336/93‑1 (MQ=255)
                                             tGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                            >  1:646333/1‑93 (MQ=255)
                                             cggATGCTCGTGGGTTACGGTGTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                            <  2:6586896/92‑1 (MQ=255)
                                              ggATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCg                                           >  1:126888/1‑93 (MQ=255)
                                                aTGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCggg                                         <  1:6100732/93‑1 (MQ=255)
                                                 tGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTATACCTTCCTGAAAATCGGGa                                        >  2:3588731/1‑93 (MQ=255)
                                                  gCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGAc                                       <  1:4154893/93‑1 (MQ=255)
                                                  gCTGGTAGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGAc                                       <  1:5927861/93‑1 (MQ=255)
                                                    tggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTa                                     <  2:8049504/93‑1 (MQ=255)
                                                     ggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAt                                    >  1:7720751/1‑93 (MQ=255)
                                                      gtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtc                                   >  1:4540992/1‑93 (MQ=255)
                                                      gtggGTTACGGTTTTGTGGCGGCGGAGCTGCCGGTGCGGTTACTGCCGGCCCCGGGTGGCTACCTCTCTACCTTCCTGAAAATCGGGACTAtc                                   >  2:8397178/1‑93 (MQ=255)
                                                       tggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcg                                  <  1:8951627/93‑1 (MQ=255)
                                                       tggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcg                                  >  1:5071074/1‑93 (MQ=255)
                                                       tggGTTACGGTTTTGTGGCGGCGGTGCTGCCCGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcg                                  >  2:9452764/1‑93 (MQ=255)
                                                        tggTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcgt                                 <  1:1875532/92‑1 (MQ=255)
                                                        gggTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcgt                                 <  1:3103003/93‑1 (MQ=255)
                                                        gggTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcgt                                 <  2:9520294/93‑1 (MQ=255)
                                                        gggTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcgt                                 <  1:6749384/93‑1 (MQ=255)
                                                         ggTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcgtc                                <  1:8162497/93‑1 (MQ=255)
                                                         ggTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAtcgtc                                <  2:3455741/93‑1 (MQ=255)
                                                           ttACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                              <  2:5680627/93‑1 (MQ=255)
                                                            tACGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGt                             >  1:6230424/1‑93 (MQ=255)
                                                              cGGTTTTGTGGCGTCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCt                           <  2:2514596/93‑1 (MQ=255)
                                                              cGGTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCt                           <  2:6515511/93‑1 (MQ=255)
                                                               ggTTTTGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGCCGGGCTg                          >  2:64194/1‑93 (MQ=255)
                                                                 ttttGTGGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGc                        >  2:6891884/1‑93 (MQ=255)
                                                                  tttGTGGCGGCGGTGCTGCCGGGGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATAGGGACTATCGTCGGTCTGGCg                       >  1:2377793/1‑93 (MQ=255)
                                                                     gtgGCGGCGGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGtg                    >  1:2786100/1‑92 (MQ=255)
                                                                        gcggcgGTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGc                 <  1:6892822/93‑1 (MQ=255)
                                                                         gggcggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCa                <  2:4480945/92‑1 (MQ=255)
                                                                         cggcggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCa                <  2:1446579/93‑1 (MQ=255)
                                                                           gcggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAtt              <  2:6176377/93‑1 (MQ=255)
                                                                            cggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAttt             <  1:9405044/93‑1 (MQ=255)
                                                                             ggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAtttt            <  1:658628/93‑1 (MQ=255)
                                                                             ggTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAtttt            >  1:873116/1‑93 (MQ=255)
                                                                              gTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTg           >  1:7406151/1‑93 (MQ=255)
                                                                              gTGCTGCCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGAATTTTg           >  2:1577697/1‑93 (MQ=255)
                                                                                  tgtcGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTa       <  2:4975870/90‑1 (MQ=255)
                                                                                  tgcCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTa       <  2:5396066/93‑1 (MQ=255)
                                                                                  tgcCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGGAGGCATTTTGATTa       >  1:3818487/1‑93 (MQ=255)
                                                                                   gcCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTAt      >  1:1942432/1‑93 (MQ=255)
                                                                                   gcCGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTAt      <  2:7242375/93‑1 (MQ=255)
                                                                                    ccGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATg     <  1:213260/93‑1 (MQ=255)
                                                                                    ccGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATg     <  1:3210044/93‑1 (MQ=255)
                                                                                     cGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATgc    <  2:864323/93‑1 (MQ=255)
                                                                                     cGGTGCGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATgc    <  1:102661/93‑1 (MQ=255)
                                                                                       gtgcGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATgcgc  >  1:4325362/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGC  >  NC_000913/629843‑630022

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