Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,382,115 C→T 5.6% P170S (CCG→TCG)  rbn → RNase BN, tRNA processing enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,382,1150CT 5.6% 281.0 / 4.7 89P170S (CCG→TCG) rbnRNase BN, tRNA processing enzyme
Reads supporting (aligned to +/- strand):  ref base C (49/35);  new base T (2/3);  total (51/38)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

CGCTGGCGAAATTCTCGATGATGGCCTGCGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGCGCTCGCTATTTTCGGCGATACCGGCCCCTGCGATGC  >  NC_000913/2381973‑2382260
                                                                                                                                              |                                                                                                                                                 
cgcTGGCGAAATTCTCGATGATGGCCTGCGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTg                                                                                                                                             >  8:479541/1‑149 (MQ=255)
  cTGGCGAAATTCTCGATGATGGCCTGCGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGc                                                                                                                                           <  4:372884/149‑1 (MQ=255)
       gAAATTCTCGATGATGGCCTGCGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACt                                                                                                                                      >  8:486970/1‑149 (MQ=255)
       gAAATTCTCGATGATGGCCTGCGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCACAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACt                                                                                                                                      >  8:297079/1‑149 (MQ=255)
             ctcGATGATGGCCTGCGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGttt                                                                                                                                  >  2:243766/1‑147 (MQ=255)
              tcGATGATGGCCTGCGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAg                                                                                                                               >  3:418013/1‑149 (MQ=255)
              tcGATGATGGCCTGCGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAg                                                                                                                               >  6:119463/1‑149 (MQ=255)
                gatgatGGCCTGCGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCa                                                                                                                                <  6:204483/146‑1 (MQ=255)
                 atgatgGCCTGAGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGaa                                                                                                                            <  3:263889/149‑1 (MQ=255)
                  tgatgGCCTGCGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGTATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTc                                                                                                                                 >  1:313274/1‑143 (MQ=255)
                  tgatgGCCTGCGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGTATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTc                                                                                                                                 <  2:313274/143‑1 (MQ=255)
                    atgGCCTGCGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAAtt                                                                                                                          <  1:428794/148‑1 (MQ=255)
                                 aGTAACCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGTCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCaaaa                                                                                                            <  2:104377/149‑1 (MQ=255)
                                     aCCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAAt                                                                                                        >  5:103006/1‑149 (MQ=255)
                                       cGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCa                                                                                                      >  4:364909/1‑149 (MQ=255)
                                       cGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCa                                                                                                      >  7:17764/1‑149 (MQ=255)
                                            aTCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGc                                                                                                   <  3:364909/147‑1 (MQ=255)
                                               cGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTgg                                                                                                >  4:436906/1‑147 (MQ=255)
                                                    gAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGAt                                                                                           <  7:486970/147‑1 (MQ=255)
                                                     aaCACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGa                                                                                        <  5:119463/149‑1 (MQ=255)
                                                      acacCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATgg                                                                                         >  1:385447/1‑147 (MQ=255)
                                                      acacCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGaa                                                                                       <  7:297079/149‑1 (MQ=255)
                                                       cacCCACTGGAATGTTATTGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGACCATGCATTAAAATCTGCTGGCGTCCCGCCTGGCACACTGTTTCAGGAATAAAAAGCGTGCAAAACAATCACGCTTGAAGATgtcag                                                                                      >  2:82564/1‑145 (MQ=255)
                                                       cacCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAg                                                                                      <  8:188690/149‑1 (MQ=255)
                                                          ccACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCa                                                                                   >  1:503571/1‑149 (MQ=255)
                                                           cACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAACCCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGTCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCACACCAATCACGCTGGAAGATGGAAgg                                                                                     >  6:429290/1‑146 (MQ=255)
                                                           cACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGTCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAgg                                                                                     >  6:429289/1‑146 (MQ=255)
                                                           cACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAgg                                                                                     >  8:226427/1‑146 (MQ=255)
                                                            aCTGTAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAgg                                                                                     >  2:290620/1‑145 (MQ=255)
                                                            aCTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGTCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGa                                                                                 <  5:429289/149‑1 (MQ=255)
                                                            aCTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGaa                                                                                                                            <  5:483714/106‑1 (MQ=255)
                                                            aCTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGaa                                                                                                                            >  6:483714/1‑106 (MQ=255)
                                                            aCTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGGGTCGCCTGGCCCACTGTTTCAGGAATTAAACGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGa                                                                                 <  5:429290/149‑1 (MQ=255)
                                                             cTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGAt                                                                                <  1:243766/149‑1 (MQ=255)
                                                              tGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGAt                                                                                >  8:23016/1‑148 (MQ=255)
                                                                gAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGAtt                                                                               <  3:50254/147‑1 (MQ=255)
                                                                  aTGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTaaa                                                                            <  1:290620/148‑2 (MQ=255)
                                                                  aTGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACg                                                                           <  2:503571/149‑1 (MQ=255)
                                                                       aTGGCTATAGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGgcg                                                                        >  8:140683/1‑147 (MQ=255)
                                                                        tGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGc                                                                                                   <  1:18102/119‑1 (MQ=255)
                                                                        tGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGc                                                                                                   >  2:18102/1‑119 (MQ=255)
                                                                           cTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAAtt                                                                                                                          <  3:54268/93‑1 (MQ=255)
                                                                           cTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAAtt                                                                                                                          >  4:54268/1‑93 (MQ=255)
                                                                              tCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCTCCTGGCCCACTGTTTCAGTAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACt                                                               <  6:199105/149‑1 (MQ=255)
                                                                              tCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACt                                                               >  1:411236/1‑149 (MQ=255)
                                                                                gTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACtt                                                              <  7:479541/148‑1 (MQ=255)
                                                                                gTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACt                                                               >  7:277323/1‑147 (MQ=255)
                                                                                 ttttGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTTGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTa                                                             <  1:427675/146‑1 (MQ=255)
                                                                                 tATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAg                                                            >  5:171671/1‑149 (MQ=255)
                                                                                 tATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTACAGGAATTAAAAGCGGGCAAAACAATCACGCTCGAAGATCGAAGGCAGATTAACGGCGAAGATTACt                                                               >  2:427675/1‑146 (MQ=255)
                                                                                     gaagaaCATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACaa                                                                                                         <  4:379734/100‑1 (MQ=255)
                                                                                     gaagaaCATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACaa                                                                                                         <  2:461160/100‑1 (MQ=255)
                                                                                     gaagaaCATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACaa                                                                                                         >  1:461160/1‑100 (MQ=255)
                                                                                     caagaaCATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACaa                                                                                                         >  3:379734/2‑100 (MQ=255)
                                                                                      aagaaCATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACg                                                                           <  4:432139/129‑1 (MQ=255)
                                                                                      aagaaCATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACg                                                                           >  3:432139/1‑129 (MQ=255)
                                                                                      aagaaCATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAgctgct                                                       >  2:69757/1‑149 (MQ=255)
                                                                                      aagaaCATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGAGTACTTAgctgct                                                       >  3:105607/1‑149 (MQ=255)
                                                                                        gaaCATGTTAAACCGGGTGCAATAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCCCCTGGCCCACTTTTTCATGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAg                                                                                      <  5:485436/116‑1 (MQ=255)
                                                                                           cATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGt                                                  >  3:229358/1‑149 (MQ=255)
                                                                                            aTGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTg                                                 >  6:252859/1‑149 (MQ=255)
                                                                                             tGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGc                                                >  4:15710/1‑149 (MQ=255)
                                                                                              gATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGc                                                >  4:40164/1‑148 (MQ=255)
                                                                                              gATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGc                                                >  8:440768/1‑148 (MQ=255)
                                                                                               aTAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGcc                                               <  8:17764/148‑1 (MQ=255)
                                                                                                tAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAg                                             >  4:338133/1‑149 (MQ=255)
                                                                                                 aaaCCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGgc                                                                         >  3:424070/1‑120 (MQ=255)
                                                                                                 aaaCCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGgc                                                                         <  4:424070/120‑1 (MQ=255)
                                                                                                           cATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCgg                                                                                                                  >  4:186996/1‑69 (MQ=255)
                                                                                                           cATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCgg                                                                                                                  <  3:186996/69‑1 (MQ=255)
                                                                                                           cATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGcgct                                   <  1:69757/148‑1 (MQ=255)
                                                                                                            aTTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGcgctcg                                 >  3:312604/1‑149 (MQ=255)
                                                                                                               aaaTGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGcgctcgct                               <  5:252859/148‑1 (MQ=255)
                                                                                                                    cccAGGCATTAAAAGCTTCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACCATCCCGCTGGAAGATGGAAGGCAGATTAACGCCGCAGAGTACTTCGCTTCTCCAGTGCCAGGTAAAGCGCTCGCTATTTTc                         <  1:82564/149‑1 (MQ=255)
                                                                                                                     ccAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAgc                                                           <  5:186197/114‑1 (MQ=255)
                                                                                                                     ccAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAgc                                                           >  6:186197/1‑114 (MQ=255)
                                                                                                                     ccAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGCGCTCGCTATTTTCg                        >  7:140602/1‑149 (MQ=255)
                                                                                                                      cAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAGAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGCGCTCGCTATTTTCgg                       >  7:84157/1‑149 (MQ=255)
                                                                                                                      cAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGCGCTCGCTATTTTCgg                       <  8:277323/149‑1 (MQ=255)
                                                                                                                           aTTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGCGCTCGCTATTTTCGGCGATa                  >  6:450784/1‑149 (MQ=255)
                                                                                                                            ttAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGCGCTCGCTATTTTCGGCGATAc                 >  2:493551/1‑149 (MQ=255)
                                                                                                                              aaaaGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTaaa                                        >  7:488948/1‑124 (MQ=255)
                                                                                                                              aaaaGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTaaa                                        <  8:488948/124‑1 (MQ=255)
                                                                                                                               aaaGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGCGCTCGCTATTTTCGGCGATACCgg              >  4:329928/1‑149 (MQ=255)
                                                                                                                                    tgctgGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAgg                                            <  3:220546/114‑1 (MQ=255)
                                                                                                                                    tgctgGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAgg                                            >  4:220546/1‑114 (MQ=255)
                                                                                                                                    tgctgGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGCGCTCGCTATTTTCGGCGATACCGGCCCCt         >  4:309957/1‑149 (MQ=255)
                                                                                                                                    tgctgGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGCGCTCGCTATTTTCGGCGATACCGGCCCCt         >  3:72683/1‑149 (MQ=255)
                                                                                                                                           gTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGCGCTCGCTATTTTCGGCGATACCGGCCCCTGCGATGc  <  3:40164/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
CGCTGGCGAAATTCTCGATGATGGCCTGCGCAAAGTAACCGCTTATCCGCTGGAACACCCACTGGAATGTTATGGCTATCGTATTGAAGAACATGATAAACCGGGTGCATTAAATGCCCAGGCATTAAAAGCTGCTGGCGTGCCGCCTGGCCCACTGTTTCAGGAATTAAAAGCGGGCAAAACAATCACGCTGGAAGATGGAAGGCAGATTAACGGCGCAGATTACTTAGCTGCTCCAGTGCCAGGTAAAGCGCTCGCTATTTTCGGCGATACCGGCCCCTGCGATGC  >  NC_000913/2381973‑2382260

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