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.6% 803.9 / ‑4.9 246M158V (ATG→GTG) yidEputative transport protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (135/110);  total (135/111)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.45e-01

Only 100 of 253 total aligned reads displayed.

TCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCAC  >  NC_000913/3863339‑3863518
                                                                                          |                                                                                         
 ctcGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCAcccc                                                                                        <  2:1258334/93‑1 (MQ=255)
  tcGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACcccc                                                                                       >  2:5866502/1‑93 (MQ=255)
    gACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCAt                                                                                     >  2:2560300/1‑93 (MQ=255)
    gACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCAt                                                                                     >  1:143769/1‑93 (MQ=255)
     aCATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATc                                                                                    >  2:5252677/1‑93 (MQ=255)
     aCATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATc                                                                                    >  2:6253267/1‑93 (MQ=255)
     aCATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATc                                                                                    <  1:5455110/93‑1 (MQ=255)
       aTTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTg                                                                                  >  1:3392458/1‑93 (MQ=255)
       aTTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGGAACTCACCCCCATCTg                                                                                  >  1:1798741/1‑93 (MQ=255)
        ttGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGa                                                                                 >  1:3189986/1‑93 (MQ=255)
          gACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATc                                                                               >  1:5030259/1‑93 (MQ=255)
          gACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATc                                                                               <  1:5552386/93‑1 (MQ=255)
           aCGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTCACTCACCCCCATCTGATCg                                                                              >  1:4536461/1‑93 (MQ=255)
            cgcgGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGa                                                                             >  2:3318404/1‑93 (MQ=255)
             gcgGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGAc                                                                            <  1:1053463/93‑1 (MQ=255)
               ggAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCa                                                                          <  2:822237/93‑1 (MQ=255)
               ggAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCa                                                                          <  1:2609056/93‑1 (MQ=255)
                gAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCAt                                                                         >  1:3309287/1‑93 (MQ=255)
                 aaaaTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCAtt                                                                        <  1:6133164/93‑1 (MQ=255)
                    aTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTcc                                                                     >  1:2381860/1‑93 (MQ=255)
                    aTAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTcc                                                                     >  1:1685681/1‑93 (MQ=255)
                     tAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCa                                                                    >  2:6158091/1‑93 (MQ=255)
                      aaCCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCAt                                                                   <  2:4540087/93‑1 (MQ=255)
                       aCCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCAtt                                                                  >  2:1289734/1‑93 (MQ=255)
                       aCCCGCAACATCCACATGGAGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCAtt                                                                  >  1:233054/1‑93 (MQ=255)
                        cccGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTg                                                                 <  2:6709972/93‑1 (MQ=255)
                         ccGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTgg                                                                <  2:2541951/93‑1 (MQ=255)
                          cGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCGCCCCCATCTGATCGACCATTTCCATTGgt                                                               <  1:4267545/93‑1 (MQ=255)
                           gCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGgtg                                                              <  2:3114544/93‑1 (MQ=255)
                             aaCATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTa                                                            >  1:6375870/1‑93 (MQ=255)
                             aaCATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTa                                                            >  1:4765454/1‑93 (MQ=255)
                              aCATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAc                                                           >  1:3021810/1‑93 (MQ=255)
                              aCATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAc                                                           >  1:5934153/1‑93 (MQ=255)
                               cATCCACATGGTGAACAAAATCCCGGAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAcc                                                          >  1:864376/1‑93 (MQ=255)
                               cATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAcc                                                          >  1:707787/1‑93 (MQ=255)
                                aTCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAccc                                                         <  2:1555331/93‑1 (MQ=255)
                                aTCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTAccc                                                         <  2:2433767/93‑1 (MQ=255)
                                 tCCACATGGTGAACAAAATCCCGCAAACGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCa                                                        <  2:6431103/93‑1 (MQ=255)
                                   cacaTGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAgg                                                      >  1:5122984/1‑93 (MQ=255)
                                   cacaTGCCGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTCACTCACCCCCATCTGATCGCCCATTTCCATTGGTGTACCCAgg                                                      <  1:1258071/93‑1 (MQ=255)
                                     caTGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTc                                                    >  2:3685723/1‑93 (MQ=255)
                                      aTGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTcg                                                   >  1:3871676/1‑93 (MQ=255)
                                       tGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTcgc                                                  >  1:6336626/1‑93 (MQ=255)
                                        ggTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTcgcg                                                 <  2:5423517/93‑1 (MQ=255)
                                        ggTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTcgcg                                                 >  1:1139901/1‑93 (MQ=255)
                                        ggTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGCTcgcg                                                 >  1:4690037/1‑93 (MQ=255)
                                         gTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCgcgc                                                >  2:5472102/1‑93 (MQ=255)
                                              cAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAAt                                           <  1:887113/93‑1 (MQ=255)
                                               aaaaTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATc                                          >  2:3237337/1‑93 (MQ=255)
                                               aaaaTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATc                                          <  1:3946715/93‑1 (MQ=255)
                                               aaaaTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATc                                          <  1:1907776/93‑1 (MQ=255)
                                               aaaaTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATc                                          >  2:1509062/1‑93 (MQ=255)
                                               aaaaTCCCGCAAATGCCGAATGGATACACCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAAta                                          >  2:2291413/1‑92 (MQ=255)
                                                aaaTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATct                                         >  1:6281643/1‑93 (MQ=255)
                                                aaaTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATct                                         >  1:4960537/1‑93 (MQ=255)
                                                 aaTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCACAATctg                                        >  2:1134782/1‑93 (MQ=255)
                                                   tCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgct                                      >  2:1393613/1‑93 (MQ=255)
                                                    cccGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATctgctg                                     >  1:2155585/1‑93 (MQ=255)
                                                      cGCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgcc                                   >  2:4966731/1‑93 (MQ=255)
                                                       gCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccc                                  <  2:2838355/93‑1 (MQ=255)
                                                       gCAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccc                                  <  2:4267659/93‑1 (MQ=255)
                                                        cAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccct                                 >  2:3147450/1‑93 (MQ=255)
                                                        cAAATGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccct                                 >  1:5845760/1‑93 (MQ=255)
                                                          aaTGCCGAATTGATACGCCATCGCGTAACTCCCCCCCATCTGATCGACCATTCCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccctgc                               <  2:4432479/93‑1 (MQ=255)
                                                           aTGCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctgccctgcc                              >  2:5113441/1‑93 (MQ=255)
                                                             gCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGcccc                            <  2:1023832/93‑1 (MQ=255)
                                                             gCCGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGcccc                            >  1:1616455/1‑93 (MQ=255)
                                                              ccGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGccccc                           <  2:4648823/93‑1 (MQ=255)
                                                               cGAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCa                          <  2:3362991/93‑1 (MQ=255)
                                                                gAATGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAg                         >  1:373589/1‑93 (MQ=255)
                                                                 aaTGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAgc                        >  1:6591862/1‑93 (MQ=255)
                                                                 aaTGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATATGCTGCCCTGCCCCCAgc                        >  1:2902410/1‑93 (MQ=255)
                                                                  aTGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAgcg                       >  2:2031580/1‑93 (MQ=255)
                                                                   tGGATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAgcgc                      <  1:1382779/93‑1 (MQ=255)
                                                                     gATACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTg                    >  2:582570/1‑93 (MQ=255)
                                                                      aTACGCCATCGCGTAACTCTCCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTgg                   <  1:801035/93‑1 (MQ=255)
                                                                       tACGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGc                  <  2:5027879/93‑1 (MQ=255)
                                                                        aCGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCg                 <  1:4115165/93‑1 (MQ=255)
                                                                         cGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGt                >  2:4216520/1‑93 (MQ=255)
                                                                         cGCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAACCTGCTGCCCTGCCCCCAGCGCTGGCGt                <  2:537391/93‑1 (MQ=255)
                                                                          gCCATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTa               <  1:4140757/93‑1 (MQ=255)
                                                                           ccATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTAt              >  1:4438859/1‑93 (MQ=255)
                                                                           ccATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTAt              >  2:3980755/1‑93 (MQ=255)
                                                                            cATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTAtt             <  1:3161819/93‑1 (MQ=255)
                                                                            cATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTAtt             >  2:1683936/1‑93 (MQ=255)
                                                                            cATCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTAtt             >  2:2057217/1‑93 (MQ=255)
                                                                             aTCGCGTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTg            >  2:2738854/1‑93 (MQ=255)
                                                                               cgcgTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGt          <  1:2243616/93‑1 (MQ=255)
                                                                               cgcgTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGt          <  1:1068385/93‑1 (MQ=255)
                                                                                gcgTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTa         >  2:2226644/1‑93 (MQ=255)
                                                                                 cgTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTaa        >  2:2431227/1‑93 (MQ=255)
                                                                                  gTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCCGCGCTGGCGTATTGGTAAc       <  1:6525854/93‑1 (MQ=255)
                                                                                  gTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAAc       <  2:5193468/93‑1 (MQ=255)
                                                                                  gTAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAAc       <  1:3729776/93‑1 (MQ=255)
                                                                                   tAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAAcc      <  1:4201670/93‑1 (MQ=255)
                                                                                   tAACTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAAcc      <  2:178514/93‑1 (MQ=255)
                                                                                    aaCTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCg     >  2:4855439/1‑93 (MQ=255)
                                                                                    aaCTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCg     >  2:6228241/1‑93 (MQ=255)
                                                                                     aCTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGc    >  2:6124445/1‑93 (MQ=255)
                                                                                     aCTCACCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGc    >  2:3552261/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCAC  >  NC_000913/3863339‑3863518

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