Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,431,203 T→C intergenic (‑6/‑124) smf ← / → def putative Rossmann‑fold nucleotide‑binding protein/peptide deformylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,431,2030TC99.5% 682.7 / NA 214intergenic (‑6/‑124)smf/defputative Rossmann‑fold nucleotide‑binding protein/peptide deformylase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (116/97);  total (116/98)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

ACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATA  >  NC_000913/3431115‑3431291
                                                                                        |                                                                                        
aCCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTcc                                                                                      >  1:4372639/1‑93 (MQ=255)
 cccAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTccc                                                                                     <  1:3699858/93‑1 (MQ=255)
  ccAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCt                                                                                    >  2:3389316/1‑93 (MQ=255)
   cAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTg                                                                                   <  1:2002318/93‑1 (MQ=255)
      tGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCa                                                                                <  2:866982/93‑1 (MQ=255)
         gcaaTACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATAAAACGCAGCAAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATaa                                                                             <  2:3257160/90‑1 (MQ=255)
         gCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATaa                                                                             >  1:3145189/1‑93 (MQ=255)
           gATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAagc                                                                           <  2:3064438/93‑1 (MQ=255)
           gATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAagc                                                                           <  1:2847824/93‑1 (MQ=255)
           gATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAagc                                                                           >  2:4448868/1‑93 (MQ=255)
            aTACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAagca                                                                          >  2:638890/1‑93 (MQ=255)
             tACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAagcag                                                                         >  2:809966/1‑93 (MQ=255)
             tACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAagcag                                                                         >  1:937996/1‑93 (MQ=255)
             tACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAagcag                                                                         <  2:4324459/93‑1 (MQ=255)
              aCGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAagcagc                                                                        <  2:1538233/93‑1 (MQ=255)
               cGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGcc                                                                       <  1:1545645/93‑1 (MQ=255)
                  aCCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTa                                                                    <  2:207164/93‑1 (MQ=255)
                   ccATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAg                                                                   >  2:706439/1‑93 (MQ=255)
                      tatCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCaa                                                                <  1:485808/93‑1 (MQ=255)
                        tcatcGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATc                                                              >  2:2520595/1‑93 (MQ=255)
                        tcatcGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATc                                                              >  2:4098630/1‑93 (MQ=255)
                        tcatcGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATc                                                              >  2:1130317/1‑93 (MQ=255)
                          atcGCCGTACAAGCTTCTGATTCGCATTAACCGCAGCCAAATTTCTGTATCGCCCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCtt                                                            <  2:3186235/93‑1 (MQ=255)
                          atcGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCtt                                                            <  1:3752668/93‑1 (MQ=255)
                            cGCCGTACAAGCTGCTGATACTCATTAACCGCGGCCAAATTTCTGTATGGACCATCCTTCCCTCCCTGCCATAAGCAGCCTTAGCAATCTTTg                                                          <  1:4037205/93‑1 (MQ=255)
                            cGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTg                                                          >  1:3997324/1‑93 (MQ=255)
                               cGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGGCTTAGCAATCTTTGCGa                                                       >  1:2263757/1‑93 (MQ=255)
                                gTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGAt                                                      >  2:2462298/1‑93 (MQ=255)
                                gTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGAt                                                      >  1:3784294/1‑93 (MQ=255)
                                 tACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTCTCGACCATCCTTATCCCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGAtt                                                     <  2:3513882/93‑1 (MQ=255)
                                 tACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGAtt                                                     >  2:163269/1‑93 (MQ=255)
                                 tACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGAtt                                                     <  1:2398091/93‑1 (MQ=255)
                                    aaGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGt                                                  <  1:875830/93‑1 (MQ=255)
                                    aaGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGt                                                  <  2:2147321/93‑1 (MQ=255)
                                     aGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTc                                                 >  1:3600024/1‑93 (MQ=255)
                                     aGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTc                                                 >  1:2460234/1‑93 (MQ=255)
                                      gctgctGATACTCATTAAACCCAGCGAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCa                                                <  2:1516219/93‑1 (MQ=255)
                                       ctgctaATACTCATTAAATGAAGCGAAATTTCTGTATCAACCTTCCTTCCCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAg                                               <  2:2248198/93‑1 (MQ=255)
                                        tgctgATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGt                                              >  1:2971723/1‑93 (MQ=255)
                                          ctgATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGa                                            <  1:2497487/93‑1 (MQ=255)
                                           tgATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGAt                                           >  1:1752687/1‑93 (MQ=255)
                                             ttACTCATTAACCGCAGCCAATTTTCTGTATCGCCCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGc                                         <  2:2318322/92‑1 (MQ=255)
                                             aTACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGc                                         <  2:3576719/93‑1 (MQ=255)
                                             aTACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGc                                         <  1:156344/93‑1 (MQ=255)
                                             aTACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGc                                         <  1:3196179/93‑1 (MQ=255)
                                              tACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCt                                        >  2:2816546/1‑93 (MQ=255)
                                              tACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCt                                        >  1:575198/1‑93 (MQ=255)
                                              tACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCt                                        >  1:3353358/1‑93 (MQ=255)
                                               aCTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTg                                       >  2:2006713/1‑93 (MQ=255)
                                                cTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGt                                      >  1:2779345/1‑93 (MQ=255)
                                                cTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGt                                      <  2:2397426/93‑1 (MQ=255)
                                                cTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGt                                      >  2:4269208/1‑93 (MQ=255)
                                                cTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGt                                      >  2:1481798/1‑93 (MQ=255)
                                                 tCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCACTCTTTGCGATTGGTCAGTGATGCTGTc                                     >  2:1187717/1‑93 (MQ=255)
                                                 tCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTc                                     >  1:1392201/1‑93 (MQ=255)
                                                 tCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTc                                     <  1:1002893/93‑1 (MQ=255)
                                                 tCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTc                                     >  2:1027010/1‑93 (MQ=255)
                                                 tCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTc                                     >  1:4447629/1‑93 (MQ=255)
                                                  cATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTCGCAATCTTTGCGATTGGTCAGTGATGCTGTCa                                    >  2:1710153/1‑93 (MQ=255)
                                                   aTTAAACGCAGCCCAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCaa                                   >  2:4042042/1‑93 (MQ=255)
                                                   aTTAAACGCAGCCAAATTTCTGTATCGACCATCCTTCCCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCaa                                   <  1:1849413/93‑1 (MQ=255)
                                                   aTTAAACGCAGCCAAATTTCTGTATCGACCATCCTTCCCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCaa                                   <  1:3731471/93‑1 (MQ=255)
                                                   aTTAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCaa                                   >  2:2754173/1‑93 (MQ=255)
                                                    ttAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAAt                                  >  1:4145309/1‑93 (MQ=255)
                                                    ttAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAAt                                  >  1:1255023/1‑93 (MQ=255)
                                                    ttAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAAt                                  >  2:4211297/1‑93 (MQ=255)
                                                     tAAACGCAGCCAAATTTCTGTATCGACCATCCTTCCCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATc                                 <  1:3177125/93‑1 (MQ=255)
                                                     tAAACGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATc                                 >  1:1812908/1‑93 (MQ=255)
                                                      aaaCGCAGCCAAATTTCTGTATCGACCATCCTTCCCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCa                                <  1:3244398/93‑1 (MQ=255)
                                                        aCGCAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCaga                              <  2:165249/93‑1 (MQ=255)
                                                         cGCACCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCagag                             <  1:3273935/93‑1 (MQ=255)
                                                           cAGCCAAATTTCTGTATCGACTATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATTTTTGCGATTGGTCAGTGATGCTGTCAATCAGAggg                           <  1:2149404/93‑1 (MQ=255)
                                                           cAGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAggg                           >  1:122630/1‑93 (MQ=255)
                                                           cAGCCAAATTTCTGTATCGACCATCCTGCCCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAggg                           <  2:364634/93‑1 (MQ=255)
                                                            aGCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAgggg                          >  1:3667416/1‑93 (MQ=255)
                                                             gCCAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAggggg                         <  2:2209809/93‑1 (MQ=255)
                                                               cAAATTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGAt                       >  1:756112/1‑93 (MQ=255)
                                                                aaaTTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGAtt                      >  1:1072774/1‑93 (MQ=255)
                                                                aaaTTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGAtt                      >  1:301059/1‑93 (MQ=255)
                                                                 aaTTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGAttt                     >  1:1147725/1‑93 (MQ=255)
                                                                 aaTTTCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGAttt                     >  1:2597276/1‑93 (MQ=255)
                                                                   tttCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGt                   >  1:211700/1‑93 (MQ=255)
                                                                    ttCTGTATCGCCCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTc                  <  2:4354194/93‑1 (MQ=255)
                                                                    ttCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTc                  <  1:461310/93‑1 (MQ=255)
                                                                    ttCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTc                  >  2:3841443/1‑93 (MQ=255)
                                                                     tCTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCt                 >  2:273441/1‑93 (MQ=255)
                                                                      ttGTATCGACCATCCTTCCCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCta                <  1:2905091/92‑1 (MQ=255)
                                                                      cTGTATCGACCTTCCTTACCACCCTGCCATAAGCATCATTAGCAATCTTTCCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCta                <  2:767003/93‑1 (MQ=255)
                                                                      cTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCta                <  2:994878/93‑1 (MQ=255)
                                                                      cTGTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCta                >  2:4482972/1‑93 (MQ=255)
                                                                        gTATCGACCATCCTTCCCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCtaga              <  1:1102226/93‑1 (MQ=255)
                                                                        gTATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCtaga              <  1:4312414/93‑1 (MQ=255)
                                                                         tATCGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCtagaa             <  1:2372623/93‑1 (MQ=255)
                                                                           tCGACCACCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCtagaata           <  2:620729/93‑1 (MQ=255)
                                                                            cGACCATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCtagaatag          <  2:781116/93‑1 (MQ=255)
                                                                                cATCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATagaaga      >  2:2497852/1‑93 (MQ=255)
                                                                                  tCCTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGaaa    <  2:4217699/93‑1 (MQ=255)
                                                                                    cTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAaata  >  1:3913934/1‑93 (MQ=255)
                                                                                    cTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAaata  >  1:4079631/1‑93 (MQ=255)
                                                                                    cTTACCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAaata  >  2:1847133/1‑93 (MQ=255)
                                                                                        |                                                                                        
ACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATA  >  NC_000913/3431115‑3431291

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