Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,712,130 T→A 8.3% Q212L (CAG→CTG)  yfiC ← tRNA m(6)A37 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,712,1300TA 8.3% 252.1 / 10.9 84Q212L (CAG→CTG) yfiCtRNA m(6)A37 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (41/36);  new base A (4/3);  total (45/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 = 5.35e-01

TGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCC  >  NC_000913/2711984‑2712279
                                                                                                                                                  |                                                                                                                                                     
tGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCg                                                                                                                                                    >  2:1397012/1‑150 (MQ=255)
 gCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTAGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCgg                                                                                                                                                   >  1:603245/1‑150 (MQ=255)
  cTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCggg                                                                                                                                                  <  1:622825/150‑1 (MQ=255)
   tCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCgggg                                                                                                                                                 <  2:286340/150‑1 (MQ=255)
    cACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGga                                                                                                                                                >  2:179629/1‑150 (MQ=255)
     aCCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGgag                                                                                                                                               <  1:1334515/150‑1 (MQ=255)
          ttCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATGc                                                                                                                                          <  1:1261382/150‑1 (MQ=255)
          ttCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATGc                                                                                                                                          <  1:48252/150‑1 (MQ=255)
              ggTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                      >  2:767941/1‑150 (MQ=255)
               gTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATGctttag                                                                                                                                     <  1:856996/150‑6 (MQ=255)
                tCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagcag                                                                                                                                    >  1:406183/1‑150 (MQ=255)
                      gATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagca                                                                                                                                     <  2:325231/143‑1 (MQ=255)
                      gATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagca                                                                                                                                     >  1:325231/1‑143 (MQ=255)
                        tCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGAt                                                                                                                            <  1:197633/150‑1 (MQ=255)
                        tCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGAt                                                                                                                            <  1:892977/150‑1 (MQ=255)
                         ccTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATg                                                                                                                           <  1:1127781/150‑1 (MQ=255)
                            ttCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGcc                                                                                                                                         >  2:375103/1‑133 (MQ=255)
                            ttCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGcc                                                                                                                                         <  1:375103/133‑1 (MQ=255)
                            ttCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCgg                                                                                                                        <  1:808007/150‑1 (MQ=255)
                             tCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGc                                                                                                                       <  1:395094/150‑1 (MQ=255)
                             tCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGc                                                                                                                       <  1:502172/150‑1 (MQ=255)
                                cccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGt                                                                                                                    >  1:517407/1‑150 (MQ=255)
                                 ccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTc                                                                                                                   >  1:836539/1‑150 (MQ=255)
                                      aaTCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCt                                                                                                              <  2:1123436/150‑1 (MQ=255)
                                                tAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCg                                                                                                    >  2:948973/1‑150 (MQ=255)
                                                    cAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCcac                                                                                                <  1:319632/150‑1 (MQ=255)
                                                    cAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCcac                                                                                                <  1:948973/150‑1 (MQ=255)
                                                        tAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTAGCGCTTCGTTTTCCGCCACAtcc                                                                                            >  2:690734/1‑150 (MQ=255)
                                                         aaaaaGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtccg                                                                                           >  2:1012066/1‑150 (MQ=255)
                                                          aaaaGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtccgt                                                                                          >  2:584556/1‑150 (MQ=255)
                                                          aaaaGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCACCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtccgg                                                                                          >  2:968778/1‑149 (MQ=255)
                                                           aaaGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtccgtc                                                                                         >  1:360024/1‑150 (MQ=255)
                                                              gCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCACCACATCCGTCCGc                                                                                      <  1:936856/150‑1 (MQ=255)
                                                               ccTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGCCTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCCGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCGCATCCGTCCGCa                                                                                     >  2:1268289/1‑150 (MQ=255)
                                                               ccTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCa                                                                                     <  2:1109242/150‑1 (MQ=255)
                                                                cTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGCATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCaa                                                                                    >  2:1170487/1‑150 (MQ=255)
                                                                    gTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGt                                                                                >  2:254797/1‑150 (MQ=255)
                                                                        gcgcCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATTACTAAGCGATCGCTAAAGCATTCTCCCGCATGCGGGGAGAATGACAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTATCCGCCACCAACGTCCGCAAACGTAAAt                                                                            >  2:652175/1‑150 (MQ=255)
                                                                        gcgcCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAAt                                                                            >  2:916127/1‑150 (MQ=255)
                                                                           ccGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGcc                                                                         <  1:78865/150‑1 (MQ=255)
                                                                             gTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAACCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGGTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAg                                                                       <  2:428394/150‑1 (MQ=255)
                                                                              tataCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCaa                                                                                    >  1:588910/1‑136 (MQ=255)
                                                                              tataCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCaa                                                                                    <  2:588910/136‑1 (MQ=255)
                                                                              tataCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGc                                                                      >  2:498545/1‑150 (MQ=255)
                                                                              tataCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGc                                                                      >  2:807286/1‑150 (MQ=255)
                                                                              tataCGCGTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGCATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCCGc                                                                      >  1:377639/1‑150 (MQ=255)
                                                                               ataCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGcc                                                                     <  1:1368484/150‑1 (MQ=255)
                                                                                 aCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCa                                                                   >  1:623888/1‑150 (MQ=255)
                                                                                   gCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATg                                                                 >  1:377698/1‑150 (MQ=255)
                                                                                     ttCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCt                                                               <  2:517407/150‑1 (MQ=255)
                                                                                      tCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCtt                                                              >  1:954335/1‑150 (MQ=255)
                                                                                       cGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCTCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGTTGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTa                                                             <  2:915818/150‑1 (MQ=255)
                                                                                       cGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTa                                                             <  1:821697/150‑1 (MQ=255)
                                                                                        ggAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTaa                                                            <  1:632215/150‑1 (MQ=255)
                                                                                            tAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGcc                                                        <  2:406183/150‑1 (MQ=255)
                                                                                                tctggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGGAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCt                                                               <  2:155349/139‑1 (MQ=255)
                                                                                                tctggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCt                                                               >  1:155349/1‑139 (MQ=255)
                                                                                                 ctggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTc                                                   >  2:378530/1‑150 (MQ=255)
                                                                                                   ggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCg                                                 >  2:1273058/1‑150 (MQ=255)
                                                                                                     tctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTa                                               <  1:1012066/150‑1 (MQ=255)
                                                                                                        ggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAAc                                           <  1:641888/150‑1 (MQ=255)
                                                                                                        ggCCCACGAATGACTAAGCGATCGCTAAAGCATACTCACGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGATTAATGCCAGCTCCGTaaac                                            >  1:970527/1‑149 (MQ=255)
                                                                                                             aCGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGACAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCAtt                                       <  2:759399/150‑1 (MQ=255)
                                                                                                                aaTGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTAcc                                    <  2:214718/150‑1 (MQ=255)
                                                                                                                aaTGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTAcc                                    <  1:254797/150‑1 (MQ=255)
                                                                                                                           gATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTg                         >  2:300804/1‑150 (MQ=255)
                                                                                                                           gATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTg                         >  1:334777/1‑150 (MQ=255)
                                                                                                                            aTCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTgg                        >  2:326899/1‑150 (MQ=255)
                                                                                                                             tCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgc                       <  1:300804/150‑1 (MQ=255)
                                                                                                                              cGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgca                      >  1:983144/1‑150 (MQ=255)
                                                                                                                              cGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCACATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgca                      >  1:815977/1‑150 (MQ=255)
                                                                                                                                cTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcagc                    >  1:843452/1‑150 (MQ=255)
                                                                                                                                cTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcagc                    <  1:179629/150‑1 (MQ=255)
                                                                                                                                  aaaGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCAc                  <  1:767941/150‑1 (MQ=255)
                                                                                                                                    ggCATTCTCCCCGCAGGCGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTa                >  1:1078142/2‑150 (MQ=255)
                                                                                                                                     gCATTCTCCCCGCAGGCGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTAc               >  1:556264/1‑150 (MQ=255)
                                                                                                                                      cATTCTCCCCGCAGGCGGGAGAATGCCAGCAGCACCCGATTCTGCAGTCTCGCTTAGTTTTAAGCCACATCCGTCCGCAAACGAAAATGCCAGACCATGATTAATGCCATCTACCTAAAACCATTACCAATCTGCACTGGAAGCACTaag              >  2:1056152/1‑148 (MQ=255)
                                                                                                                                      cATTCTCCCCGCAGGCGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACg              >  2:299758/1‑150 (MQ=255)
                                                                                                                                       aTTCTCCCGCCTGCGGGGAGAATGCCAGCATCACCCGATGCGGCAGGCGCGCTTCGTTTTACGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCTTTACCAATCTGCTCTGTCAGCACTACGc             <  2:377639/150‑1 (MQ=255)
                                                                                                                                       aTTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGc             >  1:991592/1‑150 (MQ=255)
                                                                                                                                          ctcCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGa          >  1:1228692/1‑150 (MQ=255)
                                                                                                                                            cccGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGaaa        >  1:689186/1‑150 (MQ=255)
                                                                                                                                             ccGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGaaaa       >  2:506124/1‑150 (MQ=255)
                                                                                                                                                ccTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATc    <  2:865849/150‑1 (MQ=255)
                                                                                                                                                  tGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATccc  >  1:372616/1‑150 (MQ=255)
                                                                                                                                                  tGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATccc  <  1:916127/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCC  >  NC_000913/2711984‑2712279

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: