Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,582,212 T→C intergenic (+22/‑178) b4547 → / → yrhB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,582,2120TC85.4% 321.1 / 44.4 123intergenic (+22/‑178)b4547/yrhBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/15);  new base C (17/88);  total (20/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

Only 100 of 127 total aligned reads displayed.

CTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAACCAGCATTTAATAATA  >  NC_000913/3582121‑3582304
                                                                                           |                                                                                            
cTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGtt                                                                                             <  1:2571910/93‑1 (MQ=255)
 tATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTc                                                                                            <  2:5482261/93‑1 (MQ=255)
  atgatgATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGctct                                                                                           <  2:6572445/93‑1 (MQ=255)
  atgatgATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCt                                                                                           <  1:850016/93‑1 (MQ=255)
   tgatgaTTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTg                                                                                          >  1:1256792/1‑93 (MQ=255)
   tgatgaTTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTg                                                                                          <  2:3493916/93‑1 (MQ=255)
     atgatTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGaa                                                                                        >  2:201630/1‑93 (MQ=255)
       gatTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGaaaa                                                                                      <  1:2086303/93‑1 (MQ=255)
        atTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGaaaaa                                                                                     <  2:2915467/93‑1 (MQ=255)
        atTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGaaaaa                                                                                     >  1:1904324/1‑93 (MQ=255)
        atTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGaaaaa                                                                                     <  2:4267981/93‑1 (MQ=255)
        atTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGaaaaa                                                                                     <  2:1850285/93‑1 (MQ=255)
        atTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGaaaaa                                                                                     <  2:1699115/93‑1 (MQ=255)
         ttttAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAg                                                                                    >  2:1193540/1‑93 (MQ=255)
         ttttAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAg                                                                                    <  2:5846068/93‑1 (MQ=255)
         ttttAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAg                                                                                    <  1:5906970/93‑1 (MQ=255)
            tAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGaaa                                                                                 <  1:5748224/93‑1 (MQ=255)
            tAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGaaa                                                                                 >  2:6609099/1‑93 (MQ=255)
            tAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGaaa                                                                                 <  2:3217137/93‑1 (MQ=255)
            tAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGaaa                                                                                 <  1:3628745/93‑1 (MQ=255)
             aaTGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAAt                                                                                >  2:2702069/1‑93 (MQ=255)
              aTGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATc                                                                               <  1:6089266/93‑1 (MQ=255)
               tGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCa                                                                              <  1:1464305/93‑1 (MQ=255)
               tGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCa                                                                              <  1:2158212/93‑1 (MQ=255)
                gACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCAt                                                                             <  2:3857338/93‑1 (MQ=255)
                gACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCAt                                                                             <  1:5857355/93‑1 (MQ=255)
                gACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCAt                                                                             <  2:3113842/93‑1 (MQ=255)
                gACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAACCAt                                                                             <  2:649843/93‑1 (MQ=255)
                  cATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATaa                                                                           <  2:2472989/93‑1 (MQ=255)
                  cATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATaa                                                                           <  2:3057479/93‑1 (MQ=255)
                       tGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAatat                                                                      >  2:6218916/1‑93 (MQ=255)
                        gATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATAtt                                                                     >  2:1747469/1‑93 (MQ=255)
                         atttatttAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTa                                                                    >  1:3840570/1‑93 (MQ=255)
                         atttatttAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTa                                                                    <  2:2872651/93‑1 (MQ=255)
                          tttatttaGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAg                                                                   <  1:3393340/93‑1 (MQ=255)
                              tttaGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGa                                                               <  1:6243219/93‑1 (MQ=255)
                              tttaGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGa                                                               <  1:1070705/93‑1 (MQ=255)
                              tttaGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGa                                                               <  1:4806010/93‑1 (MQ=255)
                              tttaGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGa                                                               <  1:436285/93‑1 (MQ=255)
                                taGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATa                                                             >  2:6541934/1‑93 (MQ=255)
                                taGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATa                                                             <  1:3279947/93‑1 (MQ=255)
                                taGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATa                                                             <  2:1367934/93‑1 (MQ=255)
                                 aGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATGTTAGCAGATAc                                                            <  2:4028004/93‑1 (MQ=255)
                                 aGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATAc                                                            <  1:6012873/93‑1 (MQ=255)
                                 aGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATAc                                                            <  2:3707014/93‑1 (MQ=255)
                                  gAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACt                                                           <  1:856308/93‑1 (MQ=255)
                                   aaTTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACtt                                                          <  2:3923050/93‑1 (MQ=255)
                                    aTTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTa                                                         <  2:408280/93‑1 (MQ=255)
                                       ttACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTAttt                                                      >  2:876176/1‑93 (MQ=255)
                                        tACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCATTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTc                                                     <  1:6555304/93‑1 (MQ=255)
                                         aCCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCg                                                    <  2:4057595/93‑1 (MQ=255)
                                         aCCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCg                                                    <  2:4301001/93‑1 (MQ=255)
                                         aCCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCg                                                    >  2:3454393/1‑93 (MQ=255)
                                          cccGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGt                                                   <  2:2531796/93‑1 (MQ=255)
                                          cccGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGt                                                   <  1:2783818/93‑1 (MQ=255)
                                          cccGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGt                                                   <  1:5271830/93‑1 (MQ=255)
                                          cccGAATCATAAGATTCTGATAACCTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGt                                                   <  2:2169440/93‑1 (MQ=255)
                                            cGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGa                                                 <  2:2980903/93‑1 (MQ=255)
                                              aaTCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACt                                               >  2:3845096/1‑93 (MQ=255)
                                              aaTCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACt                                               >  2:3845095/1‑93 (MQ=255)
                                               aTCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACta                                              <  2:329125/93‑1 (MQ=255)
                                               aTCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACta                                              <  1:837137/93‑1 (MQ=255)
                                                tCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACtat                                             <  2:2198802/93‑1 (MQ=255)
                                                tCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACtat                                             <  2:4843637/93‑1 (MQ=255)
                                                  aTAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTatat                                           <  1:3944847/93‑1 (MQ=255)
                                                  aTAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTatat                                           <  2:1402514/93‑1 (MQ=255)
                                                   tAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATc                                          <  1:6156483/93‑1 (MQ=255)
                                                       aTTCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTAt                                      <  1:349114/93‑1 (MQ=255)
                                                        ttCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATc                                     <  2:906490/93‑1 (MQ=255)
                                                         tCTGATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATcc                                    <  1:428596/93‑1 (MQ=255)
                                                            gATATCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGAt                                 <  1:3589480/93‑1 (MQ=255)
                                                             atatCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATg                                <  1:3589722/93‑1 (MQ=255)
                                                                tCTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAAc                             <  1:5392877/93‑1 (MQ=255)
                                                                 cTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACt                            <  1:1215265/93‑1 (MQ=255)
                                                                 cTTAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACt                            <  1:1733796/93‑1 (MQ=255)
                                                                   tAAATAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACttt                          <  1:1416118/93‑1 (MQ=255)
                                                                     aaTAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACttttt                        <  2:5543640/93‑1 (MQ=255)
                                                                      aTAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACtttttt                       <  2:854373/93‑1 (MQ=255)
                                                                      aTAAGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACtttttt                       <  1:1427387/93‑1 (MQ=255)
                                                                        aaGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCa                     <  1:1302512/93‑1 (MQ=255)
                                                                         aGAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCaa                    <  1:2236563/93‑1 (MQ=255)
                                                                          gAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCaaa                   <  2:3575472/93‑1 (MQ=255)
                                                                          gAAAATCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCaaa                   <  1:2709089/93‑1 (MQ=255)
                                                                          aaaaaCCCCGTTCGACACTCTGAAAAAGAAATCATAAAAGAGTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCaaa                   <  1:6619176/92‑1 (MQ=255)
                                                                           aaaaTCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCaaaa                  <  2:3027686/93‑1 (MQ=255)
                                                                           aaaaTCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCaaaa                  <  2:4503153/93‑1 (MQ=255)
                                                                            aaaTCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAAc                 <  1:6218916/93‑1 (MQ=255)
                                                                             aaTCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAAcc                <  1:2324407/93‑1 (MQ=255)
                                                                             aaTCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAAcc                <  2:2153469/93‑1 (MQ=255)
                                                                             aaTCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAAcc                <  1:3180286/93‑1 (MQ=255)
                                                                              aTCCCGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAACCa               <  1:2674937/93‑1 (MQ=255)
                                                                                 ccGTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAACCAGCa            <  2:5109234/93‑1 (MQ=255)
                                                                                   gTTCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAACCAGCAtt          <  2:4904898/93‑1 (MQ=255)
                                                                                    ttCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAACCAGCAttt         <  2:4218270/93‑1 (MQ=255)
                                                                                     tCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAACCAGCATTta        <  1:5385692/93‑1 (MQ=255)
                                                                                     tCGACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAACCAGCATTta        <  1:5501723/93‑1 (MQ=255)
                                                                                      cGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGATCTTTTTTCAAAACCAGCATTtaa       <  2:2740057/93‑1 (MQ=255)
                                                                                       gACGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAACCAGCATTtaat      >  1:694660/1‑93 (MQ=255)
                                                                                        aCGCTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAACCAGCATTtaata     <  1:2061404/93‑1 (MQ=255)
                                                                                           ttCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAACCAGCATTTAataata  <  1:236374/93‑1 (MQ=255)
                                                                                           |                                                                                            
CTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATCCGATGAACTTTTTTCAAAACCAGCATTTAATAATA  >  NC_000913/3582121‑3582304

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