Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,935,569 A→C 5.2% E94D (GAA→GAC rbsC → ribose ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,935,5690AC 5.2% 244.1 / 4.9 77E94D (GAA→GACrbsCribose ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (35/38);  new base C (2/2);  total (37/40)
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.15e-01

CTTATTCAATATTCTCCAGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTGCGCGGCGTGACCATGGTTT  >  NC_000913/3935422‑3935717
                                                                                                                                                   |                                                                                                                                                    
cTTATTCAATATTCTCCAGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGt                                                                                                                                                    >  1:299480/1‑150 (MQ=255)
 ttATTCAATATTCTCCAGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTc                                                                                                                                                   <  1:1046415/150‑1 (MQ=255)
  tATTAATTATTCTCCAGCAACCCTTAGTGAACGCCATTATGGCGGTCGGGTTGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCa                                                                                                                                                  <  1:714869/150‑1 (MQ=255)
      cAATATTCTCCAGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATgc                                                                                                                                              >  1:1508641/1‑150 (MQ=255)
      cAATATTCTCCAGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATgc                                                                                                                                              >  1:445157/1‑150 (MQ=255)
        atatTCTCCAGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATgcgc                                                                                                                                            <  1:737293/150‑1 (MQ=255)
         tatTCTCCAGCAAACCTCAGTGAACGCCATTATGGCGGTCGTGATGACGCTGGTGATCCTGACGTCGGTCATCGACTTATCGCTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTACATCTATCGTCGGCATTGACTTCAATgccga                                                                                                                                           <  1:176542/150‑4 (MQ=255)
          atTCTCCAGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGACTTCAATgccgac                                                                                                                                          <  1:292228/150‑5 (MQ=255)
           ttCTCCAGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCtgg                                                                                                                                         >  2:1054686/1‑150 (MQ=255)
           ttCTCCAGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCtgg                                                                                                                                         <  2:916165/150‑1 (MQ=255)
             ctcCAGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCtggtg                                                                                                                                       <  2:1605878/150‑1 (MQ=255)
              tcCGGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCtggtgg                                                                                                                                      <  1:725942/150‑1 (MQ=255)
              tcCAGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCtggtgg                                                                                                                                      <  1:295249/150‑1 (MQ=255)
               ccAGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGc                                                                                                                                     <  2:481278/150‑1 (MQ=255)
                 aGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGAATTGAAGTCAATGCGCTGGTGGCtt                                                                                                                                   >  1:1650748/1‑149 (MQ=255)
                  gCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGt                                                                                                                                  >  1:312674/1‑150 (MQ=255)
                       ccTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCgctg                                                                                                                             >  1:463522/1‑150 (MQ=255)
                       ccTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCgctg                                                                                                                             <  1:1277708/150‑1 (MQ=255)
                        cTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCgctgc                                                                                                                            <  2:1030668/150‑1 (MQ=255)
                          cAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCgctgctg                                                                                                                          >  2:474516/1‑150 (MQ=255)
                           aGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCgctgctgc                                                                                                                         <  2:4242/150‑1 (MQ=255)
                           aGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCgctgctgc                                                                                                                         <  2:10438/150‑1 (MQ=255)
                            gTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCgctgctgct                                                                                                                        <  1:1228803/150‑1 (MQ=255)
                            ccgaccGCCATAATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCAACTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCgctgctgct                                                                                                                        >  1:1471743/6‑150 (MQ=255)
                              gAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGaa                                                                                                                                                      <  1:296038/118‑1 (MQ=255)
                              gAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGaa                                                                                                                                                      >  2:296038/1‑118 (MQ=255)
                               aaCGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCtctc                                                                                                                     >  2:64442/1‑150 (MQ=255)
                                 cGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTcgc                                                                                                                   <  1:234968/150‑1 (MQ=255)
                                     aTTATGGCGTTCGGGATGAGGCTGGTGATCCTGACGTCGGGAATCGACTTATAGGTAGTTTATCTGTTGGCGCTGACCGGCTCATTTGCTGCATCTATCGTCGGGTTTGTAGTAAATGCGGTGGTGGCTGTCGCTGCTGCTCTCGCGTTa                                                                                                               <  2:1003850/150‑1 (MQ=255)
                                     aTTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTa                                                                                                               >  1:492857/1‑150 (MQ=255)
                                           gCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGcc                                                                                                         <  2:1382018/150‑1 (MQ=255)
                                           gCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTGGGCGCGGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGcc                                                                                                         <  2:1387716/150‑1 (MQ=255)
                                             ggTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGc                                                                                                       <  2:33800/150‑1 (MQ=255)
                                                 gggATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAAtt                                                                                                   <  2:299480/150‑1 (MQ=255)
                                                  ggATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTg                                                                                                  >  1:1213956/1‑150 (MQ=255)
                                                     tGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTGCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTg                                                                                               <  1:897127/150‑1 (MQ=255)
                                                      gACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGc                                                                                              <  2:1213956/150‑1 (MQ=255)
                                                        cGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCgg                                                                                            >  1:711207/1‑150 (MQ=255)
                                                          cTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTa                                                                                          <  1:1036396/150‑1 (MQ=255)
                                                           tGGTGATCCAGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTaa                                                                                         >  2:1185963/1‑150 (MQ=255)
                                                             gTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAAcc                                                                                       >  2:827027/1‑150 (MQ=255)
                                                                   cTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTACCCGGGGTg                                                                                 >  2:328243/1‑150 (MQ=255)
                                                                     gACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGAt                                                                               <  2:445157/150‑1 (MQ=255)
                                                                          cGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCGGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTat                                                                          <  2:969108/150‑2 (MQ=255)
                                                                                tCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAg                                                                    >  1:452224/1‑150 (MQ=255)
                                                                                    cTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTcg                                                                <  1:1043277/150‑1 (MQ=255)
                                                                                      tATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTcgcg                                                              >  1:689695/1‑150 (MQ=255)
                                                                                       aTCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGt                                                             >  2:1334626/1‑150 (MQ=255)
                                                                                         cGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTg                                                                             <  1:1007801/132‑1 (MQ=255)
                                                                                         cGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTg                                                                             >  2:1007801/1‑132 (MQ=255)
                                                                                         cGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTcc                                                           <  1:802317/150‑1 (MQ=255)
                                                                                          ggTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCgctgctgc                                                                                                                         <  1:462695/87‑1 (MQ=255)
                                                                                          ggTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCgctgctgc                                                                                                                         >  2:462695/1‑87 (MQ=255)
                                                                                          ggTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAg                                                                          >  2:748003/1‑134 (MQ=255)
                                                                                          ggTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAg                                                                          <  1:748003/134‑1 (MQ=255)
                                                                                                   tctGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTAt                                                 >  1:1105920/1‑150 (MQ=255)
                                                                                                            gcTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCt                                        <  1:64442/150‑1 (MQ=255)
                                                                                                              tGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTgg                                      >  1:46648/1‑150 (MQ=255)
                                                                                                              tGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTgg                                      <  2:1656279/150‑1 (MQ=255)
                                                                                                                 ccGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTa                                   >  2:168535/1‑150 (MQ=255)
                                                                                                                 ccGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTa                                   >  1:91678/1‑150 (MQ=255)
                                                                                                                   ggCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTatg                                 >  1:1107910/1‑150 (MQ=255)
                                                                                                                      gcAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGTTTatgatg                              <  1:63062/150‑1 (MQ=255)
                                                                                                                        aGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCt                            >  2:595302/1‑150 (MQ=255)
                                                                                                                         gTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCtt                           >  2:1553504/1‑150 (MQ=255)
                                                                                                                          ttGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTTGTGGCTGTCGCTGCTGATCTCGCGTTAGGTGCCGCAATTGGTGAGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCttt                          <  2:34168/150‑1 (MQ=255)
                                                                                                                          ttGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCttt                          <  1:977170/150‑1 (MQ=255)
                                                                                                                          ttGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCttt                          <  2:312674/150‑1 (MQ=255)
                                                                                                                             ctgcATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTAc                       >  1:1207780/1‑150 (MQ=255)
                                                                                                                               gcATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTg                     >  2:1336730/1‑150 (MQ=255)
                                                                                                                                  tCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTgcgc                  >  1:513112/1‑150 (MQ=255)
                                                                                                                                       cgtcgGCATTGAAGTCAATGCGCTGATGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTGCGCGGCGt             >  1:1122341/1‑150 (MQ=255)
                                                                                                                                       ccagcgCATTGACTTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTGCGCGGCGt             >  2:865652/6‑150 (MQ=255)
                                                                                                                                        gtcgGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTGCGCGGCGTg            <  2:338478/150‑1 (MQ=255)
                                                                                                                                         agcgCATTGACTTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTGCGCGGCGTGa           >  1:220633/4‑150 (MQ=255)
                                                                                                                                          cgGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTGCGCGGCGTGAc          <  2:711207/150‑1 (MQ=255)
                                                                                                                                           ggCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTGCGCGGCGTGAcc         <  2:1471743/150‑1 (MQ=255)
                                                                                                                                            gCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTGCGCGGCGTGACCa        >  1:1330581/1‑150 (MQ=255)
                                                                                                                                             cATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTGCGCGGCGTGACCAt       >  2:636426/1‑150 (MQ=255)
                                                                                                                                                  aaGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTGCGCGGCGTGACCATGGttt  <  2:689695/150‑1 (MQ=255)
                                                                                                                                                  aaGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTGCGCGGCGTGACCATGGttt  <  2:1107910/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTTATTCAATATTCTCCAGCAAACCTCAGTGAACGCCATTATGGCGGTCGGGATGACGCTGGTGATCCTGACGTCGGGCATCGACTTATCGGTAGGTTCTCTGTTGGCGCTGACCGGCGCAGTTGCTGCATCTATCGTCGGCATTGAAGTCAATGCGCTGGTGGCTGTCGCTGCTGCTCTCGCGTTAGGTGCCGCAATTGGTGCGGTAACCGGGGTGATTGTAGCGAAAGGTCGCGTCCAGGCGTTTATCGCTACGCTGGTTATGATGCTTTTACTGCGCGGCGTGACCATGGTTT  >  NC_000913/3935422‑3935717

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: