Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,994,371 T→C 8.1% pseudogene (433/503 nt) pbl → transglycosylase domain‑containing protein Pbl

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,994,3710TC 8.1% 252.4 / 8.2 86pseudogene (433/503 nt)pbltransglycosylase domain‑containing protein Pbl
Reads supporting (aligned to +/- strand):  ref base T (39/40);  new base C (2/5);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

TTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGT  >  NC_000913/2994273‑2994468
                                                                                                  |                                                                                                 
ttttATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAat                                                                                                   >  1:476451/1‑99 (MQ=255)
ttttATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAat                                                                                                   >  2:1101225/1‑99 (MQ=255)
 tttATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAata                                                                                                  <  1:362847/99‑1 (MQ=255)
 tttATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAata                                                                                                  <  2:663862/99‑1 (MQ=255)
 tttATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAata                                                                                                  <  2:2037931/99‑1 (MQ=255)
 tttATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAata                                                                                                  <  2:1851891/99‑1 (MQ=255)
   tatataACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATg                                                                                                <  1:1222027/99‑1 (MQ=255)
   tatataACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATg                                                                                                <  2:233830/99‑1 (MQ=255)
   tatataACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATg                                                                                                <  1:602204/99‑1 (MQ=255)
   tatataACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATg                                                                                                <  1:733880/99‑1 (MQ=255)
    atataACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAACATGc                                                                                               <  2:599956/99‑1 (MQ=255)
     tataACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGcc                                                                                              <  2:595352/99‑1 (MQ=255)
     tataACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGcc                                                                                              <  1:1779573/99‑1 (MQ=255)
     tataACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGcc                                                                                              >  1:1612857/1‑99 (MQ=255)
      ataACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGccc                                                                                             <  2:995329/99‑1 (MQ=255)
      ataACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGccc                                                                                             <  2:146197/99‑1 (MQ=255)
       taACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGcccc                                                                                            >  1:820083/1‑99 (MQ=255)
        aaCCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCa                                                                                           >  1:100429/1‑99 (MQ=255)
         aCCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCaa                                                                                          >  1:751495/1‑99 (MQ=255)
         aCCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCaa                                                                                          >  2:1248377/1‑99 (MQ=255)
           cATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCaaaa                                                                                        >  2:841101/1‑99 (MQ=255)
           cATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCaaaa                                                                                        >  2:1029995/1‑99 (MQ=255)
           cATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCaaaa                                                                                        >  1:1425288/1‑99 (MQ=255)
           cATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAACATGCCCCaaaa                                                                                        <  2:90188/99‑1 (MQ=255)
            aTTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAat                                                                                       <  2:1540454/99‑1 (MQ=255)
             tttttCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAata                                                                                      >  1:1575089/1‑99 (MQ=255)
              ttttCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAatat                                                                                     >  1:1240991/1‑99 (MQ=255)
                  cACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAatatatat                                                                                 <  1:1613155/99‑1 (MQ=255)
                   aCGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATAtt                                                                                >  1:1091488/1‑99 (MQ=255)
                   aCGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATAtt                                                                                <  2:1383201/99‑1 (MQ=255)
                   aCGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATAtt                                                                                <  1:870262/99‑1 (MQ=255)
                   aCGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATAtt                                                                                >  1:680344/1‑99 (MQ=255)
                    cGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTg                                                                               <  1:344263/99‑1 (MQ=255)
                     gTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGt                                                                              >  1:643068/1‑99 (MQ=255)
                      ttGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGta                                                                             <  2:512890/99‑1 (MQ=255)
                      ttGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGta                                                                             >  2:1509902/1‑99 (MQ=255)
                      ttGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGta                                                                             <  2:1466478/99‑1 (MQ=255)
                        gtgGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGtata                                                                           >  1:218264/1‑99 (MQ=255)
                         tgGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGtatat                                                                          <  1:795340/99‑1 (MQ=255)
                              tAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAg                                                                     >  2:880293/1‑99 (MQ=255)
                              tAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAg                                                                     >  2:445341/1‑99 (MQ=255)
                                aCGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGc                                                                   >  2:649377/1‑99 (MQ=255)
                                  gTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCtt                                                                 <  1:630966/99‑1 (MQ=255)
                                  gTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAACATGCCCCAAAATATATATTGTATATACCAGGCtt                                                                 >  2:623417/1‑99 (MQ=255)
                                   tGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTa                                                                >  1:697688/1‑99 (MQ=255)
                                   tGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTa                                                                <  1:1612674/99‑1 (MQ=255)
                                    ggCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTaa                                                               <  2:1220572/99‑1 (MQ=255)
                                     gCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAAt                                                              <  1:1495113/99‑1 (MQ=255)
                                       aaTGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGa                                                            <  1:237075/99‑1 (MQ=255)
                                        aTGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGaa                                                           <  1:1880887/99‑1 (MQ=255)
                                          gCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACt                                                         >  1:1835069/1‑99 (MQ=255)
                                          gCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACt                                                         <  1:383559/99‑1 (MQ=255)
                                                ggACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACaa                                                   <  1:73552/99‑1 (MQ=255)
                                                ggACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACaa                                                   >  1:1496497/1‑99 (MQ=255)
                                                 gACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAAt                                                  <  2:408789/99‑1 (MQ=255)
                                                 gACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAAt                                                  >  1:1334954/1‑99 (MQ=255)
                                                   cTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCg                                                <  2:770142/99‑1 (MQ=255)
                                                       taATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAg                                            >  1:270038/1‑99 (MQ=255)
                                                         aTGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGc                                          <  2:1347690/99‑1 (MQ=255)
                                                          tGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCt                                         <  1:141303/99‑1 (MQ=255)
                                                            cGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTc                                       <  1:1014980/99‑1 (MQ=255)
                                                            cGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTc                                       >  1:1587475/1‑99 (MQ=255)
                                                             gggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAACATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCa                                      <  2:1285030/99‑1 (MQ=255)
                                                              ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAg                                     <  1:687803/99‑1 (MQ=255)
                                                               gTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGc                                    >  2:422239/1‑99 (MQ=255)
                                                                   tGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCaaaa                                <  2:1147720/99‑1 (MQ=255)
                                                                   tGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCaaaa                                <  2:906936/99‑1 (MQ=255)
                                                                       aTGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAg                            >  2:1537818/1‑99 (MQ=255)
                                                                       aTGGATAACCAAAAAAAGAGACAGCAATACGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAg                            >  1:1054715/1‑99 (MQ=255)
                                                                        tGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGc                           >  2:209810/1‑99 (MQ=255)
                                                                        tGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGc                           <  1:1705774/99‑1 (MQ=255)
                                                                            tAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTaa                       >  1:419178/1‑99 (MQ=255)
                                                                            tAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTaa                       >  2:983577/1‑99 (MQ=255)
                                                                             aaCCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTaaa                      <  2:1788419/99‑1 (MQ=255)
                                                                              aCCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAg                     >  2:460385/1‑99 (MQ=255)
                                                                               ccAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACATTCGTAAGGCTTCAGCAAAATGAGCCTAAAGc                    >  1:285987/1‑99 (MQ=255)
                                                                               ccAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGc                    >  2:730928/1‑99 (MQ=255)
                                                                                  aaaaaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc                 >  1:1908651/1‑99 (MQ=255)
                                                                                  aaaaaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc                 >  2:1529981/1‑99 (MQ=255)
                                                                                    aaaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtt               >  2:304150/1‑99 (MQ=255)
                                                                                    aaaaGAGACAGCAACATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtt               >  1:405441/1‑99 (MQ=255)
                                                                                     aaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCttt              <  2:891922/99‑1 (MQ=255)
                                                                                      aaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGGGCCTAAAGCCTCtttt             >  2:775259/1‑99 (MQ=255)
                                                                                      aaGAGACAGCAACATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttt             <  1:1942608/99‑1 (MQ=255)
                                                                                        gagaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttt           <  2:905250/99‑1 (MQ=255)
                                                                                        gagaCAGCAACATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttt           <  1:909232/99‑1 (MQ=255)
                                                                                         agaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCttttttt          >  2:1367335/1‑99 (MQ=255)
                                                                                         agaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCttttttt          <  2:93774/99‑1 (MQ=255)
                                                                                             aGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtata      <  2:262876/99‑1 (MQ=255)
                                                                                             aGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtata      <  2:1478762/99‑1 (MQ=255)
                                                                                               cAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtatata    <  1:1806076/99‑1 (MQ=255)
                                                                                                aaTATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAg   >  2:214227/1‑99 (MQ=255)
                                                                                                 atatGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGt  >  2:1110656/1‑99 (MQ=255)
                                                                                                  |                                                                                                 
TTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGT  >  NC_000913/2994273‑2994468

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