Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,915 C→T 5.1% pseudogene (112/198 nt) ypaA → pseudogene, C‑terminal fragment, Transposase_31 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,9150CT 5.1% 265.3 / 3.1 79pseudogene (112/198 nt)ypaApseudogene, C‑terminal fragment, Transposase_31 family
Reads supporting (aligned to +/- strand):  ref base C (47/27);  new base T (2/2);  total (50/29)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATG  >  NC_000913/2357767‑2358060
                                                                                                                                                    |                                                                                                                                                 
ccAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgcg                                                                                                                                                  <  2:1639667/150‑1 (MQ=255)
 cAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgcgc                                                                                                                                                 >  2:1336372/1‑150 (MQ=255)
   aCTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCAt                                                                                                                                               <  2:437952/150‑1 (MQ=255)
    cTCACTGATGTTGATATTCCGAACAGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATg                                                                                                                                              <  2:1535842/150‑1 (MQ=255)
    cTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATg                                                                                                                                              <  2:602307/150‑1 (MQ=255)
         tGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGa                                                                                                                                         >  2:1007411/1‑150 (MQ=255)
          gATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGac                                                                                                                                        >  2:747396/1‑149 (MQ=255)
          gATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGaa                                                                                                                                        <  1:737962/150‑1 (MQ=255)
          gATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGaa                                                                                                                                        <  1:856387/150‑1 (MQ=255)
          gATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGaa                                                                                                                                        >  2:47061/1‑150 (MQ=255)
              cTGATATTCCGAACTGTCCTTAAAAGGAGAGATTAAGGACTATTGCTGACCGGCGGCGGCAGGAAGGACCTCAAATTGGCTGGCAGGAAGGTACCGTAGAAGGTTTGCATGACCAAGCCAGTAAAATTGCTTTGCGCCTGCTGGAACAgg                                                                                                                                    >  2:456715/1‑150 (MQ=25)
              cTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAgg                                                                                                                                    >  2:316369/1‑150 (MQ=255)
               tGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAggg                                                                                                                                   >  2:1073334/1‑150 (MQ=255)
                      ccGAACTGTCATTAAAATGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAATTTGGCTGGCAGTAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGAt                                                                                                                            <  2:1197884/150‑1 (MQ=255)
                      ccGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGAt                                                                                                                            <  1:1520665/150‑1 (MQ=255)
                              tCATTAAAAGGAGAGATTAATGACTATTGCTGTTCGGCTGCGGCAGGAAGGACCTCAAATTGGCTGGCAGGAAGGTATATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCa                                                                                                                    <  1:1626073/150‑1 (MQ=255)
                              tCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCa                                                                                                                    <  1:771609/150‑1 (MQ=255)
                              actacaaaaGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCa                                                                                                                    <  2:912501/145‑1 (MQ=255)
                                 ttAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAATGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGAGCATGCTGGAACAGGGCTTTGATCGTGACCAGGt                                                                                                                 >  2:552474/1‑150 (MQ=255)
                                 taataaGGAGAGATTAATGTCTAGTGCTGATCGGCTGCGGCAGGCTGGACATCACTTTGGGGGGCAGGAAGGTAACTTTGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGt                                                                                                                 <  2:57924/146‑1 (MQ=25)
                                    aaaGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCt                                                                                                              >  1:912163/1‑150 (MQ=255)
                                    aaaGGAGAGATAAATGACTAATGCAGAACGGCTGCGGCAGGACGGACATCAAATAGGCTGCCAGGAAGGTAAATTAGAAGGTTAGCATGAACAAGCCATAAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGAGACCACGAGCt                                                                                                              >  1:912180/1‑150 (MQ=25)
                                           agaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGACAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGcc                                                                                                       >  1:917268/1‑150 (MQ=255)
                                            gaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCa                                                                                                      >  1:1079014/1‑150 (MQ=255)
                                                  tGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTACAATTGCTTTGCGCATGCTGGAACAGGGCTTTGTTCGTGACCTGGTGCTCGCGGCCCCCCAGc                                                                                                >  1:1573421/1‑150 (MQ=255)
                                                  tGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGc                                                                                                >  1:1573413/1‑150 (MQ=255)
                                                   gACTATTGCTGAACGGCTTCGGCATGTTGGACATCAACTTGGCTGTCAGGAAGGTACCTTTGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCt                                                                                               <  2:1653075/150‑1 (MQ=25)
                                                        ttGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTTCGCATGCTGGAACAGTGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCg                                                                                          >  2:864382/1‑150 (MQ=255)
                                                        ttGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCg                                                                                          >  1:497954/1‑150 (MQ=255)
                                                          gCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGaa                                                                                        <  1:72927/150‑1 (MQ=255)
                                                              aaCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCg                                                                                    >  1:294350/1‑150 (MQ=255)
                                                               aCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGa                                                                                   >  1:824728/1‑150 (MQ=255)
                                                                   cTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTg                                                                               >  2:1488045/1‑150 (MQ=255)
                                                                    tGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTgg                                                                              <  2:109236/150‑1 (MQ=255)
                                                                     gCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGc                                                                             >  1:606112/1‑150 (MQ=255)
                                                                      cGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGTGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTaa                                                                                             >  1:1581907/1‑133 (MQ=255)
                                                                      cGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGTGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTaa                                                                                             >  1:1581901/1‑133 (MQ=255)
                                                                      cGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGTGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTaa                                                                                             <  2:1581901/133‑1 (MQ=255)
                                                                      cGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGTGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTaa                                                                                             <  2:1581907/133‑1 (MQ=255)
                                                                            gaaggaCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAAt                                                                      <  2:252091/150‑1 (MQ=255)
                                                                              aggaCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGACTaa                                                                    >  2:181701/1‑150 (MQ=255)
                                                                              aggaCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATaa                                                                    >  2:1416505/1‑150 (MQ=255)
                                                                              aggaCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATaa                                                                    >  2:1416506/1‑150 (MQ=255)
                                                                                 aCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCa                                                                 >  1:407468/1‑150 (MQ=255)
                                                                                              tggcAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAg                                                    >  2:1291115/1‑150 (MQ=255)
                                                                                                 gtggAAGGTTAATTAGCAGGTTTGCATGAACATGCCTTTATTAGTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACGAGGTGCTCGCGGCGACCCAGGTTAGCGATGCCGATCTGGCTGCGAATATCCTCTAATTAACACAGGcc                                                 <  2:912180/148‑1 (MQ=255)
                                                                                                 gaGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGcc                                                 <  2:912163/149‑1 (MQ=255)
                                                                                                   ggAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCca                                               <  1:792493/150‑1 (MQ=255)
                                                                                                    gAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCcac                                              >  1:514217/1‑150 (MQ=255)
                                                                                                    gAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCcac                                              >  2:434535/1‑150 (MQ=255)
                                                                                                             tttgAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATcc                                     <  2:536894/147‑1 (MQ=255)
                                                                                                             ttAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATcc                                     >  2:302063/1‑150 (MQ=255)
                                                                                                             ttAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCAACAGCCGATcc                                     >  2:1262124/1‑150 (MQ=255)
                                                                                                                    gTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATggg                              >  1:537025/1‑150 (MQ=255)
                                                                                                                    gTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACAAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATggg                              >  1:537049/1‑150 (MQ=255)
                                                                                                                      ttGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGcc                            <  1:662662/150‑1 (MQ=255)
                                                                                                                       tGCATGAACATGCCATTAAAATTGCTTTGCGCATTCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCTGCGAATAACCACTAATTATCACAGGCGCACAGCCGATGCCCATGGGCCt                           >  2:768626/1‑150 (MQ=38)
                                                                                                                       tGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGTCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCt                           >  2:1245749/1‑150 (MQ=255)
                                                                                                                         cATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCttt                         >  2:1493983/1‑150 (MQ=255)
                                                                                                                         cATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCttt                         >  2:646211/1‑150 (MQ=255)
                                                                                                                            gAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGat                      <  2:1696594/150‑1 (MQ=255)
                                                                                                                             aaCAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGata                     >  2:246927/1‑150 (MQ=255)
                                                                                                                               cAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATc                   <  1:829519/150‑1 (MQ=255)
                                                                                                                               cAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATc                   <  1:829517/150‑1 (MQ=255)
                                                                                                                                aaGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCa                  <  2:1619653/150‑1 (MQ=255)
                                                                                                                                aaGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCa                  <  2:1615025/150‑1 (MQ=255)
                                                                                                                                     aTTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAAc             >  1:534307/1‑150 (MQ=255)
                                                                                                                                     aTTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAAc             >  1:534276/1‑150 (MQ=255)
                                                                                                                                     aTTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAAc             <  2:1443750/150‑1 (MQ=255)
                                                                                                                                     aTTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAAc             <  1:969536/150‑1 (MQ=255)
                                                                                                                                      ttAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACt            >  2:1518504/1‑150 (MQ=255)
                                                                                                                                      ttAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACt            >  2:217831/1‑150 (MQ=255)
                                                                                                                                      ttAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACt            >  2:1518489/1‑150 (MQ=255)
                                                                                                                                      ttAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACt            >  2:1518459/1‑150 (MQ=255)
                                                                                                                                          aaTTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCt        >  1:200129/1‑150 (MQ=255)
                                                                                                                                                tttGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATg  >  2:26612/1‑150 (MQ=255)
                                                                                                                                                tttGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATg  >  2:1678652/1‑150 (MQ=255)
                                                                                                                                                tttGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATg  >  2:11997/1‑150 (MQ=255)
                                                                                                                                                tttGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATg  >  2:11993/1‑150 (MQ=255)
                                                                                                                                                tttGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATg  >  2:11978/1‑150 (MQ=255)
                                                                                                                                                tttGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATg  >  1:733114/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATG  >  NC_000913/2357767‑2358060

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: