Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,970 A→G 7.5% intergenic (+368/+246) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,269,9700AG 7.5% 215.2 / 7.3 81intergenic (+368/+246)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (33/41);  new base G (3/3);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGT  >  NC_000913/3269826‑3270106
                                                                                                                                                |                                                                                                                                        
aaaCCGCACTCAACAATCTGATACCCGCCGGAACGGCTTATTACCTACATCAGGTAAAACTGACCGACAAGCCGCTAGCTTTTGGGTATAGTGGCGGGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                      <  4:356005/149‑1 (MQ=18)
  aGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATATTGTCGTGGACCGGCAGTGATCTTGAGTCCCCTACAAAAGTAAAAACACGCCGAAGCGGGt                                                                                                                                    >  1:354275/1‑149 (MQ=18)
       aCTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                               <  8:413446/149‑1 (MQ=255)
        cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                              <  4:184133/149‑1 (MQ=255)
        cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                              <  5:54849/149‑1 (MQ=255)
        cTCAACAAGCTCAAACCCGCCGCACCGGTTTATTACGTAAAACAGGTTAAACGGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGTGGGTTTTTAc                                                                                                                              <  5:290055/149‑1 (MQ=11)
           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAACCCCGCCGAAGCGGGTTTTTACGTa                                                                                                                           <  6:431218/149‑1 (MQ=255)
           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                      >  2:330262/1‑138 (MQ=255)
           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                      <  1:330262/138‑1 (MQ=255)
           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                            <  5:242530/148‑1 (MQ=255)
           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                            >  6:388357/1‑148 (MQ=255)
           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATACTTCTGGCCCACCAATAGAAAAAACACGCGGAAGCGGTttttgtatt                                                                                                                            >  1:152699/1‑143 (MQ=17)
             caaTCTGAAACCCTCCGTAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGATTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                      <  3:274282/136‑1 (MQ=255)
              aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                        <  1:294231/149‑1 (MQ=255)
              aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                        <  5:388357/149‑1 (MQ=255)
              aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                        >  7:108007/1‑149 (MQ=255)
              aaTCTCAAACCCGCCGCACCGGTTTATTACGTACAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaac                                                                                                                        <  1:423112/149‑2 (MQ=11)
                  tGAAACCCGCGGGAACGGTTTATTAGGTACATCAGGTAAAAGTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGACCCTCCTAAAGTAAAAACCCGCGGATGCGGGTTTTTACGTAAAACAg                                                                                                                     <  2:267426/148‑1 (MQ=17)
                  tGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                      >  1:467647/1‑131 (MQ=255)
                  tGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                      <  2:467647/131‑1 (MQ=255)
                    aaaCCCTCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCTATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                                                                                                    <  5:325040/147‑1 (MQ=255)
                    aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                                                                                                    <  1:391813/147‑1 (MQ=255)
                    aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                                                                                                    <  5:105352/147‑1 (MQ=255)
                       cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                               <  2:321389/149‑1 (MQ=255)
                       cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                               <  4:397537/149‑1 (MQ=255)
                        ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAAc                                                                                                                       <  6:436040/140‑1 (MQ=255)
                        ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACaagatcgt                                                                                                               >  5:436040/1‑141 (MQ=255)
                        ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACATGGGaa                                                                                                                >  1:179747/1‑147 (MQ=255)
                            cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATCAGCCGCTTTCTTATGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                          <  5:182939/149‑1 (MQ=255)
                            cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                          <  8:108007/149‑1 (MQ=255)
                                                   cggTACAACTGACATATAACCCGTTTTCTTATGCTTTTAGTGTAGTGGACAGTCATTCTATTTTCTGGCCCTCAAAAATTATAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTggg                                                                                    <  2:354275/147‑1 (MQ=2)
                                                                cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                                                                                                      >  3:334845/1‑101 (MQ=35)
                                                                cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                                                                                                      <  4:334845/101‑1 (MQ=35)
                                                                cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                   >  1:166135/1‑104 (MQ=35)
                                                                cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                   >  1:172595/1‑104 (MQ=35)
                                                                cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                   <  2:172595/104‑1 (MQ=35)
                                                                cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                   <  2:166135/104‑1 (MQ=35)
                                                                      gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTTACCGATAAGCCGctttct                                                                                          >  1:201596/1‑123 (MQ=25)
                                                                      gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                                                                                          <  2:201596/123‑1 (MQ=35)
                                                                               ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                               <  1:202162/93‑1 (MQ=37)
                                                                               ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                               >  2:202162/1‑93 (MQ=37)
                                                                                      tataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                       >  7:451626/1‑110 (MQ=35)
                                                                                      tataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                       <  8:451626/110‑1 (MQ=35)
                                                                                          gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                         <  4:88839/120‑1 (MQ=39)
                                                                                          gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                         >  3:88839/1‑120 (MQ=39)
                                                                                          aaaagaaagcggcttatcggtcagtttcacctgttttacgtaaAAACCCGCTTCGGCGGGCTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                         <  1:213839‑M2/77‑1 (MQ=255)
                                                                                          aaaagaaagcggcttatcggtcagtttcacctgttttacgtaaAAACCCGCTTCGGCGGGCTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                         >  2:213839‑M2/44‑120 (MQ=255)
                                                                                                tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                   <  6:312712/72‑1 (MQ=39)
                                                                                                tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                   >  5:312712/1‑72 (MQ=39)
                                                                                                  gACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                     >  4:315440/1‑100 (MQ=18)
                                                                                                       tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                                         <  8:246501/107‑1 (MQ=17)
                                                                                                       tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                                         >  7:246501/1‑107 (MQ=17)
                                                                                                       tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTg                                                                        >  6:387507/1‑108 (MQ=40)
                                                                                                       tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACACGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTg                                                                        <  5:387507/108‑1 (MQ=255)
                                                                                                         attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattc                                                             >  5:32709/1‑117 (MQ=17)
                                                                                                         attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattc                                                             <  6:32709/117‑1 (MQ=17)
                                                                                                         atcggtcagtttcacctggtttacgtaaAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                                                                                         >  4:305286‑M2/29‑73 (MQ=255)
                                                                                                         atcggtcagtttcacctggtttacgtaaAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                                                                                         <  3:305286‑M2/45‑1 (MQ=255)
                                                                                                            catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGcttt                                                                                            >  5:416761/1‑83 (MQ=255)
                                                                                                             atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctct                                                                                            >  2:211758/1‑80 (MQ=255)
                                                                                                             atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGGCGGGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAc                                 <  3:325904/141‑1 (MQ=255)
                                                                                                             atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattca                                                            >  8:124020/1‑114 (MQ=17)
                                                                                                             atCTTTCTGCCCCTCCAAAAGAAAAAACCAGCACAAGCGAGGTTATACGTAAAACAGGTTAAACAGACCGATAAGCCGCTTTCTTTTGGGGATAGTGTATTGTAAAATAATTCGTCTTTGTGGCCCTGCAAAAGGAAAGAccagatcgg                         >  4:325904/1‑142 (MQ=2)
                                                                                                                   cTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTg                                                                        >  1:252048/1‑96 (MQ=40)
                                                                                                                    tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                     <  6:366344/82‑1 (MQ=38)
                                                                                                                      ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                                         >  7:14919/1‑92 (MQ=18)
                                                                                                                      ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTg                                                                        >  3:206/1‑93 (MQ=40)
                                                                                                                      ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTg                                                                        <  4:206/93‑1 (MQ=40)
                                                                                                                      acctgttttacgtaaAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTcc                                             <  7:58777‑M2/105‑1 (MQ=255)
                                                                                                                      acctgttttacgtaaAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTcc                                             >  8:58777‑M2/16‑120 (MQ=255)
                                                                                                                           ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                      <  6:432817/58‑1 (MQ=255)
                                                                                                                           ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                      >  5:432817/1‑58 (MQ=255)
                                                                                                                           ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                   <  6:296684/109‑1 (MQ=40)
                                                                                                                           ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                   >  5:296684/1‑109 (MQ=40)
                                                                                                                             aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                                                                                                  >  8:299683/1‑60 (MQ=255)
                                                                                                                                 gtaaAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGCCATTCATCTTTCTGc                                                   >  7:31837‑M2/5‑103 (MQ=255)
                                                                                                                                 gtaaAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGCCATTCATCTTTCTGc                                                   <  8:31837‑M2/99‑1 (MQ=255)
                                                                                                                                  taaAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                                                                                        >  1:78640‑M2/4‑65 (MQ=255)
                                                                                                                                  taaAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                                                                                        <  2:78640‑M2/62‑1 (MQ=255)
                                                                                                                                  taaAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTAAACCAg    >  6:257232‑M2/4‑149 (MQ=255)
                                                                                                                                    aAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGt  <  5:257232‑M2/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                        
AAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGT  >  NC_000913/3269826‑3270106

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