Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,991,754 (A)7→6 intergenic (‑135/+61) yecA ← / ← leuZ UPF0149 family protein YecA/tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,991,7480A.99.6% 1296.2 / ‑2.6 243intergenic (‑129/+67)yecA/leuZUPF0149 family protein YecA/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (168/74);  minor base T (1/0);  total (169/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 279 total aligned reads displayed.

CAGCACCTATCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGACTTGAACCCGCA  >  NC_000913/1991648‑1991843
                                                                                                    |                                                                                               
cAGCACCTATCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGcaca                                                                                                 >  2:3643993/1‑101 (MQ=255)
  gCACCTATCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACaaa                                                                                               <  2:3643386/101‑1 (MQ=255)
   cACCTATCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACaaaa                                                                                              >  2:3643998/1‑101 (MQ=255)
   cACCTATCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACaaaa                                                                                              <  1:3643505/101‑1 (MQ=255)
   cACCTATCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACaaaa                                                                                              <  2:3643491/101‑1 (MQ=255)
    aCCTATCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACaaaaa                                                                                             >  1:3644000/1‑101 (MQ=255)
     ccTATCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACaaaaaa                                                                                            <  1:3643450/101‑1 (MQ=255)
     ccTATCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGAACaaaaaa                                                                                            <  1:3643581/101‑1 (MQ=255)
        aTCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACaaaaaacca                                                                                         <  2:3643686/101‑4 (MQ=255)
        aTCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACaaaaaacca                                                                                         <  2:3643535/101‑4 (MQ=255)
         tCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACaaaaaaccac                                                                                        <  1:3643324/101‑5 (MQ=255)
         tCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACaaaaaaccac                                                                                        >  2:3644005/1‑97 (MQ=255)
         tCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACaaaaaaccac                                                                                        <  2:3643870/101‑5 (MQ=255)
         tCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACaaaaaaccac                                                                                        <  2:3643786/101‑5 (MQ=255)
             aCATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGa                                                                                   >  2:3644008/1‑101 (MQ=255)
              cATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGaa                                                                                  <  2:3643508/101‑1 (MQ=255)
              cATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGaa                                                                                  >  1:3644009/1‑101 (MQ=255)
                tGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAgg                                                                                >  2:3644014/1‑101 (MQ=255)
                   cGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAggtgg                                                                             >  1:3644016/1‑101 (MQ=255)
                    ggTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGaa                                                                                  >  2:3644020/1‑95 (MQ=255)
                    ggTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAggtggt                                                                            >  1:3644018/1‑101 (MQ=255)
                        cGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCa                                                                        >  2:3644024/1‑101 (MQ=255)
                         gTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCAc                                                                       >  1:3644026/1‑101 (MQ=255)
                          tATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACg                                                                      >  2:3644029/1‑101 (MQ=255)
                             cAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGaca                                                                   >  1:3644032/1‑101 (MQ=255)
                              agagATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGacac                                                                  >  2:3644036/1‑101 (MQ=255)
                               gagaTACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACt                                                                 >  1:3644039/1‑101 (MQ=255)
                               gagaTACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACt                                                                 >  2:3644041/1‑101 (MQ=255)
                                agaTACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTg                                                                >  2:3644044/1‑101 (MQ=255)
                                agaTACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTg                                                                >  2:3644047/1‑101 (MQ=255)
                                 gaTACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGGAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGc                                                               <  1:3643421/101‑1 (MQ=255)
                                  aTACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCt                                                              <  2:3643596/101‑1 (MQ=255)
                                  aTACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCt                                                              >  2:3644051/1‑101 (MQ=255)
                                   tACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCtt                                                             >  1:3644055/1‑101 (MQ=255)
                                   tACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGGGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCtt                                                             >  1:3644053/1‑101 (MQ=255)
                                      tttttGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTAtt                                                          >  2:3644059/1‑101 (MQ=255)
                                       ttttGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTg                                                         >  1:3644061/1‑101 (MQ=255)
                                       ttttGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTg                                                         >  2:3644062/1‑101 (MQ=255)
                                          tGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCtt                                                      >  1:3644070/1‑101 (MQ=255)
                                          tGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCtt                                                      >  2:3644067/1‑101 (MQ=255)
                                           gAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCttt                                                     >  2:3644071/1‑101 (MQ=255)
                                                gCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGAttt                                                >  2:3644075/1‑101 (MQ=255)
                                                 cTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGAtttt                                               <  1:3643513/101‑1 (MQ=255)
                                                 cTTTGCTGGTGATTAAAAATTAAGGAGGGGGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGAtttt                                               >  1:3644077/1‑101 (MQ=255)
                                                  tttGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTa                                              >  2:3644081/1‑101 (MQ=255)
                                                  ttagCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTa                                              >  1:3644084/4‑101 (MQ=255)
                                                   ttGCTTGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGAtttt                                               <  2:3643337/99‑1 (MQ=255)
                                                   ttGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGAtt                                                 <  1:3643545/97‑1 (MQ=255)
                                                   ttGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTAt                                             >  2:3644086/1‑101 (MQ=255)
                                                    tGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTAtt                                            >  2:3644089/1‑101 (MQ=255)
                                                    tGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTAtt                                            <  1:3643379/101‑1 (MQ=255)
                                                     gCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTc                                           >  1:3644090/1‑101 (MQ=255)
                                                        ggTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTa                                        >  2:3644095/1‑101 (MQ=255)
                                                         gTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTAt                                       <  2:3643549/101‑1 (MQ=255)
                                                         gTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTAt                                       <  1:3643482/101‑1 (MQ=255)
                                                          tGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATc                                      <  1:3643993/101‑1 (MQ=255)
                                                           gATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCt                                     >  2:3644099/1‑101 (MQ=255)
                                                           gATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCt                                     >  1:3644097/1‑101 (MQ=255)
                                                           gATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCt                                     <  1:3643816/101‑1 (MQ=255)
                                                             ttAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCttt                                   >  1:3644103/1‑101 (MQ=255)
                                                             ttAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCttt                                   <  1:3643671/101‑1 (MQ=255)
                                                               aaaaaTTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTcc                                 <  2:3643579/101‑1 (MQ=255)
                                                                aaaaTTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTccc                                >  2:3644107/1‑101 (MQ=255)
                                                                aaaaTTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTccc                                >  2:3644109/1‑101 (MQ=255)
                                                                   aTTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATg                             >  1:3644112/1‑101 (MQ=255)
                                                                    ttAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATgg                            <  1:3643638/101‑1 (MQ=255)
                                                                    ttAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATgg                            >  1:3644117/1‑101 (MQ=255)
                                                                    ttAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATgg                            >  1:3644115/1‑101 (MQ=255)
                                                                     tAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGt                           <  1:3643697/101‑1 (MQ=255)
                                                                      aaGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTa                          >  2:3644121/1‑101 (MQ=255)
                                                                      aaGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTa                          >  1:3644120/1‑101 (MQ=255)
                                                                      aaGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTa                          <  1:3643559/101‑1 (MQ=255)
                                                                         gaggGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTAccc                       >  2:3644126/1‑101 (MQ=255)
                                                                         gaggGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTAccc                       <  1:3643752/101‑1 (MQ=255)
                                                                          aggGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCg                      >  2:3644130/1‑101 (MQ=255)
                                                                          aggGTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCg                      >  1:3644129/1‑101 (MQ=255)
                                                                            ggTGTAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGa                    >  1:3644133/1‑101 (MQ=255)
                                                                             gtgtAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAg                   >  1:3644136/1‑101 (MQ=255)
                                                                             gtgtAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAg                   <  1:3643948/101‑1 (MQ=255)
                                                                               gtAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCg                 >  2:3644139/1‑101 (MQ=255)
                                                                                tAACGACAGGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCgg                >  1:3644141/1‑101 (MQ=255)
                                                                                tAACGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCgg                <  1:3643247/101‑1 (MQ=255)
                                                                                  aCGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGa              >  1:3644145/1‑101 (MQ=255)
                                                                                   cGACAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGAc             >  2:3644146/1‑101 (MQ=255)
                                                                                    gACAAGTTGCAGGCACAAAAAACCACCCGAAGGGGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCtt                                    <  2:3643481/77‑1 (MQ=255)
                                                                                     aCAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGACtt           >  1:3644150/1‑101 (MQ=255)
                                                                                      cAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGa                    >  1:3644155/1‑91 (MQ=255)
                                                                                      cAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGACTTg          <  2:3643520/101‑1 (MQ=255)
                                                                                      cAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGACTTg          >  2:3644153/1‑101 (MQ=255)
                                                                                      cAAGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGACTTg          <  2:3643471/101‑1 (MQ=255)
                                                                                       aaGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGACTTGa         >  2:3644160/1‑101 (MQ=255)
                                                                                       aaGTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGACTTGa         >  1:3644157/1‑101 (MQ=255)
                                                                                        aGTTGCAGGCACAAAAAACCACCCGAAGGGGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGACTTGaa        >  2:3644165/1‑101 (MQ=255)
                                                                                         gTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGACTTGAAc       <  2:3643647/101‑1 (MQ=255)
                                                                                         gTTGCAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGACTTGAAc       >  1:3644166/1‑101 (MQ=255)
                                                                                             cAGGCACAAAAACCCCCCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCt                                          <  1:3643701/62‑1 (MQ=255)
                                                                                             cAGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGACTTGAACCCGc   >  1:3644170/1‑101 (MQ=255)
                                                                                              aGGCACAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGACTTGAACCCGca  >  1:3644173/1‑101 (MQ=255)
                                                                                                  caaaaaaaaCCACCCGAAGGTGGTTTCACGACACTGCTTATTg                                                         <  1:3643766/41‑1 (MQ=255)
                                                                                                   aaaaaaacCCCCCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGAt                                                  <  1:3643623/48‑1 (MQ=255)
                                                                                                    |                                                                                               
CAGCACCTATCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGTACCCGGAGCGGGACTTGAACCCGCA  >  NC_000913/1991648‑1991843

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