Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 461,311 (T)7→8 intergenic (+69/‑140) lon → / → hupB Lon protease/DNA‑binding protein HU‑beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913461,3041.T99.6% 1096.9 / ‑3.7 252intergenic (+62/‑147)lon/hupBLon protease/DNA‑binding protein HU‑beta
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (129/122);  total (129/123)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

AATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGATATAACTGCTGCGCGT  >  NC_000913/461204‑461398
                                                                                                     |                                                                                              
 aTGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCt                                                                                               <  2:787433/101‑1 (MQ=255)
 aTGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCt                                                                                               >  2:787829/1‑101 (MQ=255)
   gAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttt                                                                                             >  1:787832/1‑101 (MQ=255)
     aCCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttt                                                                                           <  2:787489/101‑1 (MQ=255)
       cGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttttt                                                                                         >  2:787835/1‑101 (MQ=255)
        gTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCtttttttt                                                                                        >  1:787837/1‑100 (MQ=255)
         tCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttttttg                                                                                       <  2:787396/101‑3 (MQ=255)
         tCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttttttg                                                                                       >  1:787840/1‑99 (MQ=255)
          cTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttttttgt                                                                                      <  2:787333/101‑4 (MQ=255)
          cTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttttttgt                                                                                      >  2:787845/1‑98 (MQ=255)
          cTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttttttgt                                                                                      >  2:787842/1‑98 (MQ=255)
            ggTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttttttgtct                                                                                    <  2:787324/101‑6 (MQ=255)
            ggTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttttttgtct                                                                                    >  1:787848/1‑96 (MQ=255)
            ggTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttttttgtct                                                                                    <  1:787393/101‑6 (MQ=255)
               aTGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGtctc                                                                                    >  2:787853/1‑99 (MQ=255)
               aTGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGc                                                                                  >  1:787851/1‑101 (MQ=255)
                tGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCt                                                                                 <  1:787677/101‑1 (MQ=255)
                 gCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTa                                                                                <  2:787257/101‑1 (MQ=255)
                 gCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTa                                                                                <  2:787304/101‑1 (MQ=255)
                  cAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTaa                                                                               >  2:787856/1‑101 (MQ=255)
                  cAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTaa                                                                               <  2:787270/101‑1 (MQ=255)
                   aGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAg                                                                              <  2:787610/101‑1 (MQ=255)
                     gTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGtt                                                                            >  1:787858/1‑101 (MQ=255)
                        gtgACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGa                                                                         >  1:787862/1‑101 (MQ=255)
                         tgACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGAt                                                                        >  2:787864/1‑101 (MQ=255)
                         tgACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGAt                                                                        <  2:787264/101‑1 (MQ=255)
                           aCTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATgg                                                                      >  2:787867/1‑101 (MQ=255)
                             tGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCg                                                                    >  1:787871/1‑101 (MQ=255)
                             tGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCg                                                                    <  2:787445/101‑1 (MQ=255)
                              gCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCgg                                                                   <  2:787348/101‑1 (MQ=255)
                              gCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCgg                                                                   >  2:787873/1‑101 (MQ=255)
                                aaaaTAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGAt                                                                 >  1:787876/1‑101 (MQ=255)
                                 aaaTAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATc                                                                <  2:787474/101‑1 (MQ=255)
                                  aaTAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCg                                                               <  1:787028/101‑1 (MQ=255)
                                  aaTAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCg                                                               <  1:787467/101‑1 (MQ=255)
                                  aaTAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCg                                                               <  1:787359/101‑1 (MQ=255)
                                   aTAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCgg                                                              <  1:787381/101‑1 (MQ=255)
                                   aTAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCgg                                                              <  1:787077/101‑1 (MQ=255)
                                     aGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGc                                                            >  1:787883/1‑101 (MQ=255)
                                     aGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGc                                                            >  2:787881/1‑101 (MQ=255)
                                     aGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGc                                                            >  2:787884/1‑101 (MQ=255)
                                      gTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCt                                                           >  2:787887/1‑101 (MQ=255)
                                        gACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGc                                                                     >  1:787895/1‑89 (MQ=255)
                                        gACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTg                                                         >  1:787893/1‑101 (MQ=255)
                                        gACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTg                                                         >  2:787892/1‑101 (MQ=255)
                                          ccTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCgg                                                                   <  2:787154/89‑1 (MQ=255)
                                           cTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGccc                                                      >  1:787899/1‑101 (MQ=255)
                                            tCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCt                                                     <  1:787403/101‑1 (MQ=255)
                                             cgcgCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCtt                                                    <  2:787048/101‑1 (MQ=255)
                                               cgcAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGa                                                                  >  2:787901/1‑85 (MQ=255)
                                               cgcAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCttat                                                  <  1:787793/101‑1 (MQ=255)
                                                gcAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCttatt                                                 <  1:787365/101‑1 (MQ=255)
                                                gcAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCttatt                                                 <  1:787565/101‑1 (MQ=255)
                                                  aaaaTGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTaa                                               >  2:787905/1‑101 (MQ=255)
                                                   aaaTGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAg                                              >  2:787909/1‑101 (MQ=255)
                                                   aaaTGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAg                                              >  1:787910/1‑101 (MQ=255)
                                                    aaTGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttttttg                                                                                       >  2:787914/1‑56 (MQ=255)
                                                     aTGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAggg                                            <  1:787282/101‑1 (MQ=255)
                                                      tGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAgggg                                           >  1:787917/1‑101 (MQ=255)
                                                       gCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGgt                                          <  2:787224/101‑1 (MQ=255)
                                                        cACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAgggg                                           >  2:787921/1‑99 (MQ=255)
                                                         aCTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCg                                                               >  2:787925/1‑78 (MQ=255)
                                                         aCTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGgtgt                                        >  1:787923/1‑101 (MQ=255)
                                                           taataaAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGtgttg                                      >  1:787929/1‑101 (MQ=255)
                                                            aataaAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAgggg                                           >  2:787933/1‑95 (MQ=255)
                                                            aataaAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGtgttgt                                     <  1:787223/101‑1 (MQ=255)
                                                               aaaaaCAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAg                                  <  2:787462/101‑1 (MQ=255)
                                                                aaaaCAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAgg                                 >  2:787936/1‑101 (MQ=255)
                                                                 aaaCAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAggg                                >  1:787939/1‑101 (MQ=255)
                                                                  aaCAGGGTTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAgggg                               <  2:787290/101‑1 (MQ=255)
                                                                   aCAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAggg                                >  1:787944/1‑99 (MQ=255)
                                                                   aCAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGa                              >  2:787941/1‑101 (MQ=255)
                                                                    cAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGAt                             <  2:787526/101‑1 (MQ=255)
                                                                     aGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGAtg                            <  1:787446/101‑1 (MQ=255)
                                                                      gggCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGgtg                                         >  2:787948/1‑87 (MQ=255)
                                                                      gggCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGAtgg                           <  1:787498/101‑1 (MQ=255)
                                                                        gCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGggttgt                                     >  2:787951/1‑89 (MQ=255)
                                                                         cTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGAtggctg                        <  1:787544/101‑1 (MQ=255)
                                                                         cTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGAtggctg                        >  1:787952/1‑101 (MQ=255)
                                                                          tGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGAtggctgg                       <  1:787299/101‑1 (MQ=255)
                                                                           ggCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGAtggctggc                      <  1:787243/101‑1 (MQ=255)
                                                                            gCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGcc                     <  1:787465/101‑1 (MQ=255)
                                                                              aGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGGCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTg                   >  1:787961/1‑101 (MQ=255)
                                                                              aGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTg                   <  1:787413/101‑1 (MQ=255)
                                                                              aGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTg                   >  2:787957/1‑101 (MQ=255)
                                                                              aGGCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTg                   >  1:787958/1‑101 (MQ=255)
                                                                                gCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGat                 <  2:787405/101‑1 (MQ=255)
                                                                                gCTAATTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGat                 >  1:787964/1‑101 (MQ=255)
                                                                                    aTTCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGATATaa             >  1:787968/1‑101 (MQ=255)
                                                                                     ttCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCtt                                                          <  2:787364/55‑1 (MQ=255)
                                                                                      tCGGGCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGATATAAct           <  1:787205/101‑1 (MQ=255)
                                                                                        gggCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGATATAActgc         <  2:787420/101‑1 (MQ=255)
                                                                                        gggCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGATATAActgc         >  2:787970/1‑101 (MQ=255)
                                                                                         ggCTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGATATAActgct        <  1:787186/101‑1 (MQ=255)
                                                                                           cTTGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGATATAACtgctgc      >  1:787975/1‑101 (MQ=255)
                                                                                             tGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGATATAACTGCTgcgc    <  1:787330/101‑1 (MQ=255)
                                                                                             tGCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGATATAACTGCTgcgc    <  2:787571/101‑1 (MQ=255)
                                                                                              gCCAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGATATAACTGCTGcgcg   <  2:787162/101‑1 (MQ=255)
                                                                                               ccAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGAtg                            >  1:787978/1‑75 (MQ=255)
                                                                                               ccAGCCTTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGATATAACTGCTGCGCGt  <  2:787659/101‑1 (MQ=255)
                                                                                                     |                                                                                              
AATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGCTGGCCTGATATAACTGCTGCGCGT  >  NC_000913/461204‑461398

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