Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,863,429 T→C M158V (ATG→GTG)  yidE ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,863,4290TC99.8% 1379.9 / ‑5.0 424M158V (ATG→GTG) yidEputative transport protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (236/187);  total (237/187)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

Only 100 of 445 total aligned reads displayed.

TCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCAC  >  NC_000913/3863339‑3863518
                                                                                          |                                                                                          
 ctcGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCAcccc                                                                                         >  2:8219085/1‑93 (MQ=255)
 ctcGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCAcccc                                                                                         <  1:2000321/93‑1 (MQ=255)
   cGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCa                                                                                       <  1:8201028/93‑1 (MQ=255)
     aCATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATc                                                                                     >  1:4370391/1‑93 (MQ=255)
     aCATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATc                                                                                     >  2:6557837/1‑93 (MQ=255)
     aCATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATc                                                                                     >  2:3757298/1‑93 (MQ=255)
      cATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCt                                                                                    <  2:9623546/93‑1 (MQ=255)
      cATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCt                                                                                    >  2:8747741/1‑93 (MQ=255)
       aTTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTg                                                                                   >  2:2370221/1‑93 (MQ=255)
        ttGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGa                                                                                  >  2:5853537/1‑93 (MQ=255)
        ttGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGa                                                                                  >  1:6757764/1‑93 (MQ=255)
          gACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATc                                                                                >  2:3626518/1‑93 (MQ=255)
          gACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATc                                                                                >  2:152238/1‑93 (MQ=255)
            cgcgGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCgc                                                                              >  2:2861195/1‑92 (MQ=255)
               ggAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCa                                                                           >  2:1097174/1‑93 (MQ=255)
                 aaaaTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCAtt                                                                         <  1:8737519/93‑1 (MQ=255)
                  aaaTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCAttt                                                                        <  1:1765974/93‑1 (MQ=255)
                  aaaTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCAttt                                                                        >  1:2692223/1‑93 (MQ=255)
                  aaaTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCAttt                                                                        <  1:7531089/93‑1 (MQ=255)
                   aaTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTc                                                                       >  1:5947691/1‑93 (MQ=255)
                   aaTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTc                                                                       <  1:8874327/93‑1 (MQ=255)
                   aaTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTc                                                                       >  1:3828703/1‑93 (MQ=255)
                    aTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTcc                                                                      >  1:2139656/1‑93 (MQ=255)
                    aTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTcc                                                                      >  2:5089941/1‑93 (MQ=255)
                     tAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATATGCCATCGCGTAACTTACCCCCATATGATCGACCATTTCCa                                                                     >  1:9012460/1‑93 (MQ=255)
                     tAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCa                                                                     <  2:755584/93‑1 (MQ=255)
                     tAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCa                                                                     >  1:9207208/1‑93 (MQ=255)
                      aaCCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCAt                                                                    >  2:1659922/1‑93 (MQ=255)
                        cccGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTg                                                                  >  1:9902903/1‑93 (MQ=255)
                         ccGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTgg                                                                 <  1:4307874/93‑1 (MQ=255)
                         ccGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTgg                                                                 >  2:5336119/1‑93 (MQ=255)
                          cGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGgt                                                                <  2:7815206/93‑1 (MQ=255)
                           gCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGgtg                                                               <  2:3022280/93‑1 (MQ=255)
                              aCATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAc                                                            >  1:5761378/1‑93 (MQ=255)
                               cATCCACATGGTGAACAACATCCCGCAAATGCAGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAcc                                                           <  1:1642742/93‑1 (MQ=255)
                               cATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGACCGACCATTTCCATTGGTGTAcc                                                           <  2:7768770/93‑1 (MQ=255)
                                aTCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAccc                                                          <  1:7727016/93‑1 (MQ=255)
                                aTCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAccc                                                          <  1:1471805/93‑1 (MQ=255)
                                  ccACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAg                                                        <  2:2589489/93‑1 (MQ=255)
                                  ccACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAg                                                        <  1:3527313/93‑1 (MQ=255)
                                  ccACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAg                                                        <  2:9381230/93‑1 (MQ=255)
                                   cacaTGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAgg                                                       >  1:6914627/1‑93 (MQ=255)
                                     caTGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTc                                                     >  1:22995/1‑93 (MQ=255)
                                     caTGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTc                                                     >  1:8642075/1‑93 (MQ=255)
                                       tGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTcgc                                                   <  1:6319459/93‑1 (MQ=255)
                                       tGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTcgc                                                   <  2:7020198/93‑1 (MQ=255)
                                       tGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTcgc                                                   >  1:8400771/1‑93 (MQ=255)
                                       tGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTcgc                                                   <  1:10002784/93‑1 (MQ=255)
                                        ggTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTcgcg                                                  <  2:6152567/93‑1 (MQ=255)
                                         gTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCgcgc                                                 <  2:4780994/93‑1 (MQ=255)
                                            aaCAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCaaa                                              >  2:4234603/1‑93 (MQ=255)
                                            aaCAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCaaa                                              <  1:5013358/93‑1 (MQ=255)
                                             aCAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCaaaa                                             <  1:5624688/93‑1 (MQ=255)
                                             aCAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCaaaa                                             >  1:5215396/1‑93 (MQ=255)
                                               aaaaTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATc                                           >  2:3494058/1‑93 (MQ=255)
                                               aaaaTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATc                                           <  2:10016879/93‑1 (MQ=255)
                                               aaaaTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATc                                           >  1:290310/1‑93 (MQ=255)
                                               aaaaTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATc                                           >  1:2378789/1‑93 (MQ=255)
                                                  aTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgc                                        <  2:4051778/93‑1 (MQ=255)
                                                  aTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgc                                        <  1:4064270/93‑1 (MQ=255)
                                                  aTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgc                                        >  1:7420648/1‑93 (MQ=255)
                                                   tCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgct                                       <  2:6681632/93‑1 (MQ=255)
                                                   tCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgct                                       >  1:4615476/1‑93 (MQ=255)
                                                    cccGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGAGCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgctg                                      <  2:9058489/93‑1 (MQ=255)
                                                     ccGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCtgctgc                                     <  1:7626670/93‑1 (MQ=255)
                                                       gCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccc                                   <  2:4904974/93‑1 (MQ=255)
                                                        cAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccct                                  <  2:205639/93‑1 (MQ=255)
                                                        cAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccct                                  >  1:3110877/1‑93 (MQ=255)
                                                        cAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccct                                  >  1:9765467/1‑93 (MQ=255)
                                                          aaTGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGTTCGCGCAAAATCTGctgccctgc                               >  2:5534247/1‑93 (MQ=255)
                                                          aaTGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccctgc                               <  1:1386717/93‑1 (MQ=255)
                                                           aTGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccctgcc                              >  2:8065158/1‑93 (MQ=255)
                                                            tGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCtgccctgccc                             >  1:8562119/1‑93 (MQ=255)
                                                            tGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCtgccctgccc                             <  1:1909292/93‑1 (MQ=255)
                                                             tccGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGcccc                            <  1:3274073/92‑1 (MQ=255)
                                                               cGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCACTGCCCCCa                          >  2:5805898/1‑93 (MQ=255)
                                                                 aaTGGATACGCCATCGCGTAACTCCCCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTGCCGCAAAATCTGCTCCCCTGCCCCCAgc                        >  1:5530017/1‑93 (MQ=255)
                                                                 aaTGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGccccccag                        <  1:805548/93‑4 (MQ=255)
                                                                 aaTGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAgc                        >  1:5721245/1‑93 (MQ=255)
                                                                  aTGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGccccccagc                       >  1:1171923/1‑89 (MQ=255)
                                                                  aTGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAgcg                       >  2:131122/1‑93 (MQ=255)
                                                                    ggATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCCAgcgc                      >  1:3392008/1‑93 (MQ=255)
                                                                    ggATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCt                     <  2:924853/93‑1 (MQ=255)
                                                                      aTACGCCATCGCGTAACTCTCCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTgg                   <  1:6143919/93‑1 (MQ=255)
                                                                       tACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGc                  >  1:7987435/1‑93 (MQ=255)
                                                                        aCGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCg                 >  2:2163719/1‑93 (MQ=255)
                                                                          gCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTa               <  2:6436306/93‑1 (MQ=255)
                                                                          gCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTa               >  1:9370493/1‑93 (MQ=255)
                                                                          gCCATCGCGTAACTCACCCCCATCGGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTa               >  2:7316712/1‑93 (MQ=255)
                                                                           ccATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTAt              >  1:7158459/1‑93 (MQ=255)
                                                                              tCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTgg           >  1:5410801/1‑93 (MQ=255)
                                                                               cgcgTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGt          <  1:1089626/93‑1 (MQ=255)
                                                                                gcgTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTa         >  1:487205/1‑93 (MQ=255)
                                                                                 cgTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTaa        <  1:7275835/93‑1 (MQ=255)
                                                                                    aaCTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCg     >  2:4402189/1‑93 (MQ=255)
                                                                                     aCTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTTGCGCAAAATCTGCTGCCCTGCCCCCCAGCGCTGGCGTATTGGTAACCg     >  2:274756/1‑93 (MQ=255)
                                                                                      cTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCCAGCGCTGGCGTATTGGTAACCGc    <  2:8452623/93‑1 (MQ=255)
                                                                                      cTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCa   <  2:6883376/93‑1 (MQ=255)
                                                                                      cTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCa   >  2:6176697/1‑93 (MQ=255)
                                                                                       tCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCAc  >  1:6489289/1‑93 (MQ=255)
                                                                                          |                                                                                          
TCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCAC  >  NC_000913/3863339‑3863518

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