Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,338,889 G→A 6.0% I120I (ATC→ATT)  mlaE ← intermembrane phospholipid transport system, integral membrane subunit MlaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,338,8890GA 6.0% 207.6 / 2.4 66I120I (ATC→ATT) mlaEintermembrane phospholipid transport system, integral membrane subunit MlaE
Reads supporting (aligned to +/- strand):  ref base G (25/37);  new base A (2/2);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

GCCCCAGATCCCCACGGCGACGAAAATAACCGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAGATAACCTTG  >  NC_000913/3338742‑3339032
                                                                                                                                                   |                                                                                                                                               
gCCCCAGATCCCCACGGCGACGAAAATAACCGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGAt                                                                                                                                               <  2:1249776/150‑1 (MQ=255)
       aTCCCCACGGCGACGAAAATAACCGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCg                                                                                                                                        <  2:472398/150‑1 (MQ=255)
        tCCCCACGGCGACGAAAATAACCGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACGCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCgg                                                                                                                                       <  2:758213/150‑1 (MQ=255)
        tCCCCACGGCGACGAAAATAACCGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCgg                                                                                                                                       <  2:113854/150‑1 (MQ=255)
         ccccACGGCGACGAAAATAACCGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCTTCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGt                                                                                                                                      <  1:910811/150‑1 (MQ=255)
              cGGCGACGAAAATAACCGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAgcg                                                                                                                                 <  2:366792/150‑1 (MQ=255)
                gCGACGAAAATAACCGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGcc                                                                                                                               <  2:534312/150‑1 (MQ=255)
                       aaaTAACCGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAg                                                                                                                        <  1:920132/150‑1 (MQ=255)
                       aaaTAACCGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAg                                                                                                                        <  2:1409697/150‑1 (MQ=255)
                         aTAACCGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCa                                                                                                                      <  2:275495/150‑1 (MQ=255)
                          tAACCGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCAc                                                                                                                     >  2:1087116/1‑150 (MQ=255)
                           aaCCGTCAAGAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTTGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACg                                                                                                                    <  2:736911/150‑1 (MQ=255)
                             ccGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGcc                                                                                                                  >  2:866472/1‑150 (MQ=255)
                                tCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCgg                                                                                                               >  1:223681/1‑150 (MQ=255)
                                 cAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGGAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGc                                                                                                              <  2:138467/150‑1 (MQ=255)
                                 cAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGc                                                                                                              >  1:1422474/1‑150 (MQ=255)
                                      gTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaaca                                                                                                         <  1:639916/150‑1 (MQ=255)
                                       tGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaacaa                                                                                                        <  1:1003372/150‑1 (MQ=255)
                                        ggTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaacaac                                                                                                       <  2:200922/150‑1 (MQ=255)
                                         gTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaacaaca                                                                                                      <  2:340957/150‑1 (MQ=255)
                                             tGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAAcgc                                                                                                  >  1:1179989/1‑150 (MQ=255)
                                                aaTAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGc                                                                                               <  1:241636/150‑1 (MQ=255)
                                                     ccccAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAacca                                                                                          >  2:947072/1‑150 (MQ=255)
                                                      cccagcccagAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAaccac                                                                                         <  1:724843/150‑1 (MQ=255)
                                                       ccagcccagAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAaccacc                                                                                        >  2:1080644/1‑150 (MQ=255)
                                                       ccagcccagAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCAATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAacc                                                                                           >  2:592691/1‑147 (MQ=255)
                                                       ccagcccagAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCAATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAacc                                                                                           <  1:592691/147‑1 (MQ=255)
                                                          gcccagAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGc                                                                                                                                         >  1:1335234/1‑98 (MQ=255)
                                                          gcccagAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGc                                                                                                                                         <  2:1335234/98‑1 (MQ=255)
                                                           cccagAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGcc                                                                                    <  1:1039913/150‑1 (MQ=255)
                                                                     gggAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGATCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACg                                                                          <  1:432255/150‑1 (MQ=255)
                                                                      ggAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGc                                                                         <  1:947072/150‑1 (MQ=255)
                                                                           aaTAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCAATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAg                                                                    >  1:465810/1‑150 (MQ=255)
                                                                           aaTAACACGACGCAGCGGATCCACCGACATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCCCATCAGGCCGATTACTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAg                                                                    >  2:365002/1‑150 (MQ=255)
                                                                            aTAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGc                                                                   >  1:689728/1‑150 (MQ=255)
                                                                             tAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCg                                                                  <  1:1080644/150‑1 (MQ=255)
                                                                                cccGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCCGCAAACAACACCGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATa                                                               >  2:1139635/1‑150 (MQ=255)
                                                                                cccGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCCGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATa                                                               >  1:708072/1‑150 (MQ=255)
                                                                                   gACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACg                                                            <  1:866472/150‑1 (MQ=255)
                                                                                     cGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGcc                                                          <  2:1314906/150‑1 (MQ=255)
                                                                                     cGCAGCGGATCCACCGCCAGCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTGCACGCAGTAGCGATAACGcc                                                          <  2:432911/150‑1 (MQ=255)
                                                                                        aGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGccacc                                                       >  1:1133996/1‑150 (MQ=255)
                                                                                         gCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGccacca                                                      <  1:923289/150‑1 (MQ=255)
                                                                                              tCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATa                                                 <  2:8210/150‑1 (MQ=255)
                                                                                                cACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaac                                                                                                          >  1:313147/1‑91 (MQ=255)
                                                                                                cACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAaac                                                                                                          <  2:313147/91‑1 (MQ=255)
                                                                                                     ccATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGAc                                          >  2:662416/1‑150 (MQ=255)
                                                                                                          atcTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGccacc                                                       >  2:1051888/1‑132 (MQ=255)
                                                                                                          atcTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGccacc                                                       <  1:1051888/132‑1 (MQ=255)
                                                                                                           tctcCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGttt                                    >  1:673136/1‑150 (MQ=255)
                                                                                                             tcCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTcc                                  >  1:1189160/1‑150 (MQ=255)
                                                                                                             tcCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTcc                                  >  2:1080992/1‑150 (MQ=255)
                                                                                                              ccATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCg                                 >  1:455943/1‑150 (MQ=255)
                                                                                                              ccATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCAATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCg                                 <  2:465810/150‑1 (MQ=255)
                                                                                                               cATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGc                                <  2:455943/150‑1 (MQ=255)
                                                                                                                 tACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCAc                              <  2:1422474/150‑1 (MQ=255)
                                                                                                                   cTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACta                            >  1:576391/1‑150 (MQ=255)
                                                                                                                     ggAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACtata                          <  1:783215/150‑1 (MQ=255)
                                                                                                                                 tGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAAc              >  1:471/1‑150 (MQ=255)
                                                                                                                                  gTAGCGAGAAACAGGCCGATTTATACAGTTTGCGACGAACCAGAACGCCCTGCAAACCACAACAATACAAAAAAAGGCCCCAGTTCACGCAGTAGCGATACCGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAAcc             >  2:1267906/1‑150 (MQ=255)
                                                                                                                                   tAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCa            <  2:673136/150‑1 (MQ=255)
                                                                                                                                    aGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAg           <  2:1133996/150‑1 (MQ=255)
                                                                                                                                      cgcgcATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAGAt         <  2:191177/150‑1 (MQ=255)
                                                                                                                                      cgcgcATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAGAt         >  2:790917/1‑150 (MQ=255)
                                                                                                                                          cATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAGATAAcc     >  1:528991/1‑150 (MQ=255)
                                                                                                                                             cAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAGATAACCTTg  <  2:1179989/150‑1 (MQ=255)
                                                                                                                                                   gATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAGATa        <  1:1210835/138‑1 (MQ=255)
                                                                                                                                                   gATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAGATa        >  2:1210835/1‑138 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GCCCCAGATCCCCACGGCGACGAAAATAACCGTCAACAGTGGTAATGAAATAACCCCAGCCCAGAAACGGGGAGAAATAACCCGACGCAGCGGATCCACCGCCATCATCTCCATACTGGAGAGTTGCTCTGTAGCGCGCATCAGGCCGATTTCTGCGGTTAGCGCCGAACCAGCACGCCCGGCAAACAACAACGCGGCAACCACCGGCCCCAGTTCACGCAGTAGCGATAACGCCACCAGCATACCCAGACTGGTTTCCGCACTATAAGTGGTCAGAACCAGATAACCTTG  >  NC_000913/3338742‑3339032

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑