Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,021,493 C→T A121A (GCC→GCT) 
A117A (GCC→GCT) 		 tatD →
b4556 →		 Mg‑dependent DNase
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,021,4930CT99.6% 987.7 / NA 262A121A (GCC→GCT) 
A117A (GCC→GCT) 		tatD
b4556		Mg‑dependent DNase
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (129/132);  minor base G (0/1);  total (129/133)
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 271 total aligned reads displayed.

GGCGATTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCCGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAACTGCC  >  NC_000913/4021403‑4021582
                                                                                          |                                                                                         
ggCGATTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATtgctgc                                                                                         <  1:6413809/93‑1 (MQ=255)
  cGATTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAg                                                                                       <  1:1566948/93‑1 (MQ=255)
  cGATTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAg                                                                                       <  1:2536449/93‑1 (MQ=255)
   gATTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGa                                                                                      >  2:833797/1‑93 (MQ=255)
   gATTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGa                                                                                      <  2:6678581/93‑1 (MQ=255)
    aTTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCGGCTACGCATTGCTGCAGAt                                                                                     <  2:3497950/93‑1 (MQ=255)
    aTTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGAt                                                                                     <  1:796373/93‑1 (MQ=255)
    aTTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGAt                                                                                     <  1:87394/93‑1 (MQ=255)
      tGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGAttt                                                                                   >  2:2848491/1‑93 (MQ=255)
         tGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTaaa                                                                                <  1:5206834/93‑1 (MQ=255)
         tGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTaaa                                                                                <  2:5831313/93‑1 (MQ=255)
           aaTGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACa                                                                              <  1:6651537/93‑1 (MQ=255)
              gtgGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGc                                                                           >  1:5753903/1‑93 (MQ=255)
                ggTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCg                                                                         <  1:5097408/93‑1 (MQ=255)
                 gTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCgg                                                                        <  1:4746709/93‑1 (MQ=255)
                 gTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCgg                                                                        <  1:4131949/93‑1 (MQ=255)
                 gTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCgg                                                                        >  2:1573287/1‑93 (MQ=255)
                  tctcGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGt                                                                       <  1:2416126/93‑1 (MQ=255)
                  tctcGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGt                                                                       <  1:2710416/93‑1 (MQ=255)
                   ctcGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTa                                                                      <  1:3946741/93‑1 (MQ=255)
                    tcttcTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTAt                                                                     <  2:1044949/89‑1 (MQ=255)
                    tcGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTAt                                                                     >  2:2039979/1‑93 (MQ=255)
                      gACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTAttt                                                                   >  1:3436918/1‑93 (MQ=255)
                       aCTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTa                                                                  <  1:3311087/93‑1 (MQ=255)
                        cTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTAt                                                                 >  1:3392502/1‑93 (MQ=255)
                        cTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTAt                                                                 <  2:5370588/93‑1 (MQ=255)
                        cTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTAt                                                                 <  2:2877709/93‑1 (MQ=255)
                         tttAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATg                                                                <  2:6580797/93‑1 (MQ=255)
                          ttAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGc                                                               <  1:4967925/93‑1 (MQ=255)
                            aaCCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCAc                                                             <  1:4902301/93‑1 (MQ=255)
                            aaCCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCAc                                                             <  1:5360030/93‑1 (MQ=255)
                            aaCCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCAc                                                             >  1:6240221/1‑93 (MQ=255)
                            aaCCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCAc                                                             >  1:214454/1‑93 (MQ=255)
                              ccGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTg                                                           <  2:4586930/93‑1 (MQ=255)
                              ccGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTg                                                           <  2:5131779/93‑1 (MQ=255)
                              ccGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGGATTTATGCACTg                                                           >  1:5232459/1‑93 (MQ=255)
                                gCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTc                                                         <  1:2341553/93‑1 (MQ=255)
                                 cAACTTTTTGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGtag                                                        >  2:924279/1‑91 (MQ=255)
                                 cAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTcg                                                        <  2:1456735/93‑1 (MQ=255)
                                  aaCTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTcgc                                                       >  2:4576120/1‑93 (MQ=255)
                                   aCTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTcgcg                                                      >  1:4255530/1‑93 (MQ=255)
                                    cTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCTCGa                                                     >  2:1605908/1‑93 (MQ=255)
                                     ttttCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGAt                                                    <  2:4911433/93‑1 (MQ=255)
                                     ttttCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGAt                                                    <  2:2811218/93‑1 (MQ=255)
                                      tttCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATg                                                   <  2:3587549/93‑1 (MQ=255)
                                       ttCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGc                                                  <  1:463662/93‑1 (MQ=255)
                                         cGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGccc                                                >  2:6121826/1‑93 (MQ=255)
                                         cGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGccc                                                <  2:4608148/93‑1 (MQ=255)
                                         cGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGccc                                                >  2:3128789/1‑93 (MQ=255)
                                         cGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGccc                                                <  1:6442528/93‑1 (MQ=255)
                                         cGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGccc                                                >  2:294972/1‑93 (MQ=255)
                                          gACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCa                                               <  1:4582835/93‑1 (MQ=255)
                                          gACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCa                                               <  1:299490/93‑1 (MQ=255)
                                           aCGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCAc                                              >  2:5800733/1‑93 (MQ=255)
                                           aCGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCAc                                              <  2:4417814/93‑1 (MQ=255)
                                           aCGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCAc                                              >  2:1278831/1‑93 (MQ=255)
                                            cGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACg                                             <  1:4038081/93‑1 (MQ=255)
                                             gCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGa                                            <  2:6382812/93‑1 (MQ=255)
                                                ggAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCg                                         <  2:4811131/93‑1 (MQ=255)
                                                ggAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCg                                         >  1:5275889/1‑93 (MQ=255)
                                                 gAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCgg                                        >  1:5941358/1‑93 (MQ=255)
                                                 gAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCgg                                        >  1:4545143/1‑93 (MQ=255)
                                                 gAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGGATTTATGCACTGTCGCGATGCCCACGAGCgg                                        >  2:5619877/1‑93 (MQ=255)
                                                   agagCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGtt                                      >  2:3412989/1‑93 (MQ=255)
                                                     agCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTa                                    <  2:594665/93‑1 (MQ=255)
                                                        aGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGa                                 >  1:1301190/1‑93 (MQ=255)
                                                         ggAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGAc                                >  1:6156177/1‑93 (MQ=255)
                                                         ggAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGAc                                >  2:3887938/1‑93 (MQ=255)
                                                         ggAACGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGCATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGAc                                >  1:6285947/1‑93 (MQ=255)
                                                           aaCGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACAt                              >  2:6498850/1‑93 (MQ=255)
                                                            aCGCGCTTTTGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACAtt                             >  2:2390079/1‑93 (MQ=255)
                                                               cgcTTTTGTTGCCCAGCAACGCATTGCTGCAGATTCAAACATGCCGGTATGGATGCACTGTCGCGATGCCCACGAGCGGTTTATGCCAGTGCt                          <  1:2573986/93‑1 (MQ=255)
                                                                  ttttGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGa                       >  2:1170229/1‑93 (MQ=255)
                                                                   tttGTTGCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAg                      >  2:3931835/1‑93 (MQ=255)
                                                                     tgttgCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGcc                    <  1:3839170/93‑1 (MQ=255)
                                                                       ttgCCCAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGt                  <  1:176386/93‑1 (MQ=255)
                                                                           ccAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCt              <  1:6552482/93‑1 (MQ=255)
                                                                            cAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTg             >  1:4200587/1‑93 (MQ=255)
                                                                            cAGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTg             <  1:964204/93‑1 (MQ=255)
                                                                             aGCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTgg            <  2:3673471/93‑1 (MQ=255)
                                                                              gCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGa           >  2:5343760/1‑93 (MQ=255)
                                                                              gCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGa           <  1:1938666/93‑1 (MQ=255)
                                                                              gCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGa           >  1:361043/1‑93 (MQ=255)
                                                                              gCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGa           >  1:2817806/1‑93 (MQ=255)
                                                                              gCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGa           <  1:2855947/93‑1 (MQ=255)
                                                                              gCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGa           >  2:3228857/1‑93 (MQ=255)
                                                                              gCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGa           >  2:3165771/1‑93 (MQ=255)
                                                                              gCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGa           >  2:1768082/1‑93 (MQ=255)
                                                                              gCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCCTGGCTGGa           >  1:1388989/1‑93 (MQ=255)
                                                                              gCTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGCTTATGACATTGCTGGAGCCGTGGCTGGa           >  2:4478771/1‑93 (MQ=255)
                                                                               cTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGAt          >  1:4312408/1‑93 (MQ=255)
                                                                               cTACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGAt          <  1:5335784/93‑1 (MQ=255)
                                                                                tACGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATa         >  1:1016668/1‑93 (MQ=255)
                                                                                  cGCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATaaa       <  1:1088078/93‑1 (MQ=255)
                                                                                   gCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAAc      <  1:6025017/93‑1 (MQ=255)
                                                                                   gCATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAAc      <  1:5453891/93‑1 (MQ=255)
                                                                                    cATTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAACt     >  1:6342883/1‑93 (MQ=255)
                                                                                     aTTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAACTg    >  2:5049553/1‑93 (MQ=255)
                                                                                     aTTGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAACTg    >  2:2188029/1‑93 (MQ=255)
                                                                                      ttGCTGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAACTGc   >  2:6327115/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGCGATTGGTGAATGTGGTCTCGACTTTAACCGCAACTTTTCGACGCCGGAAGAGCAGGAACGCGCTTTTGTTGCCCAGCTACGCATTGCCGCAGATTTAAACATGCCGGTATTTATGCACTGTCGCGATGCCCACGAGCGGTTTATGACATTGCTGGAGCCGTGGCTGGATAAACTGCC  >  NC_000913/4021403‑4021582

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